Runga Kutta Method of Order 4 in C

    #include<stdio.h>
    #include<math.h>
    float f(float x,float y)
    {
        return (x*x+y*y);
    }
    main()
    {
        printf("\nProgram to implement Runga Kutta Method (x-y) : \n");
        int i,n;
        float x0,y0,x,y,h,k,k1,k2,k3,k4;
        printf("\nEnter the values of x0,y0,h,x : ");
        scanf("%f%f%f%f",&x0,&y0,&h,&x);
        n=(x-x0)/h;
        x=x0,y=y0;
        for(i=0;i<=n;i++)
        {
            k1=h*f(x,y);
            k2=h*f(x+h/2,y+k1/2);
            k3=h*f(x+h/2,y+k2/2);
            k4=h*f(x+h,y+k3);
            k=(k1+(k2+k3)*2+k4)/6;
            x=x+h;y=y+k;
        }
        printf("\n\nThe required value is %f\n\n",y);
    }

Top