Matlab - ayuda con codigo matlab

 
Vista:

ayuda con codigo matlab

Publicado por Juan Guaman (1 intervención) el 06/10/2016 03:17:38
1
2
3
4
5
6
7
8
9
10
11
12
function y=newton_r(fun,x0,tol)
dfun=diff(sym(fun));
f=subs(fun,x0);
df=subs(dfun,x0);
c=1;
while abs(f)>tol
    a=x0-f/df;
    f=subs(fun,a);
    df=subs(dfun,a);
    x0=a;
    c=c+1;
end
c
a
este es mi codigo en matlab para el metodo de newton pero cuando lo corro me sale que a =4/3 y quiero q me salga 1.00000
la funcion q estoy utilizando es x^2-1

ayudenme con eso por favor
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de JOSE JEREMIAS CABALLERO
Val: 6.975
Oro
Ha mantenido su posición en Matlab (en relación al último mes)
Gráfica de Matlab

ayuda con codigo matlab

Publicado por JOSE JEREMIAS CABALLERO (5917 intervenciones) el 06/10/2016 13:37:09
1
2
3
4
5
6
7
8
9
10
11
12
13
function y=newton_r(fun,x0,tol)
dfun=diff(sym(fun));
f=subs(fun,x0);
df=subs(dfun,x0);
c=1;
while abs(f)>tol
    a=x0-f/df;
    f=subs(fun,a);
    df=subs(dfun,a);
    x0=a;
    c=c+1;
end
y=double(a);

1
2
3
>> syms x; fun=x^2-1; x0=0.33; tol=0.0001; y=newton_r(fun,x0,tol)
y =
    1.0000


Saludos.
JOSE JEREMÍAS CABALLERO
Asesoría online en Matlab
Servicios de programación matlab
[email protected]
El correo es para servicios de programación, toda ayuda gratuita es vía foro.


http://matlabcaballero.blogspot.com
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar