#
# disurf.py
# (c) Neil Gershenfeld 2/5/05
# plot sin(r)/r surface with dislin

from dislin import *
from numarray import *

def f(x,y):
   r = sqrt(x*x + y*y)
   return (sin(r)/r)

vec = arange(-15.5,15.5,0.1)
x = outerproduct(vec,ones(vec.shape))
y = transpose(x)
z = f(x,y)
setvar('X3VIEW',5)
setvar('Y3VIEW',2)
setvar('Z3VIEW',2)
surshade(z,vec,vec)
raw_input("CR to continue")