四次方程式　quartic equation 要修正！    実行数: 26

  1. ホーム
  2.  / みんなの自作式（科学）
  3.  / 化学（一般）

四次方程式
解の公式（全文）式の一部を置き換えたことにより簡略化したもの
Wikipedia の公式を使ったが結果が微妙に公式の値と異なる。
タイプミスの可能性大！　いつか修正したい！
function cubicRoot(x1)を使うことで修正してみたが
Input Errorで走らなくなった！

function cubicRoot(x1)
{
if (x1 > 0.0) {
return x1^(1/3);
} else {
return -(-x1)^(1/3);
}
}

b3 = -a3^3/(a4^3) + 4*a2*a3/a4^2 - 8*a1/a4;
b5 = a2^2 - 3*a1*a3 + 12*a0*a4;
b6 = 2*a2^3 - 9*a1*a2*a3 + 27*a0*a3^2 +27*a1^2*a4 - 72*a0*a2*a4;
b7 = -4*b5^3 + b6^2;
b8 = b6 + sqrt(b7);

crb8 = cubicRoot(b8);

b4 = a3^2/(4*a4^2) + 2^(1/3)*b5/(3*crb8*a4) +crb8/(3*2^(1/3)*a4) - 2*a2/(3*a4);
b2 = -b3/(b4)^(1/4) + a3^2/(2*a4^2) - 2^(1/3)*b5/(3*crb8*a4) -crb8/(3*2^(1/3)*a4) - 4*a2/(3*a4);
b1 = b3/b4^(1/4) + a3^2/(2*a4^2) - 2^(1/3)*b5/(3*crb8*a4) -crb8/(3*2^(1/3)*a4) - 4*a2/(3*a4);

x1 = -a3/(4*a4) - sqrt(b2)/2 - sqrt(b4)/2;
x2 = -a3/(4*a4) + sqrt(b2)/2 - sqrt(b4)/2;
x3 = -a3/(4*a4) - sqrt(b1)/2 + sqrt(b4)/2;
x4 = -a3/(4*a4) + sqrt(b1)/2 + sqrt(b4)/2;

print(b1);
print(b2);
print(b3);
print(b4);
print(b5);
print(b6);
print(b7);
print(b8);
print(crb8);
print(x1);
print(x2);
print(x3);
print(x4);
a4
    3. 炭酸水素ナトリウム pH
a3
    3. 0.000001 mol/L
a2
    3. Ka1*Ka2
a1
    3. 0.000001 mol/L
a0
    3. - Ka1*Ka2*Kw
b1
    1.  
b2
    1.  
b3
    1.  
b4
    1.  
b5
    1.  
b6
    1.  
b7
    1.  
b8
    1.  
crb8
    1.  
x1
    1.  
x2
    1.  
x3
    1.  
x4
    1.  
本ライブラリは会員の方が作成した作品です。 内容について当サイトは一切関知しません。
    四次方程式　quartic equation 要修正！
    [0-0] / 0件表示件数
    BACK NEXT
    メッセージは１件も登録されていません。

    BACK NEXT

    アンケートにご協力頂き有り難うございました。

    送信を完了しました。



    Back

    【 四次方程式　quartic equation 要修正！ 】のアンケート記入欄
    年齢
    職業

    この計算式は

    使用目的
    ご意見・ご感想・ご要望（バグ報告はこちら) バグに関する報告 (ご意見・ご感想・ご要望はこちら）
    計算バグ(入力値と間違ってる結果、正しい結果、参考資料など)
    説明バグ(間違ってる説明文と正しい説明文など)
    アンケートは下記にお客様の声として掲載させていただくことがあります。
    送信

    ブックマーク

    実行履歴