RNPL example
########################################################################
# RNPL source for axisymmetric wave-equation in cylindrical coordinates
# (rho,z). Second order form and single uniform mesh (including (0,0)
# used. Radiation conditions imposed on outer-rho and both z-boundaries.
# Regularity imposed at rho=0. Scheme seems to be stable
# phi_tt = phi_zz + phi_rhorho + phi_rho / rho
# Copyright 1996, Matthew W Choptuik, The University of Texas at Austin
########################################################################
########################################################################
# Parameters, grid and grid functions
########################################################################
Wave equation in 2D, cylindrical coordinates: