Habe deshalb mal ganz klein angefangen und versucht die cramersche regel mal in ne Funktion zu setzen.
habe das auch soweit fertig, wenn auch viel zu kompliziert.
dummerweise bekomme ich kein Ergebnis beim ausführen. ´
Code: Alles auswählen
import numpy as np
def cr(ma,ve):
import copy
ma1=copy.deepcopy(ma)
ma2=copy.deepcopy(ma)
ma3=copy.deepcopy(ma)
def ce(mat,ve,n):
mat[0,n]=ve[0]
mat[1,n]=ve[1]
mat[2,n]=ve[2]
ce(ma1,ve,0)
ce(ma2,ve,1)
ce(ma3,ve,2)
x1=np.linalg.det(ma1)/np.linalg(ma)
x2=np.linalg.det(ma2)/np.linalg(ma)
x3=np.linalg.det(ma3)/np.linalg(ma)
print(x1,x2,x3)
widerstandsmatrix=np.matrix('55 15 5;40 15 -5;25 -5 10')
loesungsvektor=[2,0,0]
cr(widerstandsmatrix,loesungsvektor)
kann mir jemand sagen wo der fehler liegt?