Singular boundary method for modified Helmholtz equations