Mathematica vs. C++ 计算结果

Question

我有两个应该相同但给出不同结果的程序。 Mathematica 中的一个（给出正确的结果）和 C++ 中的一个（不正确）。

首先是数学：

q = 0.002344;
s = 0.0266;
v = 0.0744;
a = -q*PDCx^2;
b = s*PDCx - 2*q*PCLx*PDCx - PDCz;
c = -1*(PCLz + q*PCLx^2 - s*PCLx + v);
d = b*b - (4*a*c);

t = (-b + Sqrt[d])/(2*a)

现在是 C++：

long double q = 0.002344;
long double s = 0.0266;
long double v = 0.0744;
long double a = -q * pow(PDCx, 2);
long double b = s * PDCx - 2 * q*PCLx*PDCx - PDCz;
long double c = (-1.0)*(PCLz + q * pow(PCLx, 2) - s * PCLx + v);
long double d = b * b - 4.0 * a*c;

t = (-b + sqrtf(d))/(2.0*a);

和

long double PCLx = -1.816017;
long double PCLz = 0.056013;
long double PDCx = 0.005073;
long double PDCz = -0.998134;

对于每个案例。 Mathematica 结果为 t = 0.1867646081，C++ 结果为 t = 0.124776。 这是二次方程的“加”解。 负解的不同分别是 16549276.47723365 和 16549276.539223。 我怀疑我允许错误地舍入 C++ 结果。