Implement Newton’s forward interpolation polynomial (interpolation with equidistance points)

#include

#include

void main()

{

            clrscr();

            int i,j,n,k;

            float f[30],x[30],df[30][30],xx,a,b,mul,ad;

            printf ("\t\tNewton's Forward Interpolation\n");

            printf ("\t\t------------------------------\n\n");

            printf("How many points???");

            scanf("%d",&n);

            for(i=0;i

            {

                        printf("\nEnter the value of x[%d]:",i);

                        scanf("%f",&x[i]);

                        printf("Enter the value of f[x%d]:",i);

                        scanf("%f",&f[i]);

                        df[0][i]=f[i];

            }

            for(j=1;j

            {

                        for(i=0;i

                        {

                                    df[j][i]=df[j-1][i+1]-df[j-1][i];

                        }

            }

            printf("\n Enter the value x:");

            scanf("%f",&xx);

            a=x[1]-x[0];

            b=(xx-x[0])/a;

            mul=1.0;

            ad=f[0];

            for(i=1;i

            {

                        j=i-1;

                        mul=mul*(b-j)/(float)i;

                        ad=ad+(mul*df[i][0]);

            }

            printf("\n The required value = %4f",ad);

            getch();

}