Я пытаюсь написать код, который разлагает матрицу A для решения для вектора b за 50 шагов. Мне дается вектор b0, и оттуда я должен go, пока я не достигну вектора b50. Каждая итерация должна перезаписывать старый вектор новым. Пока это мой код:
%find lower and upper of the A matrix Ax=b
% L U X = b
% d = b/L
% U X = d
% L d = b
[L U] = lu(A);
b0 = b;
d0 = L\b0;
x0 = U\d0
%%%%%
b1 = x0
d1 = L\b1;
x1 = U\d1
%%%%
b2 = x1
d2 = L\b2;
x2 = U\d2
%%%% and so on.....
, но я хочу продолжить это до b50
Как бы я мог записать это в a для l oop для пятидесяти итераций? Я знаю, что это может быть базовый c вопрос, но я новичок.
Спасибо