Skip to content

Commit

Permalink
draft
Browse files Browse the repository at this point in the history
  • Loading branch information
shenvitor committed Jun 4, 2024
1 parent 3a620ba commit e44262f
Showing 1 changed file with 43 additions and 1 deletion.
44 changes: 43 additions & 1 deletion docs/symbolics.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
"source": [
"import graphviz\n",
"import sympy as sp\n",
"from sympy import re, im, I, E, symbols, Sum"
"from sympy import re, im, I, E, symbols, Sum, IndexedBase, Function"
]
},
{
Expand Down Expand Up @@ -127,6 +127,48 @@
"tags": []
},
"outputs": [],
"source": [
"s, m_a2, Gamma_a2 = symbols(\"s m_a2 Gamma_a2\")\n",
"a_m = IndexedBase(\"a_m\")\n",
"Y2_m = Function(\"Y_2^m\")(\"Omega_1\")\n",
"m = symbols(\"m\", cls=sp.Idx)\n",
"\n",
"A12 = sp.Sum(a_m * Y2_m, (m, -2, 2)) / (s - m_a2**2 + I * m_a2 * Gamma_a2)\n",
"A12"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "8",
"metadata": {},
"outputs": [],
"source": [
"m_Delta, Gamma_Delta = symbols(\"m_Delta Gamma_Delta\")\n",
"b_m = IndexedBase(\"b_m\")\n",
"Y1_m = Function(\"Y_1^m\")(\"Omega_2\")\n",
"A23 = Sum(b_m * Y1_m, (m, -1, 1)) / (s - m_Delta**2 + I * m_Delta * Gamma_Delta)\n",
"A23"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "9",
"metadata": {},
"outputs": [],
"source": [
"c0, m_Nstar, Gamma_Nstar = symbols(\"c_0 m_N* Gamma_N*\")\n",
"A31 = c0 / (s - m_Nstar**2 + I * m_Nstar * Gamma_Nstar)\n",
"A31"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "10",
"metadata": {},
"outputs": [],
"source": []
}
],
Expand Down

0 comments on commit e44262f

Please sign in to comment.