#include<stdio.h>
#include<math.h>
float f(float x, float y)
{
return (x+y);
}
void main()
{
float x0,y0,h,x,m,m1,m2,m3,m4;
printf("Enter x0, y0, h , x : ");
scanf("%f%f%f%f",&x0,&y0,&h,&x);
while(x0<x)
{
m1=f(x0,y0);
m2=f(x0+h/2,y0+m1*h/2);
m3=f(x0+h/2,y0+m3*h/2);
m4=f(x0+h,y0+m3*h);
m=(m1+2*m2+2*m3+m4)/6;
y0+=h*m;
x0+=h;
printf("x = %2.4f \t\t y = %2.4f\n",x0,y0);
}
}
0 comments:
Post a Comment