Fixed typo in exponent of qb_star in engelund_hansen function. Changed ^(5.2) to **(5./2.)