Algorithm 572: Solution of the Helmholtz Equation for the Dirichlet Problem on General Bounded Three-Dimensional Regions [D3]
H e r e ~2, a t h r e e d i m e n s i o n a l b o u n d e d reg ion , c, a n a r b i t r a r y r e a l c o n s t a n t (posit ive , nega t ive , o r zero) , a n d t h e f u n c t i o n s f a n d g a r e spec i f i ed b y t h e user . T h e L a p l a c e o p e r a t o r A is in C a r t e s i a n c o o r d i n a t e s . A s e c o n d o r d e r a c c u r a t e f i n i t e -d i f f e rence m e t h o d is u s e d to d i s c r e t i ze t h e H e l m h o l t z e q u a t i o n . T h e r e s u l t i n g l i n e a r s y s t e m of e q u a t i o n s is r e d u c e d to a c a p a c i t a n c e m a t r i x e q u a t i o n t h a t is so lved a p p r o x i m a t e l y b y a c o n j u g a t e g r a d i e n t