g=9.81; c=12.5;
m=150;
t=[0:0.1:120];
v=(g*m/c)*(1-exp(-(c/m)*t));
plot(t,v)