<p> Numerical Inverse Laplace Transform ... [a,ns,nd] </p>

2b089dbf-0b0f-415e-963d-1f6a0d8c0ca9 3 3 12 decimal

&[DATE] &[TIME] - &[FILENAME]

in program ... user a 6 : ns 20 : nd 19 : 13 mat

s F t invLaplace t 1 el 105 -^: t.ini 0 ≡ t cases [a,ns,nd] user comman line n 1 ns nd + 1 + (range α n el a n 1 - (π * i * + : β n el e a^1 - (n^* - : 21 line for bdif 1 el nd 1 + Gamma nd 2 + 1 - Gamma / 1 Gamma / (: β2 β ns 2 + (ns 1 + nd + (range 11 range el : n 2 nd range bdif n el nd 1 + Gamma nd 2 + n - Gamma / n Gamma / (bdif n 1 - el + : for bdif bdif reverse 2 nd^/ : n 1 nd range β2 n el β2 n el bdif n el * : for β β 1 ns 1 + (range 11 range el β2 stack : β 1 el β 1 el 2 / : 41 line kt 1 pts range tt t kt el : s α tt / : βt β tt / : n 1 α length range βtF n el βt n el s n el F * : for RβtF βtF Re : ft kt el RβtF eval sum : 61 line for OUT t ft augment : 91 line :

Parametric algo style adapted from Oscar

s F 1 L / s s 2^s R L / * + 1 L C * / + / (* :

x c 10 x ≤ cases :

t XFR s F t 0 > assume s t invlaplace maple :

========================================

user range/mesh t.ini 0 : t.end 3 : pts 100 : 13 mat

in-situ incremental t t.ini t.end t.end t.ini - (pts 1 - / range :

R 1 : L 1 : C 0.01 : 13 mat

x c 10 x ≤ cases :

Export s F t invLaplace :

User XFR parameters [H, γ, T]

Export x XFR x c * 2 1 sys

Export 01 0.03 0.926 0.061 0.77 0.091 0.55 0.121 0.288 0.152 0.007 0.182 0.266 - 72 mat

R t * 2 L * / - exp 4 L * - R 2^C * + L 2^C * / - sqrt t * 2 / cos 4 L * - R 2^C * + (* 4 L * - R 2^C * + L 2^C * / - sqrt R * 4 L * - R 2^C * + L 2^C * / - sqrt t * 2 / sin * L * C * + (* 4 L * - R 2^C * + (L * /