Made first dgdlnrp calculation - its pretty slow, some interpolation might be a good idea