Numerical computation of Tricomi's psi function by the trapezoidal rule