CODIGOS MATLAB AREA 1

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

area1/.DS_Store
__MACOSX/area1/._.DS_Store
area1/exemplo_metodo_grafico.m
x=[0.1:0.1:3];
f1 = @(x) log(x);
f2 = @(x) 2-x;
figure(1)
plot(x,f1(x),'r');
hold on
plot(x,f2(x),'b');
hold off
area1/exemplo_raizes.m
f = @(x) sin(x);
df = @(x) cos(x);
x1=mat_bissecao(f,3,4,5)
x2=mat_pfalsa(f,3,4,5,0,0,0)
x3=mat_secante(f,3,4,5,0,0,0)
x4=mat_newton(f,df,3,5,0,0,0)
__MACOSX/area1/._exemplo_raizes.m
��������Mac OS X �����	���2���F�������x����MATFMATL��������������������������ATTR�������x���x��������������������
area1/mat_bissecao.m
function p=mat_bissecao(f,a,b,N)
% Metodo da Bissecao
%
% p=mat_bissecao(f,a,b,N)
fprintf('Metodo da Bissecao\n\n')
fprintf('i\ta\t\tb\t\tp\n')
for i=1:N
 
 %Ponto medio
 p=(a+b)/2;
 
 fprintf('%d\t%f\t%f\t%f\n',i,a,b,p)
 
 %Criterio de parada
 if f(p)==0 
 break 
 end
 
 %Escolha de "a" e "b" da proxima iteracao
 if f(a)*f(p)<0
 b=p;
 else
 a=p;
 end 
end
__MACOSX/area1/._mat_bissecao.m
��������Mac OS X �����	���2���F�������x����MATFMATL��������������������������ATTR�������x���x��������������������
area1/mat_newton.m
function xn1=mat_newton(f,df,xn,N,e1,e2,e3)
% Metodo de Newton
%
% xn1=mat_newton(f,df,xn,N,e1,e2,e3)
fprintf('Metodo de Newton\n\n')
fprintf('i\tx(n)\t\tx(n+1)\n')
for i=1:N
 
 %Novo ponto
 xn1=xn-df(xn)\f(xn);
 
 fprintf('%d\t%f\t%f\n',i,xn,xn1)
 
 %Criterios de parada
 if isnan(xn1)
 xn1=xn;
 return;
 end
 if abs(f(xn1))<=e1 
 break 
 end
 if abs(xn1-xn)<e2 
 break 
 end
 if abs((xn1-xn)/xn1)<e3 
 break 
 end
 
 
 xn=xn1;
end
__MACOSX/area1/._mat_newton.m
��������Mac OS X �����	���2���F�������x����MATFMATL��������������������������ATTR�������x���x��������������������
area1/mat_pfalsa.m
function p=pfalsa(f,a,b,N,e1,e2,e3)
% Metodo da Posicao Falsa
%
% p=pfalsa(f,a,b,N,e1,e2,e3)
fprintf('Metodo da Posicao Falsa\n\n')
fprintf('i\ta\t\tb\t\tp\n')
for i=1:N
 
 %Novo ponto
 p=(a*f(b)-b*f(a))/(f(b)-f(a));
 fprintf('%d\t%f\t%f\t%f\n',i,a,b,p)
 %Criterios de parada
 if abs(f(p))<=e1
 break 
 end
 if i>1
 if abs(p-panterior)<e2 
 break
 end
 if abs((p-panterior)/p)<e3 
 break
 end
 end
 
 %Escolha de "a" e "b" da proxima iteracao
 if f(a)*f(p)<0
 b=p;
 else
 a=p;
 end
 
 panterior=p; 
end
__MACOSX/area1/._mat_pfalsa.m
��������Mac OS X �����	���2���F�������x����MATFMATL��������������������������ATTR�������x���x��������������������
area1/mat_secante.m
function xn1=mat_secante(f,xn,xnm1,N,e1,e2,e3)
% Metodo da Secante
%
% xn1=mat_secante(f,xn,xnm1,N,e1,e2,e3)
fprintf('Metodo da Secante\n\n')
fprintf('i\tx(n-1)\t\tx(n)\t\tx(n+1)\n')
for i=1:N
 %Novo ponto
 xn1=(xn*f(xnm1)-xnm1*f(xn))/(f(xnm1)-f(xn));
 fprintf('%d\t%f\t%f\t%f\n',i,xnm1,xn,xn1)
 
 %Criterios de parada
 if isnan(xn1)
 xn1=xn;
 return;
 end
 if abs(f(xn1))<=e1
 break 
 end
 if abs(xn1-xn)<e2 
 break 
 end
 if abs((xn1-xn)/xn1)<e3 
 break 
 end
 
 xnm1=xn;
 xn=xn1;
end
__MACOSX/area1/._mat_secante.m
��������Mac OS X �����	���2���F�������x����MATFMATL��������������������������ATTR�������x���x��������������������