|
수정된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));
}&nb
|
|
|
적분학, 경문사, 1996
4. 송만석, 장건수 편역, 수치해석학, 김영사, 1986
5. 정상권, 하성남, 수치해석, 경문사, 1997
6. D.Conte, C.de Boor, Elementary Numerical Analysis, McGraw-Hill,Inc. 1980
7. Eusebius Doedel, Numerical Analysis, Lecture Notes, 1999
8. N. Blachman, Mathematica : A Pra
|
|
|
\n",j-1,X[i],j-1,e[i+1],j-1,E[i],j-1,DEF[i]);
j++;
i++;
}
}
//N=100 일때
#include <stdio.h>
#include <math.h>
#define e Euler_Y_Value
#define E EXACT_Euler_Value
#define X Euler_X_Value
main()
{
int i=0,j=1;
float X[101],x=1,y=1,h=0.01,e[102],E[102],c[102],DEF[102];//초기값 X=1,y=1 h=0.0
|
|
|
수치해석이란 과목을 들으면서, 어려운 부분도 많이 있었지만, 그 만큼 얻어가는 지식도 많은 것 같습니다. 너무 수고하셨습니다.^_^ 1. 서론
- 과제내용
- 문제1
- 문제2
2. 본론
- 문제1 풀이
- 문제2 풀이
3. 결론
- 결과
|
|
|
해석학에서는 순서, 규칙성이 조사된 반면, 이 시기에는 예외와 불규칙이 조사되었다. 바이어쉬트라스(Weierstrass)는 1872년에 모든 곳에서 미분 불가능한 연속함수를 제시하였다. 그는 그럴 듯 하고 널리 인정된 개념에 반례를 제시함으로써 엄
|
|
메뉴펼치기
