1

このコードの結果を、FindRoot ルーチンの解の位置にそれぞれ垂直線としてプロットしようとしました (そのようなもの: http://www.astrobio.net/albums/xsolar/ach.sized.jpg )。

    omega[a_] := 2 \[Pi] /a^(3/2);
Do[
Print[FindRoot[omega[a]/omega[5.2]==j/i+1, {a, 1}]],
{j,1,7,1},{i,1,7,1}]

誰か助けてくれませんか?

4

1 に答える 1

0

あなたの質問は十分に明確ではないので、私は推測しています:

omega[a_] := 2 Pi/a^(3/2);
ListPlot[Flatten[
  Table[{a /. FindRoot[omega[a]/omega[5.2] == j/i + 1, {a, 1}], i + j}, 
              {j, 1, 7, 1}, {i, 1, 7, 1}], 1], 
 Filling -> Axis,  PlotRange -> {{0, 6}, {0, 15}}]

Mathematica グラフィックス

于 2013-03-14T00:04:58.867 に答える