/tags/rel-1-3-24/SWIG/Examples/GIFPlot/Php4/shadow/runme.php4
PHP | 79 lines | 63 code | 12 blank | 4 comment | 4 complexity | 56fcbe29e40db615ba0896d213e735d8 MD5 | raw file
Possible License(s): LGPL-2.1, Cube, GPL-3.0, 0BSD, GPL-2.0
- <?
- # Plot a 3D function
- include("gifplot.php");
- # Here is the function to plot
- function 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 = new ColorMap("cmap");
- $frame = new FrameBuffer(500,500);
- $frame->clear(BLACK);
- $p3 = new Plot3D($frame,$xmin,$ymin,$zmin,$xmax,$ymax,$zmax);
- $p3->lookat(2*($zmax-$zmin));
- $p3->autoperspective(40);
- $p3->rotu(60);
- $p3->rotr(30);
- $p3->rotd(10);
- function drawsolid() {
- global $xmax;
- global $xmin;
- global $ymax;
- global $ymin;
- global $zmin;
- global $zmax;
- global $nxpoints;
- global $nypoints;
- global $p3;
- $p3->clear(BLACK);
- $p3->start();
- $dx = 1.0*($xmax-$xmin)/$nxpoints;
- $dy = 1.0*($ymax-$ymin)/$nypoints;
- $cscale = 240.0/($zmax-$zmin);
- $x = $xmin;
- for ($i = 0; $i < $nxpoints; $i++) {
- $y = $ymin;
- for ($j = 0; $j < $nypoints; $j++) {
- $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...\n";
- drawsolid();
- $frame->writeGIF($cmap,"image.gif");
- print "Wrote image.gif\n";
- ?>