// xcas version=0.7.3 fontsize=20 font=0
// fltk 7Fl_Tile 13 74 930 50 20 0
[
// fltk N4xcas19Multiline_Input_tabE 13 74 930 26 20 0
N:=nextprime(1000)*nextprime(2000)
,
// fltk N4xcas10Log_OutputE 13 100 930 1 20 0

,
// fltk N4xcas8EquationE 13 101 930 23 20 0
2021027
]
,
// fltk 7Fl_Tile 13 126 930 61 20 0
[
// fltk N4xcas19Multiline_Input_tabE 13 126 930 27 20 0
f(x):=x^2+1
,
// fltk N4xcas10Log_OutputE 13 153 930 1 20 0

,
// fltk N4xcas8EquationE 13 154 930 33 20 0
 (x)->x^2+1
]
,
// fltk 7Fl_Tile 13 189 930 31 20 0
[
// fltk N4xcas19Multiline_Input_tabE 13 189 930 30 20 0

,
// fltk N4xcas10Log_OutputE 13 219 930 1 20 0

]
,
// fltk 7Fl_Tile 13 222 930 369 20 0
[
// fltk N4xcas7EditeurE 13 222 930 301 20 0
263 ,
iter(a,N,f,nmax):={ // Pollard absolument pas optimise
  local y,z,j,r;
  y:=a;
  z:=a;
  for (j:=0;j<nmax;j++){
    y:=irem(f(y),N);
    z:=irem(f(f(z)),N);
    if (y==z) return "echec";
    r:=gcd(y-z,N);
    if (r!=1) return j,y,z,r;
  }
  return "j>nmax";
}:;,
// fltk N4xcas10Log_OutputE 13 523 930 45 20 0
// Parsing iter£// Success compiling iter£
,
// fltk N4xcas8EquationE 13 568 930 23 20 0
"Done"
]
,
// fltk 7Fl_Tile 13 593 930 50 20 0
[
// fltk N4xcas19Multiline_Input_tabE 13 593 930 26 20 0
iter(1,N,f,100)
,
// fltk N4xcas10Log_OutputE 13 619 930 1 20 0

,
// fltk N4xcas8EquationE 13 620 930 23 20 0
44,402447,971299,2003
]
,
// fltk 7Fl_Tile 13 645 930 36 20 0
[
// fltk N4xcas19Multiline_Input_tabE 13 645 930 35 20 0

,
// fltk N4xcas10Log_OutputE 13 680 930 1 20 0

]