Added numerical integration of B'*C*B(works) using blas.
Added predictor/corrector time int scheme(doesn't work)