Lineare Gleichungen lösen mit numpy
Hier ein kurzes Beispiel aus der numpy-Dokumentation, wie man mit Hilfe von numpy lineare Gleichungssysteme lösen kann:
Zu lösen sind folgende Gleichungen:
- 3 * x0 + 1 * x1 = 9
- 1 * x0 + 2 * x1 = 8
Die Koeffizienten kommen in die entsprechenden numpy-Arrays, dann ruft man linalg.solve
auf:
import numpy as np a = np.array([[3,1], [1,2]]) b = np.array([9,8]) x = np.linalg.solve(a, b) print(x) # gibt [ 2. 3.] |
Den Plot habe ich mit LaTeX erstellt, siehe http://uweziegenhagen.de/?p=3516.