Matlab Método de Euler en Matlab.

[lindtaylor]

Buenas, quiero resolver en Matlab (en realidad era Scilab pero luego intentaré moldear el código de Matlab a Scilab) mediante las poligonales de Euler el problema de cauchy siguiente:

$\left\{ \begin{array}{lcc}<br /> x'=-x & \\<br /> \\ x(0)=1 & \\<br /> <br /> \end{array}<br /> \right.$

Sobre el intervalo $[0,1]$ con tamaño de paso $h=0.05.$

¿Cómo sería el código? Soy nuevo en Matlab .

Desde ya muchas gracias.

[[Beto]]

Re: Método de Euler en Matlab.

hola, esto te puede servir http://mate.uprh.edu/~pnm/notas4061/odes1/odes1.htm

[lindtaylor]

Re: Método de Euler en Matlab.

Probé los códigos y no me funcionan :C. O tal vez sea porque no se usarlo bien o quizás sea porque tengo Matlab 5.6 .

¿Alguien podría poner el código para llegar y pegarlo y probar si es un problema distinto el que no me funcione? .

[lindtaylor]

Re: Método de Euler en Matlab.

function [Y, X] = Euler(f, x0, y0, xf, n)
// Metodo de Euler para la ecuacion diferencial
//
// y’ = f(x,y)
// y(x0) = y0
// en intervalo [x0, xf]
//
// n = numero de subintervalos
// Y, X seran vectores fila de n+1 elementos
// Y contendra las aproximaciones de
// y(x0) y(x0+h) y(x0+2h) ... y(xf)
// con h = (xf-x0)/n
// X contendra los valores x0 x0+h x0+2h ... xf
h = (xf-x0)/n
X = zeros(1,n+1)
Y = X
X(1) = x0
Y(1) = y0
xi = x0
yi = y0
for i=1:n
yi = yi + h*f(xi,yi)
xi = xi+h
Y(i+1) = yi
X(i+1) = xi
end
endfunction


Este código es para Scilab, alguien podría adaptarlo para mi problema .