%高斯賽德爾迭代
clear?all;
close?all;
clc;
tic
format?longe?
disp（‘請輸入參數‘）;
K=input（‘維數K=‘）;
tic
A=100*rand（K）;%?A元素是0-100
for?i=1:K
????A（ii）=sum（abs（A（i:）））+25*rand（1）;?%對角占優的量為0~25
end
b=zeros（K1）;
for?i=1:K;
????x=0;
????for?r=1:K;
????x=x+A（ir）;
????end
????b（i1）=x;
end???%產生b矩陣，b中的元素為A中對應行的和，目的是使方程解全為?1
jd=input（‘控制精度j