パイソニスタの方へ
~~２点教えて下さい。~~ 1点教えて下さい。

~~(1) 本文「sympyで(T氏の数学日記様の方法を参考に)」内の、掛け算を実行しない。~~
(2023/08/29済) 私は、expand()の使い方を間違っていました。()の追加です。ページ中盤。

(2)　本文「sympyで(T氏の数学日記様の方法を参考に)」内の、
　　表現「0<2α<π/4」の表現方法について。
参考
＞不等式をやめて、以下の処理で区間に変換しています

オリジナル

上と同じです。大学入試数学問題集成様＞【3】テキスト

xxxx ChatGPTで(できませんでした。) xxx

2通りでやってみました。

正の数α,β,a,bが2*α+β=pi/4,tan(α)=1/a,tan(β)=1/bを満たすとき,aを用いてbを表しなさい.

？？？？＞これで、aを用いてbを表す式が得られました。

Solve[{α>0,β>0,a>0,b>0,2*α+β=pi/4, tan(α)=1/a,tan(β)=1/b}, {α,β,b}]

＞You can solve this system of equations using Mathematica or a similar computational software. Here's how you can do it in Mathematica:.......

＞（DepL翻訳)この連立方程式は，Mathematica または同様の計算ソフトを使って解くことができます．Mathematica での解き方は次のとおりです．

ChatGPT先生の言われるとおりでした。

xxxx WolframAlphaで(できませんでした。) xxx

Solve[{α>0,β>0,a>0,b>0,2*α+β==pi/4, tan(α)==1/a,tan(β)==1/b}, {α,β,b}]

結果
(整数解はありません)

xxx sympyで(T氏の数学日記様の方法を参考に) xxx　できませんでした。

＞αの範囲からaの範囲を厳密に調べ，最後の結論で分母が0にならないことを確かめる...

sympy.simplify.fu.TR12(rv, first=True)
Separate sums in tan.
https://docs.sympy.org/latest/modules/simplify/fu.html#sympy.simplify.fu.TR12

(2023/08/29済) 私は、expand()の使い方を間違っていました。
eqA.lhsmy1.expand()→(eqA.lhsmy1).expand()
掛け算にカッコが必要でした。eqB=を直しました。

from sympy import *
from sympy.simplify.fu import TR12
var('α,β,a,b',real=True,nonnegative=True)
eq1=Eq(2*α+β ,pi/4)                                                                 ;print("#eq1",eq1)
eq2=Eq(tan(α), 1/a)                                                                 ;print("#eq2",eq2)
eq3=Eq(tan(β), 1/b)                                                                 ;print("#eq3",eq3)
mySubs={eq2.lhs:eq2.rhs,eq3.lhs:eq3.rhs}
# 
eq4=Eq(  tan(eq1.lhs)                      ,tan(eq1.rhs)                         )  ;print("#eq4",eq4)
eq5=Eq( TR12(eq4.lhs)                      ,eq4.rhs                              )  ;print("#eq5",eq5)
eq6=Eq(numer(eq5.lhs)                      ,eq5.rhs*denom(eq5.lhs)               )  ;print("#eq6",eq6)
eq7=Eq(eq6.lhs.expand(trig=True).simplify(),eq6.rhs.expand(trig=True).simplify() )  ;print("#eq7",eq7)
eq8=Eq(eq7.lhs*denom(eq7.lhs)              ,eq7.rhs*denom(eq7.lhs)               )  ;print("#eq8",eq8)
eq9=Eq(eq8.lhs.subs(mySubs)                ,eq8.rhs.subs(mySubs)                 )  ;print("#eq9",eq9)
eqA=Eq(eq9.lhs.expand()                    ,eq9.rhs.expand()                     )  ;print("#eqA",eqA)
my1=a**2*b
eqB  =Eq((eqA.lhs*my1).expand()            ,(eqA.rhs*my1).expand()               )  ;print("#eqA",eqA)
eqC  =Eq( eqB.lhs*(-1)                     , eqB.rhs*(-1)                        )  ;print("#eqC",eqC)
ans  =solve(eqC,b)[0]
eqD  =Eq( numer(ans)                       ,denom(ans)*b                         )  ;print("#eqD",eqD)              
# 2023/08/28
# この間は、勉強中です。Interval(0,2*α) and Interval(2*α,pi/4)
# ここで、~
# したがって、...
eqE  =Eq( b                                ,solve(eqD,b)[0]                      )  ;print("#eqE",eqE)              
eqF  =Eq( b                                ,-numer(eqE.rhs)/-denom(eqE.rhs)      )  ;print("#eqF",eqF)              
#eq1 Eq(2*α + β, pi/4)
#eq2 Eq(tan(α), 1/a)
#eq3 Eq(tan(β), 1/b)
#eq4 Eq(tan(2*α + β), 1)
#eq5 Eq((tan(2*α) + tan(β))/(-tan(2*α)*tan(β) + 1), 1)
#eq6 Eq(tan(2*α) + tan(β), -tan(2*α)*tan(β) + 1)
#eq7 Eq(((tan(α)**2 - 1)*tan(β) - 2*tan(α))/(tan(α)**2 - 1), (tan(α)**2 + 2*tan(α)*tan(β) - 1)/(tan(α)**2 - 1))
#eq8 Eq((tan(α)**2 - 1)*tan(β) - 2*tan(α), tan(α)**2 + 2*tan(α)*tan(β) - 1)
#eq9 Eq((-1 + a**(-2))/b - 2/a, -1 + 2/(a*b) + a**(-2))
#eqA Eq(-1/b - 2/a + 1/(a**2*b), -1 + 2/(a*b) + a**(-2))
#eqB Eq(-a**2 - 2*a*b + 1, -a**2*b + 2*a + b)
#eqC Eq(a**2 + 2*a*b - 1, a**2*b - 2*a - b)
#eqD Eq(-a**2 - 2*a + 1, b*(-a**2 + 2*a + 1))
#eqE Eq(b, (-a**2 - 2*a + 1)/(-a**2 + 2*a + 1))
#eqF Eq(b, (a**2 + 2*a - 1)/(a**2 - 2*a - 1))

sympyで(α,β,bの3元連立方程式で)

・がんばれば、3行で？できました。
・π/8を使っていません。おかしいですね。

from sympy import *
var('α,β,a,b',real=True,nonnegative=True)
print("#",solve([Eq(2*α+β ,pi/4),
                 Eq(tan(α), 1/a),
                 Eq(tan(β), 1/b)],[α,β,b])[0][2].expand(trig=True).simplify())
# (a**2 + 2*a - 1)/(a**2 - 2*a - 1)

上と同じです。

from sympy import *
var('α,β,a,b',real=True,nonnegative=True)
eq1=Eq(2*α+β ,pi/4)
eq2=Eq(tan(α), 1/a)
eq3=Eq(tan(β), 1/b)
ans=solve([eq1,eq2,eq3],[α,β,b])
print("# α=",ans[0][0])
print("# β=",ans[0][1])
print("# b=",ans[0][2])
print()
print("#",ans[0][2])
print("#",ans[0][2].simplify())
print("#",ans[0][2].expand(trig=True))
print("#",ans[0][2].expand(trig=True).simplify())
# α= atan(1/a)
# β= -2*atan(1/a) + pi/4
# b= tan(2*atan(1/a) + pi/4)

# tan(2*atan(1/a) + pi/4)
# tan(2*atan(1/a) + pi/4)
# 2/(a - 2*a/(a - 1/a) + 2/(a**2 - 1) - 1/a) + 1/(1 - 2/(a - 1/a))
# (a**2 + 2*a - 1)/(a**2 - 2*a - 1)

sympyの実行環境

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

参考

＞expandは式の展開をしますが、そのときtrig=Trueと指定することにより、加算公式の計算をすることができます。

以下、いつもの?おすすめです。