Resolver o PVI pelo método de Euler Modificado usando o Octave.

• Download Euler_Modificado.m
• Enviar nova versão

Esconder código-fonte

function [vetx,vety] = Euler_Modificado(funcao,a,b,m,y0);

%parametros de entrada: funcao,a,b,m,y0 -> funcao,limites inferior, superior, num. de
%intervalos e valor inicial
%parametros de saida: vetx, vety -> abcissas e solucao do PVI

h = (b-a)/m;
x=a;
y = y0;
Fxy = eval(funcao);
vetx(1) = x;
vety(1) = y;
disp('     i      x     y     Fxy  ');
for i=1:m
    yi = y;
    x = x+h/2;    
    y = y+(h/2)*Fxy;
    Fxy = eval(funcao);
    y = yi + h*Fxy; 
   
    disp([i  x    y    Fxy]);
    vetx(i+1) = x;
    vety(i+1) = y;
end
end

Octave - Método de Euler para solução de EDO

Octave - Método de Runge-Kutta

Octave - Sistemas lineares por Gauss-Seidel

Octave - Calcular raiz pelo método de Newton

Octave - Método de Euler-Melhorado para solução de EDO

Nenhum comentário foi encontrado.