perturbacionodefile.m

function dVdh = perturbacionodefile(h, V, flag, F, G, problema)
%perturbacionodefile
%   dVdh = perturbacionodefile(h, V, flag, F, G, problema)
%
%   Derivada del factor, intensidades y valores ante la perturbaci�n.
%      h          punto de integraci�n
%      V          variables V = [k; X'; Y]
%      flag       par�metro sin uso incluido para compatibilidad con Matlab
%      F          procesos de producci�n
%      G          perturbaci�n
%      problema   para autovalores 0, para VN 1, y para TE 2
%      dVdh       derivada de las variables dVdh = [dkdh; dXdh'; dYdh]
%
%   Se usa la sintaxis de Matlab; v�ase help odefile. El par�metro flag 
%   est� incluido para garantizar la compatibilidad con Matlab y puede 
%   definirse como []. V�ase tambi�n perturbacionode.m

%Necesita perturbacion.m 

[filas, columnas, ancho] = size(F);                             %determina el tama�o del problema
k = V(1);                                                       %extrae k
X = V(2:filas+1)';                                              %extrae X
Y = V(filas+2:filas+columnas+1);                                %extrae Y
[sk, sX, sY] = perturbacion(F + h.*G, G, k, X, Y, problema, 1); %calcula las series para el grado 1
dVdh = [sk(1); sX(:,1); sY(:,1)];                               %extrae dkdh, dXdh, dYdh y las almacena