function[alpha_mrq,beta_mrq,chisq] = mrqcof(funcs,x,y,sig,ndata,ma,a);

g=1.0./sig.^2;

chisq=0;

[ymod,dyda] = funcs(x,a,ma,ndata);

dy=y-ymod;

gdyda=diag(g)*dyda';
beta_mrq=dy*gdyda;
alpha_mrq=dyda*gdyda;
chisq=dy.^2*g';
% for k=1:ndata
%   for i=1:ma
%     for j=1:ma
%       alpha(i,j)=alpha(i,j)+ g(k)*dyda(i,k)*dyda(j,k);
%     end
%     beta(i)=beta(i)+g(k)*dy(k)*dyda(i,k);
%   end
%   chisq=chisq+g(k)*dy(k)^2;
% end
