Module HSTPLR in SLATECSolve the standard five-point finite difference approximation on a staggered grid to the Helmholtz equation in polar coordinates.
CALL HSTPLR (A, B, M, MBDCND, BDA, BDB, C, D, N, NBDCND, BDC, BDD, ELMBDA, F, IDIMF, PERTRB, IERROR, W)| I2b1a1a | Poisson (Laplace) or Helmholtz equation on a rectangular domain (or topologically rectangular in the coordinate system) |