FUNCTION SVD_SOLVE, a, b

at=transpose(a)

SVD, at, w, u, v

n = N_ELEMENTS(w)

wp = FLTARR(n,n)

FOR i=0,n-1 DO IF w(i) NE 0 THEN wp(i,i)=1./w(i)

RETURN,v # wp # (TRANSPOSE(u) # b)

END