function [alpha,beta,chisq] = mrqcof(x,y,sig,ndata,a,ia,ma,nalp,funcs,model)
%
% Used by mrqmin to evaluate the linearized fitting matrix alpha,
% and vector beta as in 15.5.8 (Numerical Recipies), and calculate X