Numerical methods for pricing occupation time derivatives