静力学と材料力学 発展例題1.16 p11「材料力学演習(20221021)」をsympyでやってみた。

オリジナル

発展例題1.16 p11
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=16

「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。＞sympy
https://qiita.com/mrrclb48z/items/c3274daf5f7a99cfe027#%E4%BD%9C%E6%A5%AD%E4%B8%AD%E7%A7%81%E3%81%AE%E8%A7%A3%E7%AD%94sympy

sympyで

「def myBibun」が怪しいです。アドバイスいただけると助かります。

from sympy import *
N,dN,dx,q,g,A,E,N0,x,L,P,e=symbols('N dN dx q g A E N0 x L P e',real=True)
a,b                       =symbols('a b'                       ,real=True)
def myBibun(f,dN,dx):
   return solve(f,dN)[0]/dx
myF=Eq(N+dN/dx*dx-N+q*g*A*dx,0)
Nx =integrate(myBibun(myF.lhs,dN,dx),x)+N0
Nx =Nx.subs({N0:solve(Eq(Nx.subs({x:L}),P),N0)[0]})
ex =Nx/(A*E).subs({x:e}).expand()
SGx=integrate(ex,x).expand()
print("#",SGx.subs({x:L}))
# L**2*g*q/(2*E) + L*P/(A*E)

①私の環境は,pycharmです。
②よく聞くのは、Jupyterです。
③web上で、上記のソースを「SymPy Live shell」に、コピー貼り付けでもできました。
黒背景の右上に、マウスを移動すると、コピーマークが発生します。
???　タブレット環境で、コピー貼り付けが実行できませんでした。???

Mathematicaで

勉強中

作図

勉強中