Solves a*x = b for x, using LU decomposition.
ludecomp(a,n,zero=0,one=1) Instance Public methods Performs LU decomposition
lusolve(a,b,ps,zero=0.0) Instance Public methods Solves a*x = b for x, using