Runga Kutta Method of Order 4 Ks
Runga Kutta Method of Order 4 Ks
#include<stdio.h>
#include<math.h>
#define F(x,y) (y-x)
main()
{
int i,n;
float x0,y0,k1,k2,k3,k4,h,xn,x,y,k;
printf("give the values of x0,y0,xn and n");
scanf("%f,%f,%f,%d",&x0,&y0,&xn,&n);
h=(xn-x0)/n;
x=x0;
y=y0;
for(i=0;i<=n;i++)
{
k1=h*F(x,y);
k2=h*F(x+(h/2.0),y+(k1/2.0));
k3=h*F(x+(h/2.0),y+(k2/2.0));
k4=h*F(x+h,y+k3);
k=(1/6.0)*(k1+2*k2+2*k3+k4);
printf("\n x=%f \t y=%f",x,y);
x=x+h;
y=y+k;
}
}
KALPATARU SAHOO
PHYS24-025