t=0:0.001:0.1;
n=length(t);
for i=1:n
    omega(i)=0.0739*(exp(-186*t(i)))*sin(94.2*t(i));
end
plot(t,omega), grid on, xlabel('time [sec]'), ylabel('kecepatan sudut [rad/sec]')