我尝试将此代码的结果绘制为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}]
有人能帮助我吗?
答案 0 :(得分: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}}]