RE:Newton-Raphson complejo
Publicado por
WAWAWWA (1 intervención) el 30/09/2009 00:45:58
puto
fprintf('Iteracion de Punto Fijo\n de la funcion f(x)=x3+4x2+10\n');
OP=input('Si vas a dar un intervalo escribe "1". Si vas a dar un valor escribe "2" ');
display('intervalo de forma [a,b], ');
a=input('dame a ');
b=input('dame b ');
E=input('dame el margen de error ');
v1=(a+b)/2;
Eo=1001;
pn1=v1;
n=0;
cont=0;
fprintf('\n\ng1(x)=((10-x3)^.5)/2\n');
fprintf('n\t\tpn\t\t\t\t\tE\n');
while((Eo>E)&&(cont<4)&&(pn1<b)&&(pn1>a))
if(n==0)
fprintf('%.f\t\t%.10f\t\t---------\n',n,v1);
n=n+1;
pn1_1=v1;
Eoant=1001;
else
pn1=((10-(pn1_1)^3)^0.5)/2;
Eo=abs((pn1-pn1_1)/pn1);
fprintf('%.f\t\t%.10f\t\t%.10f\n',n,pn1,Eo);
n=n+1;
if(Eoant<=Eo)
cont=cont+1;
end;
Eoant=Eo;
pn1_1=pn1;
end;
end;
if(cont==4)
fprintf('Es divergente, el error aumenta\n\n');
ST1=0;
end;
if((pn1<a)||(pn1>b))
fprintf('Fuera del intervalo\n\n');
end;
if((pn1<a)||(pn1>b))
end;
fprintf('\n\ng2(x)=((10/x)^.5)\n');
fprintf('n\t\tpn\t\t\t\t\tE\n');
n=0;
cont=0;
while((Eo<E)&&(cont<4)&&(pn1<b)&&(pn1>a))
if(n==0)
fprintf('%.f\t\t%.10f\t\t---------\n',n,v1);
n=n+1;
pn1_1=v1;
Eoant=1001;
else
pn1=(10/(pn1_1))^0.5;
Eo=abs((pn1-pn1_1)/pn1);
fprintf('%.f\t\t%.10f\t\t%.10f\n',n,pn1,Eo);
n=n+1;
if(Eoant<=Eo)
cont=cont+1;
end;
Eoant=Eo;
pn1_1=pn1;
end;
end;
if(cont==4)
fprintf('Es divergente, el error aumenta\n\n');
ST1=0;
end;
if((pn1<a)||(pn1>b))
fprintf('Fuera del intervalo\n\n');
end;
if((pn1<a)||(pn1>b))
end;