|
수정된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
|
|
|
\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
|
|
|
.0)*w[9]+(1.97-1.9)/(2.0-1.9)*w[10];
printf("w(1.97) = %f y(1.97) = %f error = %f\n",y3,y(1.97),y3-y(1.97));
}
float f(float w, float t){
return 2*w/t+t*t*exp(t);
}
float y(float t){
return t*t*(exp(t)-exp(1));
}
Result
/*
This program is Euler's Method for exercise 5.2.10
y' = exp(-0.06*pi*t)*(1.20
|
|
|
않을까?
빌라니 설명 : 아주 간단한 3차원 구 모양의 형태로 구멍이 없고 좌표가 세 개인 그런 3차원의 구
한붓그리기를 알면 문제 이해가 좀 더 쉬움 1부. 수의시작
2부. 원론
3부. 신의 숫자
4부. 움직이는세계 미적분
5부 남겨진 문제들
|
|
|
적분하면,
∴ Q=AV
즉, 관 내(또는 유수)의 단면적과 그 점에서의 유속을 곱한 것은 항상 일정하다.
그리고 이것을 Q라 쓰고 유량(discharge)이라고 한다. 1. Euler(오일러)의 운동방정식
2. Euler 의 연속방정식
1) 일반 유동에 대한 연속방정식
|
|
메뉴펼치기
