MATLAB: Vectorisation of a simple for loop

cumsumfor loopspeedvectorisation

Will it be possible to write this code vectorised?
x=ones(10,1);
a=zeros(size(x));
b=zeros(size(x));
c=zeros(size(x));
for i = 1:size(x,1)-1
    a(i+1) = a(i)+x(i)+b(i);
    b(i+1) = b(i)+a(i)+c(i);
    c(i+1) = c(i)+b(i)+c(i);
end
Any help would be greatly appreciated.

Best Answer

 A = [1 1 0 0 ; 1 1 1 0 ; 0 1 2 0; 1 0 0 1];
 Q = cell2mat(arrayfun(@(i) A^i, 1:9, 'Uni', false));
 Q = [zeros(1, 4); reshape(Q(end,:), 4, [])'];
 a = Q(:,1); b = Q(:,2); c = Q(:,3);
Related Question