Newton 1

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

double l(int n,
            double x,
            double h,
            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];
                pr=pr*(p-(x+(i-1)*h))/h/i;
                su=su+y[0]*pr;
            }
        return su;
    }

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

        printf("Dati gradul polinomului: ");
        scanf("%d",&m);
        printf("Dati pasul: ");
        scanf("%lf",&hh);
        printf("Dati abscisa: ");
        scanf("%lf",&x0);
        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,x0,hh,y,p);

        printf("Rezultatul este: %lf",r);

        getche();
        return 1;
    }