C Program for Heun's Method

Program #include<conio.h> #include<stdio.h> #define f(x,y) 2*y/x void main() { float x,y,h,xn,l; clrscr(); printf(...

Program

#include<conio.h>
#include<stdio.h>
#define f(x,y) 2*y/x
void main()
{
float x,y,h,xn,l;
clrscr();
printf("Program for Solution of Ordinary Differential Equation\nHeun's Method\n");
printf("Enter value for x and y\n");
scanf("%f%f",&x,&y);
printf("Enter value for h and last of x\n");
scanf("%f%f",&h,&xn);
while(x+h<=xn)
{
l=(h/2)*(f(x,y)+f(x+h,y+h*f(x,y)));
y=y+l;
x=x+h;

printf("y = %f\tx = %f\n",y,x);
}
 getch();
}

Output