subplot(311), plot(time, ref), grid on, ylabel('reference INPUT')
subplot(312), plot(time, omega), grid on, ylabel('omega [RPM]')
subplot(313), plot(time, teta), grid on, ylabel('teta [degree], xlabel[time [sec]/')

