Les équations du système sont:
> eq := {i1 * R*dx + U2 = U1, i1 = i2 + U2*I*C/dx*omega};
Pour obtenir 
en fonction de 
, il suffit encore
une fois de résoudre le système par rapport aux bonnes variables:
> solve(eq, {i2,U2});
Ce qui donne matriciellement
Autrement dit, 
avec:
Pour obtenir 
, il suffit d'écrire:
La relation générale est donc 
.
L'application numérique se fait donc immédiatement:
> AN := {R=10, C=10^(-7), omega=100*Pi, dx=10^(-2)};
> Man := evalm(subs(AN,M));
> evalm(Man^100);
