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

• Download Euler_Melhorado.m
• Enviar nova versão

Esconder código-fonte

function [vetx,vety] = Euler_Melhorado(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

    xa = x;

    ya = y;

    x = x+h;

    y = y+h*Fxy;

    ter = eval(funcao);

    y = ya + (h/2)*(Fxy+ter);

    x= a+i*h;

    Fxy = eval(funcao);    

    disp([i  x    y    Fxy]);

    vetx(i+1) = x;

    vety(i+1) = y;

end 

end



    

Octave - Sistemas lineares por Gauss-Seidel

Octave - Método de Runge-Kutta

Octave - Método de Gauss-Legendre

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

Cálculo da raiz de uma função pelo Método da Bisseção - Octave

Nenhum comentário foi encontrado.