clf
subplot(311),plot(time,u1),grid on
subplot(312),plot(time,u2),grid on
subplot(313),plot(time,u3),grid on