수치해석학 개량된 오일러 매소드로 적분하기

소개글

수치해석학 개량된 오일러 매소드로 적분하기에 대한 보고서 자료입니다.

본문내용

hile(E[i]-y1>0.0001)//오차 10-4
{
y1=y+h*(f(c,y)+f(x,y1))/2;//개량된 Euler 구하기
}
y=y1;
}
return;
}
// f(t,y) 서브루틴
double f(double x, double y)
{
return(x*pow(y,1/3));
}
주어진 문제 ) y ' = x · y1/3 =f(x, y)
y ( 1 ) = 1
개량된 Euler 방법에 의한 Y값 추정
N=100 H=0.01 x=1 y=1 y(1)=1
#include
#include
double f(double, double); // 함수 f(x,y) 서브루틴 선언
double y;
void main()
{
int n;
double x=1.00, h=0.01, y1=y=1,E[101],X[101],c;
for(int i=0; i<101; i++)
{
c=X[i]=x;
E[i]=pow((pow(X[i],2)+2)/3,1.5); //실제값
printf("X=%.2lf 수정된Euler=%.5lf 실제적분값=%.5lf\n", x, y,E[i]);
x=x+h;
y1=y+h*f(c,y);
y=y+h*(f(c,y)+f(x,y1))/2;// y(n+1)=y(n)+h*(f[x,y]+f[x+h,y+h*f[x,y]])/2
y1=y;
while(E[i]-y1>0.0001)
{
y1=y+h*(f(c,y)+f(x,y1))/2;
}
y=y1;
}
return;
}
// f(t,y) 서브루틴
double f(double x, double y)
{
return(x*pow(y,1/3));
}

키워드

개량된 오일러,   적분법,   수치해석,   C소스,   알고리즘,   n=100과 n=20인 경우