Interpolacja













Interpolacja
metodą Lagrange'a 
#pragma hdrstop#include
<condefs.h>#include <conio.h>#include
<iostream.h>#include <math.h>#pragma argsused
const int n=3; // stopien
wielomianu interpolujacego
// wartosci funkcji
(y[i]=f(x[i]))
double
x[n+1]={3.0,    5.0,   
6.0,    7.0};double y[n+1]={1.732,  2.236, 
2.449,  2.646};
double interpol(double z, double
x[n], double y[n])// zwraca wartosc funkcji w punkcie
'z'{double wnz=0,om=1,w;for(int i=0;i<=n;i++) 
{  om=om*(z-x[i]);  w=1.0;  for(int
j=0;j<=n;j++)  if(i!=j) w=w*(x[i]-x[j]); 
wnz=wnz+y[i]/(w*(z-x[i]));  }return wnz=wnz*om;}
int main(int argc, char
**argv){  cout<<"Interpolacja metoda
Lagrange'a";  cout << endl;  cout <<
endl;double z=4.5;cout << "Wartosc funkcji sqrt(x) w punkcie
" << z << " wynosi " << interpol(z,x,y)
<<endl;  getch();  return
0;}Powrót