Trabalho de calculo numerico

Universidade Federal de Alagoas – UFAL
Campus do Sertão
Cursos de Engenharia Civil e de Produção

Disciplina: Cálculo Numérico
Professor: Márcio André Araújo Cavalcante

Lista de Exercícios 02: Zeros de Funções

Aluno: Karmem Julya Miron Batista

Delmiro Gouveia - Alagoas
Fevereiro de 2014

a) function f=mecanismo(theta) % o valor zero da função será encontrado a partir da equação : f=(cos(theta))^3-sin(theta);

b) function f=equilibrio(theta) % valores adotados: P=2.5; L=3.5; k=10; % valor da função: f=P*L*cos(theta)-k*(theta);

c) function f=equacao_de_manning(y) % valores adotados: Qmax=24; s=0.0009; n=0.012; h=1.9; b=2.5; B=6.3;
%
% área molhada:
A=(B-b)/(2*h)*y^2+b*y;
% % perímetro molhado: p=b+2*sqrt(1+((B-b)/(2*h))^2)*y; % %raio da área molhada: R=A/p; % % valor da função: f=Qmax-A*R^(2/3)*p^(1/2)/2;

MÉTODO DA BISSEÇÃO function [z,nite]=mne_bissecao(fun,a,b,tolf,nitemax)
% Dados de Entrada: fun = função;
% [a,b] = extremos do intervalo;
% tolf = tolerância;
% nitemax => número de iterações; %
% Dados de saída: z = zero estimado;
% nite = número de iterações. %
% Objetivo: Cálculo do zero de funções utilizando o método da bisseção. %
% avaliação da função para os extremos do intervalo dado:
Fa=feval(fun,a);
Fb=feval(fun,b);

% avaliação da expressão: if Fa*Fb>0
% se a avalição for verdadeira z=0 e nite=0 disp('refefina o intervalo!'); z=[ ]; nite=[ ]; end %
% fómula do método da bisseção para encontrar