# cosconst = rootof(x = cos(x))

cosconst=\
.7390851332151606416553120876738734040134117589007574649656806357732\
84654883547594599376106931766531849801246643987163027714903691308420\
31578044057462077868852490389153928943884509523480133563127677223158\
09563537765724512043734199364335125384097800343406467004794021434780\
80271801883771136138204206631633503727799169673122323006138865820362\
17708109978970626842405880948986832618606004858989585487257367640150\
75227608180391459518101628159120096461646067544051326415171064466281\
10936082584878371383955556175141494715939006277527563258634938869730\
14083665152511520426788515302529417180365176420177086071899276016098\
74327154552267565798246297611775539616699549311158566534834953838523\
15963602527499558725250666640131318740139253888805520618698592139252\
52854154110791002998282929864052169046554736696871438735646006521225\
46891499759209699758501364249508565047324972584248371554836483437275\
83746752545335800664200478839718858489014531155060417812337047773953\
4717103451195854600726561464721419787537388023680

/*
  scale=1000;pi_2=a(1)*2
  for(scale=10;scale<=1000;scale+=5){
    eps=A^(1-scale);p=pi_2/1
    oc=0;while(abs(oc-c)>eps){
      oc=c;s=s(c+p);c=(c+s+s)/3
    }
    c
  }
*/