Efficient numerical approximation of a non-regular Fokker–Planck equation associated with first-passage time distributions