유한요소 리포트
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.,
.... |
2011.03.11
5페이지
