define f(x,y) {return  re( power(x,2) - 3*power(x,2)*y + power(y,3) -5);}

define Dfx(x,y) {return  re( 2*x - 2*x*y );}

define Dfy(x,y) {return  re(  - 3*power(x,2) + 3*power(y,2) );}



dados = fopen("dados", "w");


define listagem(iniciox, fimx,inicioy,fimy, delta, precisao) 
{
 local teste = 0;   
 local dados = fopen("dados","w");
 local teste = 0.0;
 local x = iniciox;
 local y;
 while( x < fimx)
 {	 
   y = inicioy;
	while( y < fimy)
   {
		if ( abs(f(x,y)) < precisao )
		{
			fprintf(dados, "%f  %f  \n", x,y);
			## printf("%f  %f  \n", x,y);
		}
		y = y +delta;
		}
   x = x + delta;
	} 
   fclose(dados);
}


## delta = 0.01; precisao = 0.1
listagem(-3,3, -3,3,0.1, 0.01);