grafica para interpolacion lineal
Publicado por cesc (2 intervenciones) el 19/02/2013 04:06:34
Hola, alguien sabe como puedo graficar el error en cada iteracion en funcion de las iteraciones? Tambien debo dar una grafica de cada x3 o solucion en funcion de las iteraciones...Estoy tratando y me a solo un punto. Gracias por la ayuda. Abajo esta como programe el metodo de interpolacion lineal.
function Sol=IL(x1,x2)
F=inline('x^3-71*x^2/30-61*x/30+2/5','x');
%Verificamos Bolzano
if F(x1)*F(x2)>0
display('No cumple Bolzano');
return
end
%Definimos
TolL2=10^-3;
c3=1;
TolDsv=10^-3;
Nmax=50;
eL2=1;Dsv=1;n=0;
while ((eL2>TolL2)||(Dsv>TolDsv)) && (n<Nmax)
%Definimos x3
x3=x2-(F(x2))*(x2-x1)/F(x2)-F(x1);
%Tenemos dos casos
if F(x3)*F(x2)<0
x2=x3;
end
if F(x3)*F(x1)<0
x1=x3;
end
eL2=abs((x3-c3)/x3);
c3=x3;
Dsv=F(x3);
n=n+1;
%Finalmente
Sol=x3;
end
function Sol=IL(x1,x2)
F=inline('x^3-71*x^2/30-61*x/30+2/5','x');
%Verificamos Bolzano
if F(x1)*F(x2)>0
display('No cumple Bolzano');
return
end
%Definimos
TolL2=10^-3;
c3=1;
TolDsv=10^-3;
Nmax=50;
eL2=1;Dsv=1;n=0;
while ((eL2>TolL2)||(Dsv>TolDsv)) && (n<Nmax)
%Definimos x3
x3=x2-(F(x2))*(x2-x1)/F(x2)-F(x1);
%Tenemos dos casos
if F(x3)*F(x2)<0
x2=x3;
end
if F(x3)*F(x1)<0
x1=x3;
end
eL2=abs((x3-c3)/x3);
c3=x3;
Dsv=F(x3);
n=n+1;
%Finalmente
Sol=x3;
end
Valora esta pregunta


0