INTERPOLACION CUADRATICA DE NEWTON

//Torres Martinez Gustavo Alan
//Grupo:4CM1
//Programa que realiza diferencias divididas para determinar el valor de
//de una funcion conocida en base al valor de la funcion evaluada en puntos
//equidistantes por medio de "INTERPOLACION CUADRATICA DE NEWTON"

#include<stdio.h>
#include<conio.h>
#include<dos.h>
#include<math.h>
#include<stdlib.h>
#include<iostream.h>

void main(void)
{
double fx,fx0,fx1,fx2,f2x,b0,b1,b2,x,x0,x1,x2,erv;
clrscr();
x=2;
x0=1;
x1=3;
x2=5;
fx=log(x);
fx0=log(x0);
fx1=log(x1);
fx2=log(x2);
cout<<"Calculando el ln(2) por medio del polinomio de interpolacion cuadratica"<<endl<<"de Newton con x=1,3,5."<<endl<<endl;
b0=fx0;
b1=(fx1-fx0)/(x1-x0);
b2=((fx2-fx1)/(x2-x1)-(fx1-fx0)/(x1-x0))/(x2-x0);
f2x=b0+b1*(x-x0)+b2*(x-x0)*(x-x1);
erv=(fx-f2x)/fx*100;
printf("bo=%1.8f  b1=%1.8f  b2=%1.8f\n\n",b0,b1,b2);
printf("f2(x)=%1.8f Erv=%2.2f%",f2x,fabs(erv));
cout<<endl<<endl;
textcolor(10+128);
cprintf("             Presiona cualquier tecla....");
getch();
}