blabla

rec = 1;

aa=[];
bb=[];
r=[];
ind=[];

for a =[2:10000],
  b = round( a^(e/pi) );
  d = a^(1/pi) - b^(1/e);
  aa=[aa;a];
  bb=[bb;b];
  r=[r;d];
  if (abs(d)<rec),
    ind = [ind,a-1];
    rec = abs(d);
  endif
endfor

[aa(ind),bb(ind),r(ind)]
plot(aa,r,'*')