function c = mysvdfit(A, y)
% MYQRFIT QR least squares fit of (x,y) data.
%          Return the coefficients of the linear least squres fit
%          f(x) = c(1)*phi_(x) + c(2)*phi_2(x) + ...
    [Q, R] = myqrfact(A);
    d = Q'*y';
    c = R\d;
end

