Vstupni hodnoty coil, charge, field concentrator (fc):

M = Length[Lchargeji]
Rcoil = 0.003; R = 0.002; Rfc = 0.004; 
    Rch = 0.001;
Rchargej = Table[Rch*M, {r, M}];
Cap = 0.001;
i0 = 0; ij0 = 0; ucb0 = 100;
obort = {t, 0, 10}; obortau = {tau, 0, 10};
ij[t] = Table[(Subscript[ij, r])[t], {r, M}]
podmi = i[0] == i0;
podmijT = 
        Table[(Subscript[ij, r])[0] == ij0, {r, 1, M}];
podmik = ik[0] == 0;
v[t] = 0;

Rovnice pro napajeci obvod:

ucb[t] = 0.1*Exp[-t];
uc = R*i[t] - ucb[t];
Plot[ucb[t], Evaluate[obort], PlotRange -> All]

Rovnice pro napeti na civce, se zjednodusenym zadanim indukcnosti. Ceka na vstup z q,Q:

sumLchargeji[n_] := 
        Plus @@ (Lchargeji[[n]]*
              D[ij[t], t]);
sumLchargejcoilD := 
        Plus @@ (Lchargejcoil*D[ij[t], t]);
sumDLchargejcoil := 
        Plus @@ (D[Lchargejcoil, z]*ij[t]);
rovij[s_] := 
        0 == Rchargej[[s]]*ij[t][[s]] + 
            Lchargejcoil[[s]]*D[i[t], t] + 
            D[Lchargejcoil[[s]], z]*i[t]*v[t] + sumLchargeji[s] + 
            Lchargejfc[[s]]*D[ik[t], t] + 
            D[Lchargejfc[[s]], z]*ifc[t]*
              v[t];
rovijT = Table[rovij[r], {r, 1, M}];
oij = Union[rovijT, podmijT];
rovuc = 
        uc == Lcoil*D[i[t], t] + sumLchargejcoilD + v[t]*sumDLchargejcoil + 
            Lcoilfc*D[ik[t], t] + Rcoil*i[t];
ouc = {rovuc, podmi};
sumLchargejfcD := 
        Plus @@ (Lchargejfc*D[ij[t], t]);
sumDLchargejfc := 
        Plus @@ (D[Lchargejfc, z]*ij[t]);
rovik = 
        0 == 
          Rfc*ik[t] + Lfc*D[ik[t], t] + sumLchargejfcD + v[t]*sumDLchargejk + 
            Lcoilfc*D[i[t], t];
oik = {rovik, podmik};
oucii = Union[oij, ouc, oik]

reseneveliciny = Union[{i[t]}, {ik[t]}, ij[t]]

NDSolve[oucii, reseneveliciny, obort] // 
      Timing
res2 = %[[2]];

Plot[Evaluate[Evaluate[i[t] /. res2]], Evaluate[obort], 
        PlotRange -> All, 
        PlotLabel -> "i[t]"];
Plot[Evaluate[Evaluate[ik[t] /. res2]], Evaluate[obort], 
        PlotLabel -> "ik[t]"];
Plot[Evaluate[Evaluate[(Subscript[ij, 1])[t] /. res2]], Evaluate[obort], 
        PlotLabel -> "i1[t]"];
Plot[Evaluate[Evaluate[(Subscript[ij, 2])[t] /. res2]], Evaluate[obort], 
        PlotLabel -> "i2[t]"];
Plot[Evaluate[Evaluate[reseneveliciny /. res2]], Evaluate[obort], 
        PlotLabel -> "Vsechny i", 
        PlotRange -> All];