euler.m

% Lab 5.
% Marta Trzaska 171632 SiSR

 function [u_C, i_L] =euler(x_0, t, h, C, L, E, R)
    t_h = t/h
    u_C = zeros(1, t_h + 1);
    i_L = zeros(1, t_h + 1);

    sprawdzajka_uC = zeros(1, t_h + 1); %
    sprawdzajka_iL = zeros(1, t_h + 1); %

    u_C(1) = x_0(1);
    i_L(1) = x_0(2);

    for i=2:1:(length(t_h)+1)
        [f_uC, f_iL] = funkcja(C, L, E, R, i_L(i - 1), u_C(i - 1));

        sprawdzajka_uC(i) = f_uC;
        sprawdzajka_iL(i) = f_iL;

        i_L(i) = i_L(i-1) + h*f_iL;
        u_C(i) = u_C(i-1) + h*f_uC;
    end

 end