--- trunk/esys2/doc/user/examples/helmholtztest.py	2004/12/15 07:08:39	102
+++ trunk/esys2/doc/user/examples/helmholtztest.py	2005/01/27 06:21:59	108
@@ -1,8 +1,9 @@
 # $Id$
-from mytools import *
-from esys.escript import *
-import esys.finley
+from mytools import Helmholtz
+from esys.escript import Lsup
+from esys.finley import Rectangle
 #... set some parameters ...
+kappa=1.
 omega=0.1
 eta=10.
 #... generate domain ...
@@ -11,7 +12,8 @@
 mypde=Helmholtz(mydomain)
 n=mydomain.getNormal()
 x=mydomain.getX()
-mypde.setValue(1,omega,omega*x[0],eta,n[0]+eta*x[0])
+mypde.setValue(kappa,omega,omega*x[0],eta,kappa*n[0]+eta*x[0])
 #... calculate error of the PDE solution ...
 u=mypde.getSolution()
 print "error is ",Lsup(u-x[0])
+# output should be similar to "error is 1.e-7" 
\ No newline at end of file