Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #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[n];
- pr=1;
- for (i=1;i<=n;i++)
- {
- for (j=n;j>=i;j--) y[j]=y[j]-y[j-1];
- pr=pr*(p-(x-(i-1)*h))/(h*i);
- su=su+y[n]*pr;
- }
- return su;
- }
- int main()
- {
- int i,m;
- double xn,y[20],p,r,hh;
- printf("Dati gradul polinomului: ");
- scanf("%d",&m);
- printf("Dati pasul: ");
- scanf("%lf",&hh);
- printf("Dati ultima abscisa: ");
- scanf("%lf",&xn);
- 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,xn,hh,y,p);
- printf("Rezultatul este: %lf",r);
- getche();
- return 1;
- }
Add Comment
Please, Sign In to add comment