cf_matrixexponential( A,b )

function [ v ] = cf_matrixexponential( A,b )  
% cf_matrixexponential
% computes v = exp(A)*b where A is a real, symmetric and negative semidefinite
% The vector b is real
% This function is much faster than expm(A)*b as it avoids the explicit
% computation of exp(A).
% Thomas Schmelzer, December 2008
    v = cf_matrixfunction(@exp,A,b,14);
end