Ваш текущий код работает нормально, поэтому я не уверен, в чем вопрос ... Вот немного более простая реализация:
b = (1:5).*3;
A = zeros(10,10,5);
for k=1:5
[X,Y] = ndgrid(linspace(0,b(k),10), linspace(0,30,10));
A(:,:,k) = X.*exp(-Y);
end
Если вы также хотите изменить пределы y, процесс аналогичен; у вас будет две петли, и результат A
будет четырехмерной матрицей