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

62d5566b-6fa6-405a-849a-e6b51962da92 7 3 12 decimal

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

in program ... user a 6 : ns 10 : nd 9 : 13 mat

Algo style adapted from Oscar

s F t invLaplace t 1 el 106 -^: 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 :

s F 12 H * s * (1 γ s * + (* 1 T s * + (* / :

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

t XFR γ T 1 - t T / - exp + (* + γ t γ / - exp * - 2 γ T - (* H * /

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

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 :

H 5 :

γ 0.2 :

T 0.5 :

x c 10 x ≤ cases :

XFR s F t invLaplace :

User XFR parameters

XFR x XFR x c * 2 1 sys