clear;
t=0:0.001:0.3;
n=length(t);
for i=1:n
    omega(i)=-1.52*(10^5)*exp(-3.992*t(i)) + 1.65*(10^5);
end
plot(t,omega), grid on, xlabel('time [sec]'), ylabel('kecepatan sudut [rad/sec]')