int y2,y,tmp;
for (int x=0; x<p; x++)
{
  tmp = apowb(x,3)+a*x+b;
tmp = amodb (tmp,p);
 for (int i=0; i<=p-1; i++)
 {

  y2 = amodb(tmp, p);
  if (fabs(sqrt(y2))-fabs(floor(sqrt(y2)))==0 && sqrt(y2)!=tmp2)
  {
   Memo1->Lines->Add("("+FloatToStr(x)+";"+FloatToStr(sqrt(y2))+")");
    tmp2=sqrt(y2);
  }
  tmp+=(p*i);
 }
}
}