function RC_reseni
% Reseni serioveho obvodu R,C s konstantnim zdrojem napeti
% Na vetev tvorenou seriove spojenym odporem 200 Ohm 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
C=5e-6; % C=5 mikroF
U0=100; % U0=100V
[t,uC]=ode45(@RC_dif_rov,[0,0.01],0);  
   % resi diferencialni rovnici danou v m-file RC_dif_rov
   % v casovem intervalu od 0 do 0.01 s pocatecni podminkou uC(0)=0

% proudova odezva a odezva napeti na odporu:
proud=(U0-uC)/R; % z 2. Kirchhoffova zakona: R*i + uC =  U0,
             % lze tez vypocitat ze vztahu:i=C*duC/dt
uR=U0-uC; % z 2. Kirchhoffova zakona: uR + uC =  U0,
          % lze tez vypocitat ze vztahu: uR = R*i

subplot(3,1,1)  % rozdeleni grafickeho okna na 3 casti (nasledujici graf v 1. 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(3,1,2)  % rozdeleni grafickeho okna na 3 casti (nasledujici graf v 2. casti)
plot(t,proud) % vykresleni dvourozmerneho grafu
xlabel('t [s]') % popis osy x
ylabel('i [A]') % popis osy y
legend('i [A]') % popis krivky (legenda)
subplot(3,1,3)  % rozdeleni grafickeho okna na 3 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
% napeti uC stoupa z 0 k hodnote U0=100V
% proud i klesa z hodnoty U0/R=100/200=0.5A na hodnotu 0
% napeti uR klesa z hodnoty U0=100V k 0
