problemas con euler
Publicado por marcos (2 intervenciones) el 11/05/2009 14:21:10
Hola buenas,
estoy con un programa para el metodo de Euler y no consigo que me devuelva el valor de y
la secuencia del programa es la siguiente:
% metodo de Euler.m
function [X,Y] = dEuler(f,h,y0,x0,xf)
% Metodo de Euler para resolver una ecuacion diferencial
% ordinaria de la forma:
% dy
% -- = f(x,y) con condicion incial y(x0) = y0
% dx
x = x0;
y = y0;
Y = y';
X = x;
while x<xf
y = y + h*feval(f,x,y);
x = x + h;
Y=[Y;y'];X=[X;x];
end
la ecuación diferencial la introducco primero en un fichero al que llamo myfun, da bien el valor de f y de x pero no consigo el de y.
Muchas gracias por adelantado.
estoy con un programa para el metodo de Euler y no consigo que me devuelva el valor de y
la secuencia del programa es la siguiente:
% metodo de Euler.m
function [X,Y] = dEuler(f,h,y0,x0,xf)
% Metodo de Euler para resolver una ecuacion diferencial
% ordinaria de la forma:
% dy
% -- = f(x,y) con condicion incial y(x0) = y0
% dx
x = x0;
y = y0;
Y = y';
X = x;
while x<xf
y = y + h*feval(f,x,y);
x = x + h;
Y=[Y;y'];X=[X;x];
end
la ecuación diferencial la introducco primero en un fichero al que llamo myfun, da bien el valor de f y de x pero no consigo el de y.
Muchas gracias por adelantado.
Valora esta pregunta


0