用四阶龙格库塔法求解初值问题

用四阶龙格库塔法求解初值问题
如图 取h=0.1算到x=1.5，并与精确解y=1/（x-2）相比较，保留4位小数



#include <iostream.h>
#include <iomanip.h>
#include <math.h>
double f(double );
void main()
{double k1,k2,k3,k4,y[5]={-1},Y[6]={-1},x=1,h=0.1,RY[6],z[5]={-1},L1,L2,L3,L4;
for(int i=0;i<4;i++)
{k1=z[i];
L1=f(y[i]);
k2=z[i]+h/2*L1;
L2=f(y[i]+h/2*k1);
k3=z[i]+h/2*L2;
L3=f(y[i]+h/2*k2);
k4=z[i]+h*L3;
L4=f(y[i]+h*k3);
y[i+1]=y[i]+h/6*(k1+2*k2+2*k3+k4);
z[i+1]=z[i]+h/6*(L1+2*L2+2*L3+L4);
}

for(i=0;i<5;i++)
Y[i+1]=Y[i]+h*y[i];
for(i=0;i<6;i++)
{RY[i]=1/(x-2);
x=x+0.1;}
cout<<"y        "<<"RY        "<<"误差"<<endl;