function [ x, y ] = odeEULER_CAUCHY( DGL, a,b,h,y0 ) % x, y: Vektoren mit N+1 Stützstellen, h: Schrittweite N = floor((b-a)/h); x(1) = a; y(1) = y0; for i = 1:N x(i+1) = x(i) + h; y(i+1) = y(i) + DGL(x(i), y(i))*h; end end