#+PROPERTY: header-args :session laplace_sage :kernel sage :pandoc yes :async yes

#+begin_src jupyter-python
  %display latex
  var("G, phi, gamma, delta, t, a, b, c, d, Omega, omega, T", domain=RR)
  var("z", domain=CC)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}z\]
:END:

#+begin_src jupyter-python
  W = gamma + I*delta
  alpha(t) = G * exp(-W*t - I * phi)
  alpha
#+end_src

#+RESULTS:

#+begin_src jupyter-python
  im_alpha = (imag(alpha))
  im_alpha
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}t \ {\mapsto}\ G e^{\left(-\gamma
t\right)} \sin\left(-\delta t - \phi\right)\]
:END:

#+begin_src jupyter-python
  im_alpha.laplace(t, z)
#+end_src

#+RESULTS:
:     -G⋅(δ⋅cos(φ) + γ⋅sin(φ) + z⋅sin(φ))
: t ↦ ────────────────────────────────────
:              2    2            2
:             δ  + γ  + 2⋅γ⋅z + z


#+begin_src jupyter-python
  matrix([[z, -Omega], [Omega + a, z]]).inverse().simplify_full(
)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rr}
\frac{z}{\Omega^{2} + \Omega a + z^{2}} & \frac{\Omega}{\Omega^{2} +
\Omega a + z^{2}} \\ -\frac{\Omega + a}{\Omega^{2} + \Omega a + z^{2}} &
\frac{z}{\Omega^{2} + \Omega a + z^{2}} \end{array}\right)\]
:END:


#+begin_src jupyter-python
  matrix([[0, 0], [1, 0]]) * matrix([[a, b], [c, d]])
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rr} 0 & 0 \\ a &
b \end{array}\right)\]
:END:

#+begin_src jupyter-python
  matrix([[0, 1], [-1, 0]]) * matrix([[a, b], [c, d]])
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rr} c & d \\ -a
& -b \end{array}\right)\]
:END:

#+begin_src jupyter-python
  assume(T>0)
  integrate(exp(I* omega * t) * exp(-t * a) * sin(delta * t + phi), t, 0, T, algorithm='giac').simplify_full()
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}\frac{{\left(a^{2} \delta e^{\left(T
a\right)} + \delta^{3} e^{\left(T a\right)} + 2 i \, a \delta \omega
e^{\left(T a\right)} - \delta \omega^{2} e^{\left(T a\right)} +
{\left(-i \, \omega^{3} \sin\left(T \delta\right) + {\left(\delta
\cos\left(T \delta\right) - a \sin\left(T \delta\right)\right)}
\omega^{2} + {\left(-2 i \, a \delta \cos\left(T \delta\right) +
{\left(-i \, a^{2} + i \, \delta^{2}\right)} \sin\left(T
\delta\right)\right)} \omega - {\left(a^{2} \delta + \delta^{3}\right)}
\cos\left(T \delta\right) - {\left(a^{3} + a \delta^{2}\right)}
\sin\left(T \delta\right)\right)} \cos\left(T \omega\right) +
{\left(\omega^{3} \sin\left(T \delta\right) + {\left(i \, \delta
\cos\left(T \delta\right) - i \, a \sin\left(T \delta\right)\right)}
\omega^{2} + {\left(2 \, a \delta \cos\left(T \delta\right) +
{\left(a^{2} - \delta^{2}\right)} \sin\left(T \delta\right)\right)}
\omega + {\left(-i \, a^{2} \delta - i \, \delta^{3}\right)} \cos\left(T
\delta\right) + {\left(-i \, a^{3} - i \, a \delta^{2}\right)}
\sin\left(T \delta\right)\right)} \sin\left(T \omega\right)\right)}
\cos\left(\phi\right) + {\left(a^{3} e^{\left(T a\right)} + a \delta^{2}
e^{\left(T a\right)} + a \omega^{2} e^{\left(T a\right)} + i \,
\omega^{3} e^{\left(T a\right)} + {\left(i \, a^{2} e^{\left(T
a\right)} - i \, \delta^{2} e^{\left(T a\right)}\right)} \omega +
{\left(-i \, \omega^{3} \cos\left(T \delta\right) - {\left(a \cos\left(T
\delta\right) + \delta \sin\left(T \delta\right)\right)} \omega^{2} +
{\left(2 i \, a \delta \sin\left(T \delta\right) + {\left(-i \, a^{2} +
i \, \delta^{2}\right)} \cos\left(T \delta\right)\right)} \omega -
{\left(a^{3} + a \delta^{2}\right)} \cos\left(T \delta\right) +
{\left(a^{2} \delta + \delta^{3}\right)} \sin\left(T
\delta\right)\right)} \cos\left(T \omega\right) + {\left(\omega^{3}
\cos\left(T \delta\right) + {\left(-i \, a \cos\left(T \delta\right) - i
\, \delta \sin\left(T \delta\right)\right)} \omega^{2} - {\left(2 \, a
\delta \sin\left(T \delta\right) - {\left(a^{2} - \delta^{2}\right)}
\cos\left(T \delta\right)\right)} \omega + {\left(-i \, a^{3} - i \, a
\delta^{2}\right)} \cos\left(T \delta\right) + {\left(i \, a^{2} \delta
+ i \, \delta^{3}\right)} \sin\left(T \delta\right)\right)} \sin\left(T
\omega\right)\right)} \sin\left(\phi\right)}{a^{4} e^{\left(T a\right)}
+ 2 \, a^{2} \delta^{2} e^{\left(T a\right)} + \delta^{4} e^{\left(T
a\right)} + \omega^{4} e^{\left(T a\right)} + 2 \, {\left(a^{2}
e^{\left(T a\right)} - \delta^{2} e^{\left(T a\right)}\right)}
\omega^{2}}\]
:END:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-\frac{{\left(G c \delta \omega^{2}
e^{\left(T a + T \gamma\right)} - 2 \, {\left(-i \, G a c \delta
e^{\left(T a\right)} - i \, G c \delta \gamma e^{\left(T
a\right)}\right)} \omega e^{\left(T \gamma\right)} - {\left(G a^{2} c
\delta e^{\left(T a\right)} + G c \delta^{3} e^{\left(T a\right)} + 2 \,
G a c \delta \gamma e^{\left(T a\right)} + G c \delta \gamma^{2}
e^{\left(T a\right)}\right)} e^{\left(T \gamma\right)} + {\left(G c
\gamma^{3} \sin\left(T \delta\right) + i \, G c \omega^{3} \sin\left(T
\delta\right) + {\left(G c \delta \cos\left(T \delta\right) + 3 \, G a c
\sin\left(T \delta\right)\right)} \gamma^{2} - {\left(G c \delta
\cos\left(T \delta\right) + 3 \, G a c \sin\left(T \delta\right) + 3 \,
G c \gamma \sin\left(T \delta\right)\right)} \omega^{2} + {\left(2 \, G
a c \delta \cos\left(T \delta\right) + {\left(3 \, G a^{2} c + G c
\delta^{2}\right)} \sin\left(T \delta\right)\right)} \gamma + {\left(-2
i \, G a c \delta \cos\left(T \delta\right) - 3 i \, G c \gamma^{2}
\sin\left(T \delta\right) - 2 \, {\left(i \, G c \delta \cos\left(T
\delta\right) + 3 i \, G a c \sin\left(T \delta\right)\right)} \gamma +
{\left(-3 i \, G a^{2} c - i \, G c \delta^{2}\right)} \sin\left(T
\delta\right)\right)} \omega + {\left(G a^{2} c \delta + G c
\delta^{3}\right)} \cos\left(T \delta\right) + {\left(G a^{3} c + G a c
\delta^{2}\right)} \sin\left(T \delta\right)\right)} e^{\left(i \, T
\omega\right)}\right)} \cos\left(\phi\right) + {\left(-i \, G c
\omega^{3} e^{\left(T a + T \gamma\right)} + 3 \, {\left(G a c
e^{\left(T a\right)} + G c \gamma e^{\left(T a\right)}\right)}
\omega^{2} e^{\left(T \gamma\right)} + {\left(3 i \, G a^{2} c
e^{\left(T a\right)} + i \, G c \delta^{2} e^{\left(T a\right)} + 6 i \,
G a c \gamma e^{\left(T a\right)} + 3 i \, G c \gamma^{2} e^{\left(T
a\right)}\right)} \omega e^{\left(T \gamma\right)} - {\left(G a^{3} c
e^{\left(T a\right)} + G a c \delta^{2} e^{\left(T a\right)} + 3 \, G a
c \gamma^{2} e^{\left(T a\right)} + G c \gamma^{3} e^{\left(T a\right)}
+ {\left(3 \, G a^{2} c e^{\left(T a\right)} + G c \delta^{2} e^{\left(T
a\right)}\right)} \gamma\right)} e^{\left(T \gamma\right)} + {\left(G c
\gamma^{3} \cos\left(T \delta\right) + i \, G c \omega^{3} \cos\left(T
\delta\right) + {\left(3 \, G a c \cos\left(T \delta\right) - G c \delta
\sin\left(T \delta\right)\right)} \gamma^{2} - {\left(3 \, G a c
\cos\left(T \delta\right) + 3 \, G c \gamma \cos\left(T \delta\right) -
G c \delta \sin\left(T \delta\right)\right)} \omega^{2} - {\left(2 \, G
a c \delta \sin\left(T \delta\right) - {\left(3 \, G a^{2} c + G c
\delta^{2}\right)} \cos\left(T \delta\right)\right)} \gamma + {\left(-3
i \, G c \gamma^{2} \cos\left(T \delta\right) + 2 i \, G a c \delta
\sin\left(T \delta\right) - 2 \, {\left(3 i \, G a c \cos\left(T
\delta\right) - i \, G c \delta \sin\left(T \delta\right)\right)} \gamma
+ {\left(-3 i \, G a^{2} c - i \, G c \delta^{2}\right)} \cos\left(T
\delta\right)\right)} \omega + {\left(G a^{3} c + G a c
\delta^{2}\right)} \cos\left(T \delta\right) - {\left(G a^{2} c \delta +
G c \delta^{3}\right)} \sin\left(T \delta\right)\right)} e^{\left(i \, T
\omega\right)}\right)} \sin\left(\phi\right)}{\omega^{4} e^{\left(T a +
T \gamma\right)} - 4 \, {\left(-i \, a e^{\left(T a\right)} - i \,
\gamma e^{\left(T a\right)}\right)} \omega^{3} e^{\left(T
\gamma\right)} - 2 \, {\left(3 \, a^{2} e^{\left(T a\right)} +
\delta^{2} e^{\left(T a\right)} + 6 \, a \gamma e^{\left(T a\right)} + 3
\, \gamma^{2} e^{\left(T a\right)}\right)} \omega^{2} e^{\left(T
\gamma\right)} - 4 \, {\left(i \, a^{3} e^{\left(T a\right)} + i \, a
\delta^{2} e^{\left(T a\right)} + 3 i \, a \gamma^{2} e^{\left(T
a\right)} + i \, \gamma^{3} e^{\left(T a\right)} + {\left(3 i \, a^{2}
e^{\left(T a\right)} + i \, \delta^{2} e^{\left(T a\right)}\right)}
\gamma\right)} \omega e^{\left(T \gamma\right)} + {\left(a^{4}
e^{\left(T a\right)} + 2 \, a^{2} \delta^{2} e^{\left(T a\right)} +
\delta^{4} e^{\left(T a\right)} + 4 \, a \gamma^{3} e^{\left(T a\right)}
+ \gamma^{4} e^{\left(T a\right)} + 2 \, {\left(3 \, a^{2} e^{\left(T
a\right)} + \delta^{2} e^{\left(T a\right)}\right)} \gamma^{2} + 4 \,
{\left(a^{3} e^{\left(T a\right)} + a \delta^{2} e^{\left(T
a\right)}\right)} \gamma\right)} e^{\left(T \gamma\right)}}\]
:END:
#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}G c e^{\left(-{\left(a +
\gamma\right)} t\right)} \sin\left(\delta t + \phi\right)
\sin\left(\omega t + b\right)\]
:END:



#+begin_src jupyter-python
  var('t,s,r,l,u')
  var('P_k,L_k,B_n,C_n,B_m,C_m,G_l,W_l,Gc_l,Wc_l', domain=CC)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}\left(P_{k}, L_{k}, B_{n}, C_{n},
B_{m}, C_{m}, G_{l}, W_{l}, \mathit{Gc}_{l}, \mathit{Wc}_{l}\right)\]
:END:

#+begin_src jupyter-python
  α(t) = G_l * exp(-W_l * t)
  α_conj(t) = Gc_l * exp(-Wc_l * t)
  α_dot(t) = P_k * exp(-L_k * t)
  B_1(t) = B_n * exp(-C_n * t)
  B_2(t) = B_m * exp(-C_m * t)
  α_conj
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}t \ {\mapsto}\ \mathit{Gc}_{l}
e^{\left(-\mathit{Wc}_{l} t\right)}\]
:END:

#+begin_src jupyter-python
  inner = integrate(B_1(t-r-u) * α(u), u, 0, t-r) + integrate(B_1(t-r+u) * α_conj(u), u, 0, r)
  inner
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-B_{n} G_{l}
{\left(\frac{e^{\left(C_{n} r - C_{n} t\right)}}{C_{n} - W_{l}} -
\frac{e^{\left(W_{l} r - W_{l} t\right)}}{C_{n} - W_{l}}\right)} + B_{n}
\mathit{Gc}_{l} {\left(\frac{e^{\left(C_{n} r - C_{n} t\right)}}{C_{n} +
\mathit{Wc}_{l}} - \frac{e^{\left(-\mathit{Wc}_{l} r - C_{n}
t\right)}}{C_{n} + \mathit{Wc}_{l}}\right)}\]
:END:

#+begin_src jupyter-python
  assume(C_n/L_k != -1)
  assume(W_l/L_k != -1)
  assume(-Wc_l/L_k != -1)
  assume(-Wc_l/L_k != -1)
  assume(L_k/(L_k+C_m) != -1)
#+end_src

#+RESULTS:

#+begin_src jupyter-python
  whole = (B_2(s-r) * α_dot(t-s) * inner).integrate(r, 0, s).simplify_full()
#+end_src

#+RESULTS:

#+begin_src jupyter-python :results scalar
  %display plain
  integ = whole.integrate(s, 0, t, algorithm='giac')
#+end_src

#+RESULTS:

#+begin_src jupyter-python :results scalar
  import sympy
  from sympy.utilities.codegen import codegen
  integ_s = sympy.sympify(integ)
  result = codegen(("conv_part", integ_s), "F95")
  for name, contents in result:
      with open(name, 'w') as f:
          f.write(contents)
#+end_src

#+RESULTS:

#+begin_src jupyter-python
  %display latex
  (B_2(s-r) * α_dot(t-s)).integrate(s, r, t)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-B_{m} P_{k}
{\left(\frac{e^{\left(C_{m} r - C_{m} t\right)}}{C_{m} - L_{k}} -
\frac{e^{\left(L_{k} r - L_{k} t\right)}}{C_{m} - L_{k}}\right)}\]
:END:

#+begin_src jupyter-python
  %display latex
  (B_1(t-r-u) * α(u)).integrate(u, 0, t-r)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-B_{n} G_{l}
{\left(\frac{e^{\left(C_{n} r - C_{n} t\right)}}{C_{n} - W_{l}} -
\frac{e^{\left(W_{l} r - W_{l} t\right)}}{C_{n} - W_{l}}\right)}\]
:END:

#+begin_src jupyter-python
  %display latex
  (B_1(t-r+u) * α_conj(u)).integrate(u, 0, r)
#+end_src

#+RESULTS:
:RESULTS:

\[\newcommand{\Bold}[1]{\mathbf{#1}}B_{n} \mathit{Gc}_{l}
{\left(\frac{e^{\left(C_{n} r - C_{n} t\right)}}{C_{n} +
\mathit{Wc}_{l}} - \frac{e^{\left(-\mathit{Wc}_{l} r - C_{n}
t\right)}}{C_{n} + \mathit{Wc}_{l}}\right)}\]
:END:

#+begin_src jupyter-python
  assume(C_n/C_m != -1)

  assume(L_k/C_m+C_n/C_m-1 != -1)
#+end_src

#+RESULTS:

#+begin_src jupyter-python
  ((exp(-C_m * r)) * (exp(-W_l * r))).integrate(r, 0, t)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-\frac{e^{\left(-C_{m} t - W_{l}
t\right)}}{C_{m} + W_{l}} + \frac{1}{C_{m} + W_{l}}\]
:END:

#+begin_src jupyter-python
  ((exp(-L_k * (t-r))) * (exp(-Wc_l*r)*exp(-C_n*t))).integrate(r, 0, t)
#+end_src

#+RESULTS:

#+begin_src jupyter-python
  var('t, wc')
  %display latex
  diff(1/pi * (wc/(1+I*wc*t))^2, t)
#+end_src

#+RESULTS:
:RESULTS:
\[\newcommand{\Bold}[1]{\mathbf{#1}}-\frac{2 i \, \mathit{wc}^{3}}{\pi
{\left(i \, t \mathit{wc} + 1\right)}^{3}}\]
:END:


#+begin_src jupyter-python :results none

#+end_src