load mechanical;

time = opd(1,:);
Tau_a = opd(2,:);
delta_w_r = opd(3,:);
delta = opd(4,:);

subplot(3,1,1);
plot(time,Tau_a,'k-');
axis([-0.1 2 -2.0 2.5]);
grid on;
xlabel('Time (seconds)');
ylabel('Torque (per unit)');
title('Accelerating Torque vs. Time','FontSize',12);

subplot(3,1,2);
plot(time,delta_w_r,'k-');
axis([-0.1 2 -0.05 0.05]);
grid on;
xlabel('Time (seconds)');
ylabel('delta_w_r (radians)');
title('delta_w_r vs. Time','FontSize',12);

subplot(3,1,3);
plot(time,delta,'k-');
axis([-0.1 2 -1.0 2.0]);
grid on;
xlabel('Time (seconds)');
ylabel('delta (radians)');
title('Rotor angle vs. Time','FontSize',12);