パイソニスタの方へ
1/tanの2乗の合計の計算(ページ後半)について、
pythonでfloat以外の方法を探しています。sympy.fu?
よろしくお願いします。
オリジナル
上と同じです。大学入試数学問題集成＞【４】テキスト
(参考) 2023/06/19追加 前期理系【４】
Mathematicaで
？？？ x1 = Cot ？？？ いい方法がありますか。
In:
s = Sin[7 \[Theta]]
f = Simplify[TrigExpand[s]/Sin[\[Theta]]^7] /.
Cot[\[Theta]]^2 -> x /. Cot[\[Theta]]^4 -> x^2 /.
Cot[\[Theta]]^6 -> x^3
ans = Solve[ s == 0 && 0 < \[Theta] < \[Pi]/2, \[Theta], Reals]
x1 = Cot[\[Theta]]^2 /. \[Theta] -> 1 \[Pi]/7;
x2 = Cot[\[Theta]]^2 /. \[Theta] -> 2 \[Pi]/7;
x3 = Cot[\[Theta]]^2 /. \[Theta] -> 3 \[Pi]/7;
Simplify[x1 + x2 + x3]
-1*(-35)/7
Out:
Sin[7 \[Theta]]
-1 + 21 x - 35 x^2 + 7 x^3
{{\[Theta] -> \[Pi]/7}, {\[Theta] -> (2 \[Pi])/7}, {\[Theta] -> (
3 \[Pi])/7}}
Cot[\[Pi]/7]^2 + Tan[\[Pi]/14]^2 + Tan[(3 \[Pi])/14]^2
5
WolframAlphaで
勉強中
sympyで
パイソニスタの方へ 短くなりますか？
方法がわからなかったので、subsにしました。
https://docs.sympy.org/latest/modules/simplify/fu.html
from sympy import *
θ =symbols('θ' ,real =True)
k =symbols('k' ,integer=True)
sin7θ=symbols('sin7θ',real =True)
x =symbols('x' ,real =True)
sinθ =symbols('sinθ' ,real =True)
cosθ =symbols('cosθ' ,real =True)
tanθ =symbols('tanθ' ,real =True)
k=7
sin7θ=sin(k*θ)
g=im((cos(θ)+I*sin(θ))**k)
h=g/sin(θ)**k
h=h.subs({sin(θ):sinθ,cos(θ):cosθ}).expand()
h=h.subs({cosθ**6/sinθ**6:1/tanθ**6,cosθ**4/sinθ**4:1/tanθ**4,cosθ**2/sinθ**2:1/tanθ**2})
h=h.subs({1/tanθ**2:x})
ans_θ=solveset(sin7θ,θ,Interval.open(0,pi/2))
print("#(1)",g)
print("#(2)",h)
print("#(2)",ans_θ)
print("#(3)",-h.coeff(x,2)/h.coeff(x,3) )
print("#",'-'*111)
myPol=7*((x-1/tan(ans_θ.args[0])**2)*(x-1/tan(ans_θ.args[1])**2)*(x-1/tan(ans_θ.args[2])**2))
myP =Poly(myPol.expand(),x).args[0]
myPx3=myP.coeff(x,3)
myPx2=myP.coeff(x,2)
myPx1=myP.coeff(x,1)
myPx0=myP.coeff(x,0)
print("#myPol" ," "*4,myPol)
print("#myP3 ",float(myPx3),myPx3 )
print("#myP2" ,float(myPx2),myPx2 )
print("#myP1 " ,float(myPx1),myPx1 )
print("#myP0 " ,float(myPx0),myPx0 )
#(1) -sin(θ)**7 + 21*sin(θ)**5*cos(θ)**2 - 35*sin(θ)**3*cos(θ)**4 + 7*sin(θ)*cos(θ)**6
#(2) 7*x**3 - 35*x**2 + 21*x - 1
#(2) {pi/7, 2*pi/7, 3*pi/7}
#(3) 5
# ---------------------------------------------------------------------------------------------------------------
#myPol 7*(x - 1/tan(pi/7)**2)*(x - 1/tan(2*pi/7)**2)*(x - 1/tan(3*pi/7)**2)
#myP3 7.0 7
#myP2 -35.0 -7/tan(pi/7)**2 - 7/tan(2*pi/7)**2 - 7/tan(3*pi/7)**2
#myP1 21.0 7/(tan(2*pi/7)**2*tan(3*pi/7)**2) + 7/(tan(pi/7)**2*tan(3*pi/7)**2) + 7/(tan(pi/7)**2*tan(2*pi/7)**2)
#myP0 -1.0 -7/(tan(pi/7)**2*tan(2*pi/7)**2*tan(3*pi/7)**2)
①私の環境は,pycharmです。
②よく聞くのは、Jupyterです。
③web上で、上記のソースを「SymPy Live shell」に、コピー貼り付けでもできました。
黒背景の右上に、マウスを移動すると、コピーマークが発生します。
??? タブレット環境で、コピー貼り付けが実行できませんでした。???
計算時間は、1分以上かかると思います。
参考
おすすめです。
goole検索
Mathematica
未整理