|
.1f) = %lf u2(%1.1f) = %lf error = %lf\\n\", t, w2[i], t, u2(t), u2(t)-w2[i]);
}
}
double f1(double t, double w1, double w2){
return 3*w1 + 2*w2 - (2*t*t+1)*exp(2*t);
}
double f2(double t, double w1, double w2){
return 4*w1 + w2 + (t*t + 2*t - 4)*exp(2*t);
}
double u1(double t){
return exp(5*t)/3 -
|
|
|
수치해석 과제
1.Euler법
2.Heun법
3.Runge-Kutta법
4.소스코드
Do it. 임의의 미분방정식 y\' = f(t, y)를 정하고, 아래 방식을 사용하여 구간 0 에서 4 까지 간격 1로 설정하여 적분하시오.
#11-1
▣ Euler 방식
접근 1.
접근 2.
#11-2
▣ Heun 방식
접근 1.
접근
|
|
|
매틀랩 M_FILE
Matlab Curvefitting Gragh
*Reference
Matlab의 입문과 활용, 김용수저, 높이깊이 -Runge-Kutta C언어 소스
-결과 출력값
-매틀랩으로 그래프 출력
-선형회귀분석 C언어 소스
-결과 출력값
-매틀랩으로 그래프 출력
|
|
|
/ Runge-Kutta 적분기
void rk4(double* y, double x,
void (*f)(double* dy, double* y, double x),
double h, int n);
// 운동방정식에 대한 연립 1차 미분방정식
void derv4(double* dz, double* z,double t);
#define DIM 4
void main()
{
int n;
double *y, h, t;
FILE *fp;
fp=fopen(\"problem8_2.dat\", \"w\");
/
|
|
|
Kutta-Joukowsky (쿠타-쥬콥스키)의 개념
- 마그너스 효과
? 운동량 특성에 의한 양력
4. 와류의 종류와 유도속도, 내리흐름--------------------------------------5
5. 유도항력이란------------------------------------------------------6
Ⅲ 결론----
|
|
메뉴펼치기
