Pas trop sûr de la syntaxe(j'ai pas matlab sous la main) mais comme ça je dirais:
function res = compute(A, B, X0, pas)
res(:, 1) = X0
for i = [2:nb]
res(:, i) = A * res(:, i - 1) + B
end
end
Par contre pour le pas je vois pas trop ce que tu veux puisque tu vas de t en t + 1...