初值问题

初值问题

-

2023年3月16日发(作者：进项发票)

#include""

#include""

#include""

#defineN10

floatf(floatx,floaty)

{

floats;

s=y-2*x/y;

returns;

}

main()

{

floath,T1,T2,x=,y=,a=,b=;

inti;

h=(b-a)/N;

printf("OUTPUT:N");

for(i=0;i<=N-1;i++)

{

printf("k=%dtx%d=%fty%d=%fn",i,i,x,i,y);

T1=y+h*f(x,y);

x+=h;

T2=y+h*f(x,T1);

y=(T1+T2)/2;

}

printf("k=%dtx%d=%fty%d=%fn",i,i,x,i,y);

}

求解方程为













1)0(

)10(/2

y

xyxyy

的数值解。其中区间为N=10份