function [ x, y ] = odeEULER_verb(DGL,a,b,h,y0)

N = floor((b-a)/h); x(1) = a; y(1) = y0;

for i = 1:N
    x(i+1) = x(i) + h;
	m1 = DGL(x(i), y(i));
    y(i+1) = y(i) + m1*h;

    m2 = DGL(x(i+1), y(i+1));
    y(i+1) = y(i) + (m1+m2) / 2 *h;
end 
end