jacexpfit.m

A sparse variant of the Levenberg-Marquardt algorithm implemented by levmar has been applied to bund

function jac = jacexpfit(p, data)
  n=data;
  m=max(size(p));

  for i=1:n
    jac(i, 1:m)=[exp(-p(2)*i), -p(1)*i*exp(-p(2)*i), 1.0];
  end