|
.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
|
|
|
e 5.3.8.a
y' = (y/t)-(y/t)^2 , for 1<=t<=2 with y(1)=1 and h=0.1
*/
#include<stdio.h>
float f(float w, float t);
main(){
int i, N = 10;
float t, h=0.1;
float w[11];
w[0] = 1;
for(i=0 ; i<N ; i++){
t = h*i +1;
w[i+1] = w[i]+ h*f(w[i]+1/2*h*f(w[i],t),t+1/2*h);
}
for(i=0 ; i<=N ; i++)
|
|
|
오일러법
2) 개선된 오일러법
3) Runge Kutta법
◆ 이상적인 현가장치
1) 완전한 노면점착(perfect road holding) - 딱딱한 스프링 시스템
2) 승차자의 안락감(passenge comfort) - 보다 유연한 스프링 시스템
◆ 실제차량설계시의 스프링 정수
|
|
|
printf("w(%1.1f) = %lf y(%1.1f) = %lf error = %lf\n", t, w[i], t, y(t), y(t)-w[i]);
}
}
double f(double t, double w){
return (w/t)-(w/t)*(w/t);
}
double y(double t){
return t/(1 + log(t));
}
Result 없음
|
|
|
문제의 해를 구하기 위해 사용될 수 있다.
3. 비교 및 토의
Euler’s Method는 실질적으로 사용하지는 않지만, 모든방법의 원리를 제공해 준다
Euler’s Method는 오차가 비교적 크다
Runge & Kutta는 정확하고, 도함수를 구할 필요가 없다
테일러 방법은
|
|
메뉴펼치기
