Code666 (代码贴、代码片段)

[Python] Python program to solve an equation using Gauss el →→→→→进入此内容的聊天室

来自 , 2021-03-21, 写在 Python, 查看 154 次.
URL http://www.code666.cn/view/f01acd85
  1. ''' x = gaussElimin(a,b).
  2.    Solves [a]{b} = {x} by Gauss elimination.
  3. '''
  4. from numpy import dot
  5.  
  6. def gaussElimin(a,b):
  7.     n = len(b)
  8.   # Elimination Phase
  9.     for k in range(0,n-1):
  10.         for i in range(k+1,n):
  11.            if a[i,k] != 0.0:
  12.                lam = a [i,k]/a[k,k]
  13.                a[i,k+1:n] = a[i,k+1:n] - lam*a[k,k+1:n]
  14.                b[i] = b[i] - lam*b[k]
  15.   # Back substitution
  16.     for k in range(n-1,-1,-1):
  17.         b[k] = (b[k] - dot(a[k,k+1:n],b[k+1:n]))/a[k,k]
  18.     return b
  19. #//python/7428

回复 "Python program to solve an equation using Gauss el"

这儿你可以回复上面这条便签

captcha

贴代码片段 - Code666