fixed a bug in doperator. Opeartor D() was not returning zero
in some cases. Instead, it was returning D(...). Example:
    in this case: D(x)*(L-x)/x it was returning: -L**(-1) + D(x)
    instead of: -L**(-1)
    NOW IT IS FIXED!