use blib;
use Carp;

# $SIG{__DIE__} = sub {die Carp::longmess(@_);};

use PDL;
use PDL::Graphics::TriD;
use PDL::IO::Pic;
use PDL::Graphics::TriD::Polygonize;


$orig = PDL->pdl(0,0,0)->float;

sub func1 {
	my($x,$y,$z) = map {$_[0]->slice("($_)")} 0..2;
	$r = $x**2 + 1.5*$y**2 + 0.3 * $z**2 + 5*($x**2-$y)**2;
	$res = ($r - 1) *  -1;
#	print $res;
	return $res;
}

$a = PDL::Graphics::TriD::StupidPolygonize::stupidpolygonize($orig,
	5, 50, 10,\&func1)  ;

# print $a;
imag3d $a,{Lines => 0, Smooth => 1};