hello, how to solve the problem for matrix that a is diagonal matrix and u is random column matrix. I would like to solve the equation to get u=a*u inside the loop. But , I got index in the position 2 is invalid answer. Could you please help me?

e=ones(10,1);

a=spdiags([-e 2*e -e],-1:1,10,10);

A=full(a);

u=rand(10,1);

for k=0:9

u(:,k+1)=(A*u(:,k))/norm(A*u(:,k));

end

Answer 1

[Alan Stevens]

Matlab's indexing starts at 1 not zero, so it is complaining at u(:,k) when k is zero.

[Thin Rupar Win]

Thank you very much Alan Stevens. Have a nice day.