% Método de Euler
% y'(x)=f(x,y) para x en el intervalo (x0,xn)
% y(x0)=y0

% Condición inicial
x0=0.0
y0=3.0

% Último nodo, número de nodos, paso. 
xn = 1
n = 10
h = (xn - x0)/n

% Función f
 f=@(x,y) 2*x*y

% Vector de nodos. Inicializamos Y
X=[x0:h:xn]
Y=zeros(1,n+1)

% Y(1), Y(2)
Y(1)=y0
Y(2)=Y(1) + f(X(1), Y(1)) * (X(2) - X(1))

% Bucle principal
for i= 3:n+1
Y(i)=Y(i-1) + f(X(i-1), Y(i-1)) * (X(i) - X(i-1));
end

% Resultados
[X' , Y']


% Gráfica de la poligonal
plot(X,Y)
