20/12/12 08:27:20.40 m4D75whP.net
>>940
最小二乗法じゃなくて、二次方程式の解の公式から
s=1/sqrt(ω-cos(ω)*sin(ω)+(sin(ω)/sin(ω+pi/3))^2*(2/3*pi-ω)+sin(ω)/sin(ω+pi/3)*cos(ω+pi/3)*sin(ω))
として
これを赤の面積の計算式に代入
(s*sin(ω)/sin(ω+2/3*pi))^2*(ω+2*pi/3) + s*sin(ω)/sin(ω+2/3*pi)*cos(pi-(ω+2*pi/3))*s*sin(ω) - (s^2*ω - s*cos(ω)*s*sin(ω))=2
一変数になったのでNewton-Raphson法で計算すると
> ω
[1] 0.21965843313038105
このとき
> s
[1] 2.8908840553803468
> BD=s*sin(ω)
> r=s*sin(ω)/sin(ω+pi/3)
> R=s*sin(ω)/sin(ω+2/3*pi)
> θ=ω+pi/3
> α=pi-θ
> β=θ+pi/3
> 2*s*ω + 2*r*(pi-θ) + 2*R*β
[1] 7.7048793147709809
やっと>801に辿りつけました。