C
A
R
L
O
S
ESTE PROGRAMA ESTA DESARROLLADO EN MATLAB
-------------------------------------------------------------------------------
function E_modificado(f,t0,tf,h,y0)
disp('*******************************************************************')
disp('* METODO DE EULER MODIFICADO (METODO DEL PUNTO MEDIO) *')
disp('*******************************************************************')
% Metodo del punto medio para resolver una ecuacion diferencial
% ordinaria de la forma:
%
% dy
% -- = f(x,y) con condicion incial y(t0) = y0
% dx
x = t0;
y = y0;
Y = y';
X = x;
while x<tf
yh = y + h*feval(f,x,y);
y = y + h*(feval(f,x+h/2,y+h*yh/2));
x = x + h;
Y = [Y;y'];
X = [X;x];
end
disp(' X Ye.modificado')
[X Y]
plot(X,Y,'Y'),title(' METODO DEL PUNTO MEDIO ','color','B') <>
-------------------------------------------------------------------------------
function E_modificado(f,t0,tf,h,y0)
disp('*******************************************************************')
disp('* METODO DE EULER MODIFICADO (METODO DEL PUNTO MEDIO) *')
disp('*******************************************************************')
% Metodo del punto medio para resolver una ecuacion diferencial
% ordinaria de la forma:
%
% dy
% -- = f(x,y) con condicion incial y(t0) = y0
% dx
x = t0;
y = y0;
Y = y';
X = x;
while x<tf
yh = y + h*feval(f,x,y);
y = y + h*(feval(f,x+h/2,y+h*yh/2));
x = x + h;
Y = [Y;y'];
X = [X;x];
end
disp(' X Ye.modificado')
[X Y]
plot(X,Y,'Y'),title(' METODO DEL PUNTO MEDIO ','color','B') <>
|
|
|
