-------先写function ,和后面的是分开的 .m file
function a = GS (A,b,U,L,D,x)
a = ......... <--- formula 书上有
end
------- main .m file
A = ...
b = ...
U = ...
L = ...
D = ...
err = zeros(1,20) <--------初始 设定err为一个有20个0组成的vector
for n = 1:20 <---自己随便设
N1 = x
x = GS (A,b,U,L,D,x)
N2 =x
err = norm(N1-N2) <-------- 这个算err的formula 书上应该有
end