Matlab. Решение СЛАУ итерационными методами.

[MyQwErTy]

Задача.
1. С помощью программы сгенерировать матрицу A и вектор правых частей b для СЛАУ вида Ax=b.
2. Проверить условия сходимости итерационных методов для матрицы A. При необходимости преобразовать матрицу A к виду, позволяющему вычислить решения.
3. Вычислить решения с помощью метода простой итерации и Гаусса-Зейделя, количество итерации n=20. Построить графики зависимости xi от n. Вычислить вектор невязки решения и ее норму.
4. Получить точное решение СЛАУ, вычислить вектор погрешности решений по методам простой итерации и Гаусса-Зейделя и их нормы.
5. Оценить относительную погрешность решения СЛАУ, если ∆A определяется десятипроцентным увеличением диагональных элементов, а ∆b таким же уменьшением всех элементов.


Решение

n1=input('Введите номер группы: ');
n2=input('Введите ваш номер в списке группы: ');
n=n1+n2;
rand('seed',n);
A=rand(7);
b=10*rand(7,1);
Au=triu(A);
Al=tril(A);
Ad=diag(diag(A));
su=sum(sum(abs(Au)));
sl=sum(sum(abs(Al)));
Su=sum(abs(Au));
Sl=sum(abs(Al));
sigma=5;
O=ones(7);
Ou=triu(O);
S6=diag(sigma*(Su./Sl));
At=sigma*((Al+Ou)./(Au'+Ou));
neo=rem(n2,2);
mode=rem(n2,3);
if (mode==0) Al=(sigma*su/sl)*Al; end
if (mode==1) Al=S6*Al; end
if (mode==2) Al=At.*Al; end
A=Al+Ad+Au;
ma=max(max(abs(A)));
if (ma>1000) A=0.08*A; end
A=0.5*A+2.75*diag(diag(A)), b
%простая итерация
x=b;
v=x;
m=20;
d=diag(A);
S=diag(d);
Si=inv(A);
T=S-A;
for i=1:m
    y=(T*x+b)/d;
    B=T*Si;
    mis=(norm(B))*m*(norm(x-y));   <--------------------------  Тут ошибка
    x=y;
    v=[vx];
end
%normaB=norm(B)
x
mis
figure(1)
plot(v);
grid
pause
%Зейдель
x=b;
v=x;
S=tril(A);
Si=inv(A);
T=S-A;
for i=1:m
    y1=Si*(T*x+b);
    B=T*Si;
    Bu=triu(B);
    Bl=tril(B);
    mis=(norm(y1-x))*(norm(Bu)/(1-(norm(Bl))));
    x=y1;
    v=[vx];
end
normaB=norm(B)
x
mis
figure(2)
plot(v);
grid

Косяк
Строка 39, а именно:

??? Error using ==> minus
Matrix dimensions must agree.

Error in ==> Untitled at 39
mis=(norm(B))*m*(norm(x-y));

Просьба помочь с исправлением ошибки. Заранее спасибо.

[MyQwErTy]

____________up