python用numpy解二元一次方程

#python用numpy解二元一次方程

#2025-11-05 win10+python3.9 调试 成功..

#7*x+12*y==10
#3*x-2*y==5

#a1*x + a2*y = c1
#b1*x + b2*y = c2

import numpy as np
# 系数矩阵和常数矩阵
A = np.array([[7, 12], [3, -2]])
b = np.array([10, 5])
# 高斯消元法求解
x = np.linalg.solve(A, b)
print(x)

#解是 [ 1.6 -0.1]

x1=round(x[0],2)
y1=round(x[1],2)
print (x1)
print (y1)

print (5*x1+5*y1)

#下面 是 已知 解 的 范围 后，
#[ 1.6 -0.1]
#x=1.6
#y=-0.1
#在 这个 范围 内 找解。
#x 在 1 到 5 ，y 在 -3 到 5 中 找解。

import numpy as np

for x in np.arange(1, 5.0, 0.1):
	for y in np.arange(-3, 5.0, 0.1):
		x=round(x,2)
		y=round(y,2)
		#print (x)
		#print (y)
		
		#if 7*x+12*y==10:
			#print (x)
			#print (y)
			#print ('====')
			#print (5*x+5*y)
			#print (3*x-2*y)
			#print (7*x+12*y)
		if int(3*x-2*y)==5 and 7*x+12*y==10:
			print ('here..')
			print (x)
			print (y)
	
	
	
'''
steps = [i * 0.1 for i in range(1, 10)]
for i in steps:
	print(i)
'''
print ('aaa')