Phương pháp Runge-Kutta

//Phuong phap Runge-Kutta
#include <stdio.h>
#include <conio.h>

float f(float x,float y)
{
return x+y;
}
void main()
{
float a,b,y0;
int n;
float h;
int i;
float x***91;100***93;,y***91;100***93;;
float k1,k2,k3,k4;

clrscr();
printf("Giai phuong trinh vi phan:");
printf("\ny'=x+y");
printf("\n");

printf("\na=");
scanf("%f",&a);

printf("\nb=");
scanf("%f",&b);

printf("\ny0=");
scanf("%f",&y0);

printf("\nn=");
scanf("%d",&n);

getch();
clrscr();

printf("Bang ket qua:");
printf("\nx0=%f",a);
printf("\ty0=%f",y0);
printf("\n");

h=(b-a)/n;
x***91;0***93;=a;
y***91;0***93;=y0;

for (i=1;i<=n;i++)
{
x***91;i***93;=x***91;i-1***93;+h;
k1=h*f(x***91;i-1***93;,y***91;i-1***93;);
k2=h*f(x***91;i-1***93;+h/2,y***91;i-1***93;+k1/2);
k3=h*f(x***91;i-1***93;+h/2,y***91;i-1***93;+k2/2);
k4=h*f(x***91;i-1***93;+h,y***91;i-1***93;+k3);
y***91;i***93;=y***91;i-1***93;+(k1+2*k2+2*k3+k4)/6;

printf("\nx%d=%f",i,x***91;i***93;);
printf("\ty%d=%f",i,y***91;i***93;);
printf("\n");
}

getch();
}