(MODE R = LONG LONG REAL; PROC (R)R rt = long long sqrt;
  R r17 = rt(17), a = rt(34 - 2*r17);
  print((
  (-1/16 + r17/16 + a/16 + rt(17 + 3*r17 - a - 2*rt(34 + 2*r17))/8),
  newline, long long cos(2*long long pi/17))))