Newton 3

#include<stdio.h>
#include<conio.h>
#include<math.h>


double l(int n,
            double x[],
            double y[],
            double p)
    {
        double su,pr;
        int i,j;
        su=y[0];
        pr=1;
        for (i=1;i<=n;i++)
            {
                for (j=0;j<=n-i;j++) y[j]=(y[j+1]-y[j])/(x[j+i]-x[j]);
                pr=pr*(p-x[i-1]);
                su=su+y[0]*pr;
            }
        return su;
    }

int main()
    {
        int i,m;
        double x[20],y[20],p,r;

        printf("Dati gradul polinomului: ");
        scanf("%d",&m);

        printf("Dati abscisele\n");
        for (i=0;i<=m;i++)
            {
                printf("x[%d]: ",i);
                scanf("%lf",&x[i]);
            }
        printf("Dati ordonatele:\n");
        for (i=0;i<=m;i++)
            {
                printf("y[%d]: ",i);
                scanf("%lf",&y[i]);
            }

        printf("Dati abscisa in care doriti calculul: ");
        scanf("%lf",&p);

        r=l(m,x,y,p);
        printf("Rezultatul este: %lf",r);

        getche();
        return 1;
    }