Focusing just on the problem at hand, I'd suggest
- changing line 126 to Thrust(mdotfuel+1,mach,Alt) = mdotfuel*(U8-Uo);
- adding Thrust = zeros(6,3,3); before the for loops (line 26)
- removing all plot-related commands in the loops
- adding something along the lines of
Mo=[0 1 2];
for k = 1:3
figure(k)
plot(Mo,Thrust(:,:,k),'b.-')
end
at the end.
