写一元二次方程组的Python代码

这里有两种方法可以写一元二次方程组的Python代码。

方法一：使用numpy库

import numpy as np

# 例子：解方程组 x + y = 3, x - y = 1
# 将方程组转化为矩阵形式 Ax = b
A = np.array([[1, 1], [1, -1]])
b = np.array([3, 1])

# 使用numpy库中的solve函数求解方程组
x = np.linalg.solve(A, b)

print(x)  # 输出解向量 [2. 1.]

方法二：使用sympy库

import sympy

# 定义未知数
x, y = sympy.symbols('x y')

# 例子：解方程组 x + y = 3, x - y = 1
# 将方程组转化为sympy库中的方程形式
eq1 = sympy.Eq(x + y, 3)
eq2 = sympy.Eq(x - y, 1)

# 使用sympy库中的solve函数求解方程组
sol = sympy.solve((eq1, eq2), (x, y))

print(sol)  # 输出解字典 {x: 2, y: 1}

相关问题

写一个求解一元二次方程组的Python代码吧

好的，以下是求解一元二次方程组的Python代码：

import math

def solve_quadratic_equation(a, b, c):
    delta = b ** 2 - 4 * a * c
    if delta < 0:
        return None
    elif delta == 0:
        x = -b / (2 * a)
        return x
    else:
        x1 = (-b + math.sqrt(delta)) / (2 * a)
        x2 = (-b - math.sqrt(delta)) / (2 * a)
        return x1, x2

# 例子
a1, b1, c1 = 1, 2, -3
a2, b2, c2 = 2, -3, 1

x1, x2 = solve_quadratic_equation(a1, b1, c1)
print(f"x1={x1}, x2={x2}")

x1, x2 = solve_quadratic_equation(a2, b2, c2)
print(f"x1={x1}, x2={x2}")

希望能帮到你！

1. 一元二次方程求解： 1） 通过键盘输入一元二次方程的系数 2） 调用numpy库的roots函数实施求解 3） 绘制出对应的一元二次函数曲线 python

在Python中，解决一元二次方程通常涉及用户输入方程的系数（如a、b和c），然后利用科学计算库NumPy中的`linalg.solve()`函数来找到根，虽然直接提供`numpy.roots()`也是可行的。以下是基本步骤：

1. **输入系数**:

import numpy as np

# 获取用户输入的系数
a = float(input("请输入二次项系数(a): "))
b = float(input("请输入一次项系数(b): "))
c = float(input("请输入常数项系数(c): "))

2. **求解方程**:

# 计算判别式（为了判断有几个实根）
discriminant = b**2 - 4*a*c

# 使用numpy.roots()函数求解
solutions = np.roots([a, b, c])

或者使用`linalg.solve()`（更适合于线性方程组而非一元二次方程）：

# 这里会报错，因为这不是线性方程组，但展示原理
# solutions = np.linalg.solve([[a, b], [b, c]], [0, 1])  # 错误用法，因为不是Ax=b形式

3. **绘制函数曲线**:

from scipy.optimize import fsolve
import matplotlib.pyplot as plt

# 实际上绘制一元二次函数的图形可能需要先转换成标准形式 y = ax^2 + bx + c
def quadratic_function(x, a, b, c):
    return a * x**2 + b * x + c

# 定义x的范围和步长
x_values = np.linspace(-10, 10, 400)
y_values = quadratic_function(x_values, a, b, c)

plt.plot(x_values, y_values, label="y = " + str(a) + "x^2 + " + str(b) + "x + " + str(c))
plt.axhline(0, color='black', linestyle='--')  # 绘制y轴
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.show()