[f,J]=LSMeyer(x)

function [f,J]=LSMeyer(x)
y=[34780;28610;23650;19630;16370;13720;...
    11540;9744;8261;7030;6005;...
    5147;4427;3820;3307;2872];
I=(1:16)';
t=45+5*I;w=x(2)./(t+x(3));
v=exp(w);z=x(1)*v;
f=y-x(1)*v;
J=[-v,-z./(t+x(3)),z.*w./(t+x(3))];