mirror of https://github.com/swig/swig
63 lines
1.2 KiB
Python
63 lines
1.2 KiB
Python
# Plot a 3D function
|
|
import gifplot
|
|
from math import *
|
|
|
|
# Here is the function to plot
|
|
def func(x,y):
|
|
return 5*cos(2*sqrt(x*x+y*y))*exp(-0.3*sqrt(x*x+y*y))
|
|
|
|
# Here are some plotting parameters
|
|
xmin = -5.0
|
|
xmax = 5.0
|
|
ymin = -5.0
|
|
ymax = 5.0
|
|
zmin = -5.0
|
|
zmax = 5.0
|
|
|
|
# Grid resolution
|
|
nxpoints = 60
|
|
nypoints = 60
|
|
|
|
cmap = gifplot.ColorMap("cmap")
|
|
frame = gifplot.FrameBuffer(500,500)
|
|
frame.clear(gifplot.BLACK)
|
|
|
|
p3 = gifplot.Plot3D(frame,xmin,ymin,zmin,xmax,ymax,zmax)
|
|
p3.lookat(2*max([xmax-xmin,ymax-ymin,zmax-zmin]))
|
|
p3.autoperspective(40)
|
|
p3.rotu(60)
|
|
p3.rotr(30)
|
|
p3.rotd(10)
|
|
|
|
def drawsolid():
|
|
p3.clear(gifplot.BLACK)
|
|
p3.start()
|
|
dx = 1.0*(xmax-xmin)/nxpoints
|
|
dy = 1.0*(ymax-ymin)/nypoints
|
|
cscale = 240.0/(zmax-zmin)
|
|
x = xmin
|
|
for i in xrange(0,nxpoints):
|
|
y = ymin
|
|
for j in xrange(0,nypoints):
|
|
z1 = func(x,y)
|
|
z2 = func(x+dx,y)
|
|
z3 = func(x+dx,y+dy)
|
|
z4 = func(x,y+dy)
|
|
c1 = cscale*(z1-zmin)
|
|
c2 = cscale*(z2-zmin)
|
|
c3 = cscale*(z3-zmin)
|
|
c4 = cscale*(z4-zmin)
|
|
c = (c1+c2+c3+c4)/4
|
|
if (c < 0) : c = 0
|
|
if c > 239 : c = 239
|
|
p3.solidquad(x,y,z1,x+dx,y,z2,x+dx,y+dy,z3,x,y+dy,z4,c+16)
|
|
y = y + dy
|
|
x = x + dx
|
|
|
|
print "Making a nice 3D plot..."
|
|
drawsolid()
|
|
|
|
frame.writeGIF(cmap,"image.gif")
|
|
print "Wrote image.gif"
|
|
|