% Metoda Runge-Kutta 4. reda - uporaba algoritma 6.4.1

N = 10;
h = 1/N;
x0 = 0;
y0 = 2;

alg641

z = 1 + exp(-X);
figure(2)
plot(X,Y-z)
        title('Metoda Runge-Kutta 4. reda')
        ylabel('globalna napaka')
        xlabel('x')
figure(1)
plot(X,Y,'r*',X,z,'g')
        title('Metoda Runge-Kutta 4. reda')
        xlabel('x')
        ylabel('Num. resitev (rdeca) in tocna resitev (zelena)')