function b = myridge(y, A, alpha)

%compact svd
[U D V]=svd(A,'econ');

d=diag(D);
d1=d./(d.^2+alpha^2*ones(size(A,1),1));
D1=diag(d1);
b=V*D1*U'*y;