Add Cholesky decomposition function in storm::numeric namespace.

Also adds a test for 3x3 and 4x4 decompositions with a tolerance of 1e-4