t=0:0.001:0.6;
n=length(t);
for i=1:n
    omega(i)=0.521*exp(-361.15*t(i)) - 17.188*exp(-10.95*t(i)) + 16.6667;
end
plot(t,omega), grid on, xlabel('time [sec]'), ylabel('kecepatan sudut [rad/sec