유한요소 리포트



유한요소 리포트
Gaussian Rule 예 의한 적분
제출일자 :
제 출 자 :

1. Wright Fortran or C program for integration using n = 1, 2, 3, 4 of Gaussian Rule.

2. Calculate IEX & derive local integration form.

3. Output file.

4. Table for each I. (Error = IEX - INUM)

n
IEX
INUM
Error
1

2

3

4



1. C program source 는 아래와 같다.

#include
#include
#define pi M_PI
///// Abscissas and weight factors for Gaussian interation
///// Both matrix psi and weight are 4 by 4
double psi[4][4] =
{ 0., 0., 0., 0.,
-0.577350269189626, 0.577350269189626, 0., 0.,
-0.774596669241483, 0., 0.774596669241483, 0.,
-0.861136311594053, -0.339981043584856, 0.339981043584856, 0.861136311594053 };
double weight[4][4] =
{ 2., 0., 0., 0.,
1., 1., 0., 0.,
....

[hwp/pdf]유한요소 리포트