function RLC_reseni
% Reseni serioveho obvodu R,L,C s konstantnim zdrojem napeti
% Na vetev tvorenou seriove spojenym odporem 200 Ohm, civkou 0.5H a kondenzatorem 5 mikroF pripojime 
% v okamziku t=0 zdroj napeti s konstantnim napetim U0=100V. 
% Napeti na kondenzatoru v okamziku t=0 je uC(0)=0.

R=200; % R=200 Ohm
L=0.5; % L=0,5 H
C=5e-6; % C=5 mikroF
U0=100; % U0=100V
[t,proud]=ode45(@RLC_dif_rov,[0,0.1],[0,(100/0.5)]); 
   % resi soustavu diferencialnich rovnic danou v m-file RLC_dif_rov.m
   % v casovem intervalu od 0 do 0.1 s pocatecnimi podminkami i(0)=0, di/dt(0)=U0/L=100/0.5;
   % zname fyzikalni pocatecni podminky i(0)=0, uC(0)=0,
   % z 2. Kirchhoffova zakona: uR + uL + uC =  U0 => uL = U0-uC-uR, 
   % tj. L*di/dt = U0-uC-uR => di/dt = (U0-uC-uR)/L
   % pro t = 0+ dostaneme: di/dt = (U0-uC(0))/L - uR/L = (U0-uC(0))/L - R*i(0)/L
   % po dosazeni pocatecnich podminek: di/dt = U0/L
   
   % t je sloupcovy vektor
   % proud je matice se 2 sloupci, 
   % v 1. sloupci je hledany proud, v 2.sloupci jeho derivace
   
subplot(2,2,1) % rozdeleni grafickeho okna na 4 casti (nasledujici graf v 1. casti)
plot(t,proud(:,1)) % vykresleni dvourozmerneho grafu
xlabel('t [s]') % popis osy x
ylabel('i [A]') % popis osy y
legend('i [A]') % popis krivky (legenda)

% odezva napeti na indukcnosti
uL=L*proud(:,2); % ze vztahu: uL = L*di/dt, di/dt ulozeno v  2. sloupci matice proud
% odezva napeti na odporu
uR=R*proud(:,1); % ze vztahu: uR = R*i
% odezva napeti na kondenzatoru
uC=U0-uR-uL; % z 2. Kirchhoffova zakona: uR + uL + uC =  U0

subplot(2,2,2) % rozdeleni grafickeho okna na 4 casti (nasledujici graf v 2. casti)
plot(t,uC,'r') % vykresleni dvourozmerneho grafu, r - zvolena cervena barva krivky
xlabel('t [s]') % popis osy x
ylabel('u_C [V]') % popis osy y, znak _ vytvori dolni index
legend('u_C [V]') % popis krivky (legenda), znak _ vytvori dolni index

subplot(2,2,3) % rozdeleni grafickeho okna na 4 casti (nasledujici graf v 3. casti)
plot(t,uR,'g') % vykresleni dvourozmerneho grafu, g - zvolena zelena barva krivky
xlabel('t [s]') % popis osy x
ylabel('u_R [V]') % popis osy y, znak _ vytvori dolni index
legend('u_R [V]') % popis krivky (legenda), znak _ vytvori dolni index

subplot(2,2,4) % rozdeleni grafickeho okna na 4 casti (nasledujici graf v 4. casti)
plot(t,uL,'m') % vykresleni dvourozmerneho grafu, m - zvolena fialova barva krivky
xlabel('t [s]') % popis osy x
ylabel('u_L [V]') % popis osy y, znak _ vytvori dolni index
legend('u_L [V]') % popis krivky (legenda), znak _ vytvori dolni index