|
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++)
|
|
메뉴펼치기
