You're super close...
Just make one change to the height function so it looks like this one below:
function [h] = height(a,t,h0,v0,t0)
h=h0+v0*(t-t0)+(1/2)*a*(t-t0).^2;
end
(Notice the .^ instead of ^ to signify element-wise exponentiation. This way you can pass entire vectors in place of t.)
Then, modify the main code by remove the for-loop and instead use something like
v0=0; %Initial velocity
h0=0; %Initial height
t0=0; %Initial time
g=9.81; %The gravitational constant %Input variables
F=input('Engine force (N): ');
m=input('Mass (kg): ');
dt=input('Precision/timesteps (s): '); %Step one
t = t0: dt: 0.15;
velocity1=velocity(acceleration(F,m,g),t,v0,t0);
height1=height(acceleration(F,m,g),t,h0,v0,t0);
plot(t,height1,'b',t,velocity1,'r');
Notice that t is now a vector and that I'm only calling the plot function once, the first curve (height vs. t) is blue ('b') and the second curve (velcoity1 vs t) is red. Hopefully that helps.
