C:\MATLAB6p1\work\pol.m
May 27, 2004
Page 1
12:02:26 AM
clear
clc
disp('*********************************************************************
**************')
disp(' METODO NUMERICO')
disp(' FACULDADE DE ENGENHARIA QUIMICA DE LORENA')
disp(' PROF. OSWALDO LUIZ COBRA GUIMARAES')
disp(' email [email protected]')
disp('*********************************************************************
**************')
n=input('entre com o grau do polinomio ');
a=input('entre com o extremo inferior ');
b=input('entre com o extremo superior ');
delta=(b-a)/n;
t=a:((delta)/100):b;
y=1./(1+25*t.^2);
x(1)=a;
for i=2:n+1
x(i)=a+(i-1)*delta;
end
for i=1:n+1
a(i)=1;
end
for i=1:n+1
for j=1:n+1
a(i,j)=x(i)^(j-1);
end
end
hold on
for i=1:n+1
f(i)=1./(1+25*x(i).^2);
end
c=a\f';
d=c;
j=n+1;
for i=1:n+1
c(i)=d(j);
j=j-1;
end
hold off
z=polyval(c,t);
comet(t,y)
hold on
comet (t,z)
plot(t,y,'b',t,z,'g')
legend('curva real','aproximada');
title('interpolaçao polinomial' );