Fixed fatal theoretical error: missing c-s
Changed to eqs from Callin, proof read these and added bunch of missing c-s
Added new function: ddH_p, for second derivatives of H_p