y=5.1*10**16*(1-x**2+x**4)的なグラフをgnuplotで描こうとしているわけですね。

で、一応、グラフは描けるわけですよ。

でもね、x=0で10**10のオーダーになっているわけですよ。

10**16のオーダーじゃないとおかしいよね?って話なわけですよ。

ちなみに中性子星の密度について考えています。