Giải gần đúng phương trình vi phân thường cấp I(PP Euler)

//Giai phuong trinh vi phan-PP Euler
#include <stdio.h>
#include <conio.h>

float f(float x,float y)
{
return (y-2*x)/(2*y+x);
}

void main()
{
float a,b,y0;
int n;
float h;
int i;
float x***91;100***93;,y***91;100***93;;
clrscr();
printf("Giai phuong trinh vi phan bang phuong phap Euler");
printf("\nPtr: y'=(y-2x)/(2y+x)");
printf("\nNhap so lieu:");

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:");

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

printf("\nx0=%f",a);
printf("\ty0=%f",y0);
printf("\n");

for (i=1;i<=n;i++)
{
x***91;i***93;=x***91;i-1***93;+h;
y***91;i***93;=y***91;i-1***93;+h*f(x***91;i-1***93;,y***91;i-1***93;);

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