Singular value decomposition

Hi.
I have added routines to compute the singular value decomposition of a matrix. Also included is a routine to compute gaussian quadrature rules for various weight functions.

Greetings
  Timo