Re 4次方程式を解くプログラム

 投稿者:エス・テー  投稿日:2012年 3月26日(月)11時38分30秒
  山中先生、4次方程式を解くプログラム有難うございました。今回は複素数計算の威力を知りました。なお、計算結果が複素数の場合、例えば虚数部が10^-8以下なら、虚数部=0にして、実数部のみにすることはできないのでしょうか。ご教示よろしくお願いします。
 

Re: Re 4次方程式を解くプログラム

 投稿者:山中和義  投稿日:2012年 3月26日(月)12時45分10秒
  > No.1836[元記事へ]

エス・テーさんへのお返事です。

> 計算結果が複素数の場合、例えば虚数部が10^-8以下なら、虚数部=0にして、実数部のみにすることはできないのでしょうか。

複素数xの実部と虚部は、RE(x),IM(x)関数で取得できるので、次のようにすれば良いかと思います。

 :
 :

LET x1=y1-a3/(4*a4) !x=y-a3/(4*a4)
LET x2=y2-a3/(4*a4)
LET x3=y3-a3/(4*a4)
LET x4=y4-a3/(4*a4)


IF ABS(Im(x1))<1E-8 THEN LET x1=Re(x1)
IF ABS(Im(x2))<1E-8 THEN LET x2=Re(x2)
IF ABS(Im(x3))<1E-8 THEN LET x3=Re(x3)
IF ABS(Im(x4))<1E-8 THEN LET x4=Re(x4)


!解を表示する

PRINT x1
PRINT x2
PRINT x3
PRINT x4

END

 :
 :
 

戻る