function dI=RLCobvod(t,I)
% diferencialni rovnice 2. radu s konstantnimi koeficienty
% pro reseni serioveho obvodu R,L,C s konstantnim zdrojem napeti,
% obvod je resen volanim funkce RLC_reseni.m

% 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 ___     _L _  _   C||
%    ---|___|---/ \/ \/ \---||-         podle 2. Kirchhoffova zakona:
%   |                       || |         uR  +  uL  + uC =  Uo
%   |                          |        R*i  +  L*i'+ 1/C*intgral do 0 do t z i(t)dt + uC(0)  =  Uo
%   |                          |        derivovanim dostaneme:
%   |     _                    |        R*i' + L*i'' + 1/C*i = 0, tj. dif. rovnice 2. radu
%   |____/ \________/ _________|                 i'' = (- 1/C*i - R*i') / L
%        \_/        t=0         v MATLABu:  i  oznacime I(1)
%        Uo                                 i' oznacime dI(1)
%                               Zavedeme substituci:
%                                                 i' = I(2), tj. vlastne dif. rovnice 1. radu
%                               v MATLABu: dI(1) = I(2)
%                               Dale pak    i'' oznacime dI(2)
%                               a dosadime-li do rovnice pro i'' plati:
%                                   dI(2) = (-1/C*I(1)-R*I(2))/L, tj. take dif. rovnice 1. radu
%                               Resime tedy soustavu dvou dif. rovnic 1. radu:
%                               dI(1) = I(2)
%                               dI(2) = (-1/C*I(1)-R*I(2))/L
R=200; % R=200 Ohm
L=0.5; % L=0,5 H
C=5e-6; % C=5 mikroF

dI=[I(2);(-1/C*I(1)-R*I(2))/L]; % sloupcovy vektor se 2 prvky [dI(1);dI(2)]
                                % sloupce oddeleny strednikem

% lze psat tez:
% dI=zeros(2,1); % dvouprvkovy sloupcovy vektor s nulovymi prvky, 
%                    % výstupem tohoto m-file musí být sloupcový vektor!
% dI(1)=I(2); % prvni prvek vektoru: prvni dif. rovnice 1. radu
% dI(2)=(-I(1)/C-R*I(2))/L; % druhy prvek vektoru: druha dif. rovnice 1. radu
