Calculo numerico (rotina Matlab)

% Este programa soluciona um problema de valor inicial em EDOS de 1°ordem.
% Não leva em conta a sensibilidade da função para variações nos valores de
% x, através de dx.
%A precisão (Formato) deve ser determinada antes da execução do programa.

%--------------------------------------------------------------------------
%Entrada de dados
%--------------------------------------------------------------------------

% entrada de dados fixa e projeto de movel, incluindo equação substituível
%x=0 %atribuição necessária para que a expressão seguinte à ser inserida não
%provoque erro por usar variaveis não inicializadas

% DÚVIDA: é necessário calcular y0 apartir de f(x,y)=y' ?
% não ! já é dado como condição inicial. disp(' ') disp('----------------------------------------------------------------------------------------------------------') disp('Masoller Corp.') disp('----------------------------------------------------------------------------------------------------------') disp(' ') disp('****Método geral para o Cálculo de EDO´s de 1° ordem****') disp(' ') disp('(Apresentação)') disp('----------------------------------------------------------------------------------------------------------') disp('-Esta rotina soluciona um problema de valor inicial em EDOS de 1°ordem;') disp('-Não leva em consideração a sensibilidade da função para valores de x através de dx;') disp('-A precisão (formato) deve ser determinada previamente à execução do programa.') disp(' ') clear all %para apagar os valores das variáveis. Aquí é a melhor opção de ficar em um programa. format long
F=input('Digite a função na forma f(x,y)=y'', na forma de operações elemento por elemento:','s');
%F='-y'
%F='x.^2-x.*y'; valor fixo p testes.

x(1)=input('entre com valor de xo:')
%x(1)=1

xf=input('entre com valor de xf:')
%xf=1.5 (valor de testes)

y(1)=input('entre com valor de yo:')
%y(1)=2

dx=input('Digite o valor inicial de dx:')
%dx=0.001