% Método de Euler para Resolver una Ecuación Diferencial Ordinaria
% Resuelve por el Método de Euler Simple Hacia Adelante
% La ecuación de oscilación en el problema de estabilidad transitoria
% FALLA SIN POTENCIA TRANSMITIDA
% Autor: Francisco M. González-Longatt
% Fecha: March 9, 2006
%
% Precaución: Solamente para usos de enseñanza.
t0=0;                 % Tiempo Inicial
t(1)=t0; 
tn=1.5;               % Tiempo Final de simulación
tc=0.3;               % Tiempo de despeje de la falla 
n=10;                 % Subintervalos de integración
delta_t=(tn-t0)/n     % Paso de Integración
delta0=asin(0.8/1.8); % Condición Inicial de angulo
vel0=0;               % Velocidad Inicial de la máquina
H=9.94;               % Constante de Inercia en segundos
Pm=0.8;               % Potencia mecánica de entrada INICIAL
f=60;                 % Frecuencia
Pmax1=1.8;            % Potencia máxima de la máquina ANTES
Pmax2=0;              % Potencia máxima de la maquina DURANTE
Pmax3=1.8;            % Potencia máxima de la maquina DESPUÉS
t(1)=t0; velocidad(1)=vel0;
delta(1)=delta0;
for i = 1:1:n
    t(i+1)=t(i)+delta_t;
    delta(i+1)=delta(i)+delta_t*velocidad(i);
    if t(i+1)<tc
        Pmax=Pmax2;
    end
    if (t(i+1)>=tc)  
        Pmax=Pmax3;
    end
     velocidad(i+1)=velocidad(i)+(pi*f/H*(Pm-Pmax*sin(delta(i))))*delta_t;
end