load equalarea;

time = opd(1,:);
delta = opd(2,:);
T_e = opd(3,:);

figure(1);
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);
pause

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