load EA;

time = opd(1,:);
gamma = opd(2,:);

%*********************************************************************
% Find the maximum swing of gamma
number = length(time);
maxgammaswing = max(gamma);

%*********************************************************************
% Find the value of gamma at t_fault
n=1;
while(time(n)< t_fault-0.01)
    n=n+1;
end
gamma0 = gamma(n);

%*********************************************************************
% Find the value of gamma at t_cr
n=1;
while(time(n)< t_fault+t_clear)
    n=n+1;
end
gamma_cr = gamma(n);

fprintf('\n\n t_fault = %5.4f \n', t_fault);
fprintf(' t_clear = %5.4f \n', t_clear);
fprintf(' gamma0 = %5.4f degrees \n', gamma0);
fprintf(' gamma_cr = %5.4f degrees \n', gamma_cr);
fprintf(' max. gamma swing = %5.4f degrees \n', maxgammaswing);


figure(1)
plot(time,gamma,'b-');
axis([10 11 0 180]);
grid on;
xlabel('Time (seconds)');
ylabel('Gamma (degrees)');
title('Angle vs. Time','FontSize',12);

OPF = fopen('file1.dat','w');

for n=1:number
    fprintf(OPF,'%5.3f %5.3f \n', time(n), gamma(n));
end;

fclose(OPF);