euler_method.sci
clc clear //metoda forward Euler dydx = 1-x*y, y(0) = 1 x= [0, 3] L = 3 h = 0.1 x = 0:h:L; y0 = 1 function yp=f(x,y) yp = 1-x*y endfunction n = length(x) y = zeros(1,n) y(1) = y0 //warunek poczatkowy for i=1:n-1 y(i+1) = y(i) + h*f(x(i), y(i)); end plot(x,y, 'bo') //rozwiazanie referencyjne x = 0:0.1:3; y = ode(y0, 0,x,f) plot(x,y,'r')