<< Numerico.m
<< Numerical.m

AL = {{2,1},{4,2.1}}

AC = {{2,1},{1,2}}

b  = {1.222,3.666}

ALFactor[AL,b,2,3,"u","cp","dp"]

ALCholeski[AC,b,2,3,"dp"]

ALFactorwithRefinement[AL,b,2,3,"u","sp","ni",0.01,"dp","dr"]

ALCholeskiwithRefinement[AC,b,2,3,"n1",0.01,"sp","dr"]

ALFactorwithnRefinements[AL,b,2,3,"u","sp",3,"dp","sr"]

ALCholeskiwithnRefinements[AC,b,2,3,3,"sp","sr"]