class d D2 { privatestaticreadonlydoubledsr3=math.sqrt (3； //summary ///盛金公式求解一元三次方程，设ax^3 bx^2 cx d=0。 //Summary//paramname='a '三阶系数/param ///param name='b '二阶系数/param ///param name='c '一阶系数/param//param 如果为null，则返回实际根/param ///param name='x2 '结果2/param ///param name='x3 '结果3/parampublicstaticVR x1，out double？ x2，out double？ x3，double dTol=1e-9 ) { x1=null； x2=空； x3=空； double A=b * b - 3 * a * c； double B=b * c - 9 * a * d； double C=c * c - 3 * b * d； double delta=B * B - 4 * A * C； if(IsZero(B(a，dTol ) iszero (b (dTol ) is zero ) b，dtol ) ) /盛金公式1 //方程式中有三重实根if (！ iszero(a，dTol ) ) doubleR1=-b/)3*a )； x1=x2=x3=r1； } else if (！ iszero(b，dTol ) ) { double r2=-c/b； x1=x2=x3=r2； } else if (！ iszero(c，dTol ) ) { double r3=-3 * d/c； x1=x2=x3=r3； (else(/无解) elseif(iszero(Delta，dTol ) ) /盛金公式)//方程有三个实根，其中有双根。 //A！=0 double K=B/A； x1=-b/a K； x2=x3=-K/2； (Elseif(Delta0) ) /盛金公式2 ) /方程有实根和共轭复根对(DoubleDSRDelta=math.sqrt ) Delta )； double yp=-B dsrdelta； double yn=-B - dsrdelta； double y1=A * b 3 * a/2 * yp； double y2=A * b 3 * a/2 * yn； const double dZhiShu=1d/3； doubley1Cr=math.pow(y1，dZhiShu )； doubley2Cr=math.pow(y2，dZhiShu )； double dAm3=3 * a； x1=-b-(y1Cry2Cr )； x1 /=dAm3； (else(/盛金公式4 )//方程有三个不等的实根。 //(A0，-1T1 ) doubleT=(2*A*B-3*A*B ) )2*math.sqrt ) math.pow(a，3 ) ) ) ) )； doubletheta=math.ACOS(t； double thetaOf1d3=theta/3； doublecosthetaof1D3=math.cos (theTaoF1D3 )，sintheta_1D3=math.sin ) the Tao f1 D3 )； double dx2 part=costhetaof1d3d sr3 * sin theta _1D3； double dx3 part=costhetaof1D3-DSR3* sin theta _1D3； doubledsra=math.sqrt(a )； double dAm3=3 * a； x1=-b - 2 * dsrA * cosThetaOf1d3； x2=-b dsrA * dx2part； x3=-b dsra * dx3部件； x1 /=dAm3； x2 /=dAm3； x3 /=dAm3； }//summary//num是否为0 . (由于双精度问题)//summary//param name=' num '/param//param name=' dtol ' returnspublicstaticbooliszero (双精度) }参考资料：
盛金公式(盛金公式(百度百科) ) ) ) ) ) )。