平方根法
开始

分别输入A矩阵和B矩阵的元素
aij(i,j=1,2,...,n)和bi(i-1,2,...,n)

对A进行LL┬分解计算，计算L的
第一列元素，
l11=√a11,li1=ai1/l11,i=2,3,...,n
对j=2,3,...,n做

ljj=(ajj-2^11jkjkl);lij=(aij-11jkjkikll)/ljj,
i=j+1,j+2,...,n且j不等于n

求解方程组Ly=b y1=b1/lii,yi=(bi-11ijjjyli)/lii,i=2,3,...,n
求解方程组L┬x=y
Xn=yn/lnn,xi=(yi-nijjjixl1)/lii,i=n-1,n-2,...,1

输出方程组的解x=(x1,x2,...,xn)┬
结束