14. Write a program to demonstrate Runge-Kutta fourth order Method.

        /* Program to execute Range-Kutta method of fourth order */
	#include<stdio.h>
	#include<conio.h>
	#include<math.h>
	float dif ( float x,float y )
	{
	    return   x+y*y;
	}
	void main()
	{
	    float x  ,h,x0,x1,y0,y1,k1,k2,k3,k4,k;
	    printf("Enter the values of x, h, x0, x0, y0: ");
	    scanf("%f%f%f%f", &x, &h, &x0,&y0);
	    x1=x0;
	    y1=y0;
	    while(x1 < x) {
		k1=h*dif(x1,y1);
		printf("\n\tk1=%f",k1);
		k2=h*dif(x1+h/2,y1+k1/2);
		printf("\n\tk2=%f",k2);
		k3=h*dif(x1+h/2,y1+k2/2);
		printf("\n\tk3=%f",k3);
		k4=h*dif(x1+h,y1+k3);
		printf("\n\tk4=%f",k4);
		k=(k1+2*(k2+k3)+k4)/6;
		printf("\n\tk=%f",k);
		y1 += k;
		x1 +=h;
		printf("\n Values of y=%f, when x=%f",y1,x1);
	    }
	    getch();
	}