% 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
% CAMBIO DE POTENCIA
% 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
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 ángulo
vel0=0;               % Velocidad Inicial de la máquina
H=9.94;               % Constante de Inercia en segundos
Pm1=0.9;              % Potencia mecánica de entrada FINAL
f=60;                 % Frecuencia
Pmax1=1.8;            % Potencia máxima de la máquina
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);
    Pmax=Pmax1;
    velocidad(i+1)=velocidad(i)+(pi*f/H*(Pm1-Pmax*sin(delta(i))))*delta_t;
end