forma0.m

function M = forma0(F, k)
%forma0
%   M = forma0(F, k)
%
%   Calcula la matriz num�rica de flujos actualizados F(k) para un k dado.
%      F(k) = f0 + f1/k + f2/k^2 + f3/k^3 + ...
%   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 = F(:,:,1);
    for c1 = 2:ancho
        M = M + F(:,:,c1) ./ (k^(c1-1));
    end
else
    M = feval(F, k);
end