MATLAB: Short code for the expression

MATLAB

h =2;
P = [1, 4, 5, 7];
a1 = (1-h)* 1 + h * P(:, 1);
a2 = (1-h) * a1 + h * P(:, 2);
a3 = (1-h) * a2 + h * P(:, 3);
a4 = (1-h) * a3 + h * P(:, 4);
disp(a1) = 1
disp(a2) = 7
disp(a3) = 3
disp(a4) = 11

Best Answer

a = num2cell( (1-h) * repmat( (1-h) * 1 + h * P(1), 1, 4) + h * P ) ;
[a1, a2, a3, a4] = deal( a{:} )