平方根法的matlab实现

function [ lyx ] = untitled5（ ab ）
y=[];x=[];
n=size（a1）;
for j=1:n
    c=0;
    for k=1:j-1
        c=c+l（jk）^2;
    end
    l（jj）=sqrt（a（jj）-c）;
    for i=j+1:n
        if j==1
            l（ij）=a（ij）/l（jj）;
        else
            c=0;
            for k=1:j-1
                c=c+l（ik）*l（jk）;
            end
            l（ij）=（a（ij）-c）/l（jj）;
        end