load torqueVSdelta;

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

figure(1)
plot(delta,T_e,'k-');
axis([-1 3.5 -2.0 2.5]);
grid on;
xlabel('delta (radians)');
ylabel('T_e (pu)');
title('Torque vs. Delta','FontSize',12);