forma2.m

function M = forma2(F, k)
%forma2
%   M = forma2(F, k)
%
%   Calcula la matriz num�rica de la segundas derivadas de los flujos
%   actualizados F(k) con respecto a k para un k dado.
%      F"(k) = 1*2*f1/k^3 + 2*3*f2/k^4 + 3*4*f3/k^5 + ...
%   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)*c1/(k^(c1+1))) );
    end
else
    h = 1e-5;
    M = (feval(F, k+h) - 2.*feval(F, k) + feval(F, k-h)) ./ (h^2);
end