Carlos André Vaz Junior
[email protected]
http://www.eq.ufrj.br/links/h2cin/carlosandre
EQ/UFRJ
fsolve
EQ/UFRJ
principal.m
solucao = fsolve('funcao',[2],optimset('Display','iter'))
function [y] = funcao(x)
y=2*x+1;
EQ/UFRJ
principal.m
solucao = fsolve('funcao',[2],optimset('Display','iter'))
function [y] = funcao(x)
y=2*x^2+x-2;
EQ/UFRJ
solucao = fsolve('funcao',[-1],optimset('Display','iter'))
function [y] = funcao(x)
y=8*x^2+4*x-2;
Experimente outro chute inicial:
solucao = fsolve('funcao',[5],optimset('Display','iter'))
EQ/UFRJ
close all
A=8;
B=4;
function [y] = funcao(x,A,B,C)
y=A*x^2+B*x+C;
C=-2;
figure(1)
x=-2:0.01:2;
plot(x,y,'*k')
y=A.*x.^2+B.*x+C;
drawnow
figure(1)
pause(0.1)
plot(x,y)
hold on
solucao = fsolve('funcao',[2],optimset('Display','iter'),A,B,C)
plot(solucao,0,'r*')
EQ/UFRJ
Achando todas as raízes do polinômio:
A=roots([8 4 -2])
EQ/UFRJ
A= 8;
B= 4;
C= -2;
solucao = fsolve('funcao',[-1],optimset('Display','iter'),A,B,C)
function [y] = funcao(x,A,B,C)
y=A*x^2+B*x+C;
EQ/UFRJ
solucao = fsolve('funcao',[-4 1],optimset('Display','iter'))
function [y] = funcao(x)
y(1)= ( x(1).^2 ) - 5 + x(2);
y(2)= (2*x(2))-2;
EQ/UFRJ