附下我自己编的 但是运行结果有错 求改正！

#include<stdio.h>
#include<math.h>
#define f(x,y,z) x*(1-0.1*y)
#define g(x,y,z) y*(-0.5+0.02*x) void main(void)
{
double x,y,z,x0,y0,z0,k1,k2,k3,k4,L1,L2,L3,L4,t,h;
printf("Enter x0,y0,h:");
scanf("%lf%lf%lf",&x0,&y0,&h);
x=x0;
y=y0;
for(t=0;t<=25;t=t+h){
k1=f(x,y,z);
L1=f(x,y,z);
k2=f(x+h/2,y+k1*h/2,z+L1*h/2);
L2=f(x+h/2,y+k1*h/2,z+L1*h/2);
k3=f(x+h/2,y+k2*h/2,z+L2*h/2);
L3=f(x+h/2,y+k2*h/2,z+L2*h/2);
k4=f(x+h/2,y+k3*h,z+L3*h);
L4=f(x+h/2,y+k3*h,z+L3*h);
printf("%lf\t%lf\t",x,y);
y=y+h*(k1+2*k2+2*k3+k4)/6;
x=z+h*(k1+2*k2+2*k3+k4)/6;
}
printf("t=%lf\t xn=%lf\t yn=%lf\n",t,x,y);
}