#! /usr/bin/python # # #-------------------------------------------------------------------- def main (): nn = 48 print "#VRML V2.0 utf8" print "Shape {" print "geometry ElevationGrid {" print "xDimension",nn print "zDimension",nn print "height [" proc02 (nn) print "]" print " solid FALSE" print "}" print "}" # #-------------------------------------------------------------------- def proc02 (nn): import math # center = nn / 2.0 - 0.5 for ix in range (nn): t1 = center - ix t2 = t1 * t1 for iy in range (nn): t3 = center - iy t4 = t3 * t3 rr = math.sqrt (t2 + t4) zz = 20.0 - rr print int(zz), print "\n#" #-------------------------------------------------------------------- main () #--------------------------------------------------------------------