forma1.m

function M = forma1(F, k)
%forma1
%   M = forma1(F, k)
%
%   Calcula la matriz num�rica de las derivadas de los flujos actualizados
%   F(k) con respecto a k para un k dado.
%      F'(k) = - 1 f1/k^2 - 2 f2/k^3 - 3 f3/k^4 + ...
%   F puede escribirse como recetas en tiempo discreto
%      F(:,:,1) = f0, F(:,:,2) = f1, F(:,:,3) = f2, F(:,:,4) = f3, ...
%   o como recetas en tiempo continuo, definidas con inline
%      F=inline('[-280+575/k,-12;-120,-8+20/k;-280+400/k,-12]')
%   o en un script m

[filas, columnas, ancho] = size(F);
if ancho > 1
    M = zeros(filas, columnas);
    for c1 = 2:ancho
        M = M - ( F(:,:,c1) .* ((c1-1)/(k^c1)) );
    end
else
    h = 1e-5;
    M = (feval(F, k+h) - feval(F, k-h)) ./ (h+h);
end