Question

我试图创建一个二次方程求解器但是当我输入一个大于1的系数时它似乎没有工作？代码和错误消息如下。非常感谢任何帮助。

print "Welcome to the quadratic equation solver."
print "The general quadratic equation = ax^2 + bx + c.\n"

def POS(a,b):
#This function gives the point at which the quadratic turns
    tp = float((-b)/(a*2))
    return (tp)

#This allows for the user to input the values of the variables
while True:
    s = raw_input ("Please insert a numerical value for a: ")
    try:
        a = float(s)
        break
    except ValueError:
        print "Please enter a numeric value."
print ("Well done.\n")

while True:
    s = raw_input ("Please insert a numerical value for b: ")
    try:
        b = float(s)
        break
    except ValueError:
        print "Please enter a numeric value."
print ("Well done.\n")

while True:
    s = raw_input ("Please insert a numerical value for c: ")
    try:
        c = float(s)
        break
    except ValueError:
        print "Please enter a numeric value."
print ("Well done.\n")

#This uses the function to give the co-ordinate of the turning point
print POS(a,b), "is the x-value of the turning point"
print ((a)*(POS(a,b)**2))+((b)*POS(a,b))+c, "is they y-value of the turning point. \n"

#This tells whether the quadratic is positive or negative
if a >0:
    print "The quadratic is positive.\n"

if a<0:
    print "The quadratic is negative.\n"

#This determines the root of the quadratic
root1 = (-b +((b**2) - (4*a*c))**0.5) / (2 * a)
root2 = (-b -((b**2) - (4*a*c))**0.5) / (2 * a)
print "The quadratic has a root at x =",root1
if root1 != root2:
    print "The quadratic has another root at x =",

#This uses the discriminant to determine the nature of the quadratic
if (b**2)-(4*a*c) == 0:
    print root1

elif (b**2)-(4*a*c) > 0:
    print root1 and root2

elif (b**2)-(4*a*c) < 0:
    print "The quadratic contains no roots"

for x in range(-1000,1000):
    quadratic = (a*x**2)+(b*x)+c

#Determines the derivitive and second derivitive of the quadratic
print "The derivitive of the quadratic is: ",2*a,"x", "+",b
print "The second derivitive of the quadratic is"

#Create plot
X = arange (-1000,1000,1)
plot(X, quadratic, linewidth=3, label="quadratic")

#Create title and axes label
title ("Graph of the quadratic")
xlabel ("x")
ylabel ("y")

#grid
grid(True)
legend()
show()

错误讯息：

Traceback (most recent call last):
  File "C:\Users\Peter\Downloads\test.py", line 49, in <module>
    root1 = (-b +((b**2) - (4*a*c))**0.5) / (2 * a)
ValueError: negative number cannot be raised to a fractional power

Answer 1

部分代码：

#This determines the root of the quadratic
root1 = (-b +((b**2) - (4*a*c))**0.5) / (2 * a)
root2 = (-b -((b**2) - (4*a*c))**0.5) / (2 * a)
print "The quadratic has a root at x =",root1
if root1 != root2:
    print "The quadratic has another root at x =",

#This uses the discriminant to determine the nature of the quadratic
if (b**2)-(4*a*c) == 0:
    print root1

elif (b**2)-(4*a*c) > 0:
    print root1 and root2

elif (b**2)-(4*a*c) < 0:
    print "The quadratic contains no roots"

for x in range(-1000,1000):
    quadratic = (a*x**2)+(b*x)+c

让我们用(b**2)-(4*a*c)替换D：

#This determines the root of the quadratic
D = (b**2)-(4*a*c)

root1 = (-b +(D)**0.5) / (2 * a)
root2 = (-b -(D)**0.5) / (2 * a)

print "The quadratic has a root at x =", root1
if root1 != root2:
    print "The quadratic has another root at x =",

#This uses the discriminant to determine the nature of the quadratic
if D == 0:
    print root1

elif D > 0:
    # old code:
    # print root1 and root2
    # new code:
    print root1, root2

elif D < 0:
    print "The quadratic contains no roots"

for x in range(-1000,1000):
    quadratic = (a*x**2)+(b*x)+c

这两行存在问题：

root1 = (-b +(D)**0.5) / (2 * a)
root2 = (-b -(D)**0.5) / (2 * a)

如果D小于0，则此行会引发ValueError。由于错误消息表示数字不能提升到分数幂。所以我们必须检查D是否小于0。

#This determines the root of the quadratic
D = (b**2)-(4*a*c)

# new code:
if D >= 0:
    root1 = (-b +(D)**0.5) / (2 * a)
    root2 = (-b -(D)**0.5) / (2 * a)

    print "The quadratic has a root at x =", root1
    if root1 != root2:
        print "The quadratic has another root at x =", root2


#This uses the discriminant to determine the nature of the quadratic
# We've already printed root1 and root2
# if D == 0:
#     print root1

# elif D > 0:
#     print root1, root2

# D < 0
else:
    print "The quadratic contains no roots"

for x in range(-1000,1000):
    quadratic = (a*x**2)+(b*x)+c

Answer 2

您也可以切换到Python3，以便代替ValueError您的号码implicitly cast to a complex number.

(-2)**.2
(0.9293164906031477+0.6751879523998812j)

Answer 3

#This uses the discriminant to determine the nature of the quadratic
if (b**2)-(4*a*c) == 0:
    print root1

elif (b**2)-(4*a*c) > 0:
    print root1 and root2

elif (b**2)-(4*a*c) < 0:
    print "The quadratic contains no roots"

在处理根之前执行此检查，如果它是复杂的根，则执行此检查。

>>> a=2
>>> b=4
>>> c=1
>>> delta=math.pow(b,2)-4*a*c
>>> math.sqrt(delta) #raises error for -ve integers
3.4641016151377544 
gamma=math.sqrt(delta)
>>>root1=(-b+gamma)/2/a
>>>root2=(-b-gamma)/2/a

Answer 4

我认为代码的问题是由于以下部分引起的：

#This determines the root of the quadratic
root1 = (-b +((b**2) - (4*a*c))**0.5) / (2 * a)
root2 = (-b -((b**2) - (4*a*c))**0.5) / (2 * a)
print "The quadratic has a root at x =",root1
if root1 != root2:
    print "The quadratic has another root at x =",

#This uses the discriminant to determine the nature of the quadratic
if (b**2)-(4*a*c) == 0:
    print root1

elif (b**2)-(4*a*c) > 0:

    print root1 and root2

elif (b**2)-(4*a*c) < 0:
    print "The quadratic contains no roots"

当提升到分数幂时，任何负数都将返回计算机无法计算的复数。

在(b**2) - (4*a*c))**0.5中，b**2可能小于4*a*c。在这种情况下会出现ValueError。

要防止这种情况：您应该按以下方式构建代码 -

if (b**2)-(4*a*c) == 0:
    print root1

elif (b**2)-(4*a*c) > 0:
    root1 = (-b +((b**2) - (4*a*c))**0.5) / (2 * a)
    root2 = (-b -((b**2) - (4*a*c))**0.5) / (2 * a)
    print "The quadratic has a root at x =",root1
    if root1 != root2:
       print "The quadratic has another root at x =",
    print root1 and root2

elif (b**2)-(4*a*c) < 0:
    print "The quadratic contains no roots"

二次方程求解器不起作用

4 个答案: