program Matlab gerak harmonik sederhana pada gerak jatuh bebas

%gerak jatuh bebas
yo=input('Masukkan ketinggian awal yo :');
vo=input('Masukkan kecepatan awal vo :');
h=0.001; %step size
g=10; %percepatan grafitasi
y=yo; %inisialisasi untuk y
v=vo; %inisialsisasi untuk
i=1;
fid=fopen('jatuh.txt','w');
while(y>0)
t=i*h;
v=v-g*h
 y=y+v*h;
if (y<0)
break;
end;
temp=y;

fprintf('%f %f %f\n',t,abs(v),y);
fprintf(fid,'%f %f %f\n',t,abs(v),y);
i=i+1;
end
fclose(fid)
load jatuh.txt;
t=jatuh(:,1);
y=jatuh(:,2);
v=jatuh(:,3);
plot(t,y,t,v,'lineWidth',2.5);
xlabel('t');
ylabel('y / v');
legend('y','v');

hasil Run nya :