Python program to solve an equation using Gauss elimination,equationgauss,''' x = gaus


''' x = gaussElimin(a,b).    Solves [a]{b} = {x} by Gauss elimination.'''from numpy import dotdef gaussElimin(a,b):    n = len(b)  # Elimination Phase    for k in range(0,n-1):        for i in range(k+1,n):           if a[i,k] != 0.0:               lam = a [i,k]/a[k,k]               a[i,k+1:n] = a[i,k+1:n] - lam*a[k,k+1:n]               b[i] = b[i] - lam*b[k]  # Back substitution    for k in range(n-1,-1,-1):        b[k] = (b[k] - dot(a[k,k+1:n],b[k+1:n]))/a[k,k]    return b

评论关闭