Breaking-up dr1.m

Currently the function dr1.m contains most of the code for computing first, second and third order approximation of stochastic model. Current shortcomings:

Proposed re-organization of the code

Factorized code

Solution code

Storing results