1

コードの何が問題なのか知りたいのですが。私はMathematicaで非線形方程式のシステムを解こうとしています（最初はwolframでしたが、コマンドが長すぎました）：

Reduce[Pi*(h^2 + 2*R*(R - r))/sqrt (h^2 + (R - r)^2) - 2*x*Pi/3*h*R - 
     x*Pi/3*h*r == 0 && 
   Pi*(h^2 + 2*r*(r - R))/sqrt (h^2 + (R - r)^2) + 2*Pi*r - 
     x*Pi/3*h*R - 2*x*Pi/3*h*r == 0 && 
   Pi*h*(r + R)/sqrt (h^2 + (R - r)^2) - x*Pi/3*R^2 - x*Pi/3*R*r - 
     x*Pi/3*r^2 == 0 && -Pi/3*h*(R^2 + R*r + r^2) + 1 == 0, {R, r, h, 
   x}];

それを再入力してこれらの方程式を解く方法を知っていますか？ドキュメントに従って入力しようとしましたが、明らかに間違いを犯しました...

これらは元の方程式です（LaTeXでは、ここで正しく表示されるかどうかはわかりません：

\begin{equation*}
     \frac{\partial}{\partial R} L(R, r, h, \lambda) = \frac{\pi(h^2 + 2R(R-r))}{\sqrt{h^2 + (R - r)^2}} - 2\lambda \frac{\pi}{3}hR - \lambda \frac{\pi}{3}hr= 0
\end{equation*}

\begin{equation*}
     \frac{\partial}{\partial r} L(R, r, h, \lambda) = \frac{\pi(h^2 + 2r(r-R))}{\sqrt{h^2 + (R - r)^2}} + 2\pi r - \lambda \frac{\pi}{3}hR - 2\lambda \frac{\pi}{3}hr= 0
\end{equation*}

\begin{equation*}
     \frac{\partial}{\partial h} L(R, r, h, \lambda) = \frac{\pi h(r + R)}{\sqrt{h^2 + (R - r)^2}} - \lambda \frac{\pi}{3}R^2 - \lambda \frac{\pi}{3}Rr - \lambda \frac{\pi}{3}r^2= 0
\end{equation*}

\begin{equation*}
     \frac{\partial}{\partial \lambda} L(R, r, h, \lambda) = - \frac{\pi}{3} h (R^2 + Rr + r^2) + 1 = 0
\end{equation*}

編集：

円周率を円周率に修正しましたが、評価が開始されたので、間違いだったのかもしれません...非常に長い時間がかかります...

4

1 に答える 1

0

少なくとも基本を学ぶ必要があります。[ヘルプ]->[ドキュメントセンター]に移動し、検索バーで本をクリックします。最初からすべてが説明されています。

コメントですでに指摘されているように、すべての関数と組み込み記号は大文字で始まります。したがって、あなたの呼び出しは

Reduce[Pi*(h^2 + 2*R*(R - r))/Sqrt[h^2 + (R - r)^2] - 2*x*Pi/3*h*R - 
    x*Pi/3*h*r == 0 && 
  Pi*(h^2 + 2*r*(r - R))/Sqrt[h^2 + (R - r)^2] + 2*Pi*r - 
    x*Pi/3*h*R - 2*x*Pi/3*h*r == 0 && 
  Pi*h*(r + R)/Sqrt[h^2 + (R - r)^2] - x*Pi/3*R^2 - x*Pi/3*R*r - 
    x*Pi/3*r^2 == 0 && -Pi/3*h*(R^2 + R*r + r^2) + 1 == 0, {R, r, h, 
  x}]
于 2012-11-14T08:21:23.740 に答える