The return value should divided by the steradian unit here, so:

```
return dh * (zp1 ** 2 * da ** 2) / self.efunc(z) / u.sr
```