par= LSquares(par)
[par,options]=leastsq('resid4',par,options,[],.0075, 2.58, .04, bond(:,2));
