C语言程序,根据具体题型可拓展。
#include
#include
double fxy(double xi,double yi) /*定义函数fxy*/
{
double y;
y=yi-2*xi/yi;
return(y);
}
void main()
{
double x0,y0,h,xi,yi,yi_1,xk2,yk2,xk3,yk3,xk4,yk4,k1,k2,k3,k4;
int i;
x0=0; /*赋初始值*/ y0=1;
h=0.1;
xi=x0;
yi=y0;
for(i=1;i<=10;i++) /*循环开始*/
{
k1=h*fxy(xi,yi); /*求解k1值*/
xk2=xi+0.5*h; /*求解k2的值*/ yk2=yi+0.5*k1;
k2=h*fxy(xk2,yk2);
xk3=xi+0.5*h; /*求解k3的值*/ yk3=yi+0.5*k2; k3=h*fxy(xk3,yk3); xk4=xi+h; /*求解k4的值*/ yk4=yi+k3; k4=h*fxy(xk4,yk4);
yi_1=yi+(k1+2*k2+2*k3+k4)/6; /*求解yi+1的值*/ yi=yi_1;
if(i==1)
{
printf("输出函数yi的近似值:\n"); /*输出所有的yi值*/ printf("y0 = %.10f ",y0);
}