include(`../test/calculus.txt'):

let(pi,4*atan(1)):

let(theta,smooth(<
for(<k,0,100,10>,{
sub({let(a,int([k])):
print_newtonsolve($({1 - cos([x]) - [a]/100}), .5 , 15 ),
})
})>)):

print_[theta],\n:

print_for(<k,1,len([theta]),1>,{
sub({
let(ang,elem([theta],[k])):
print_([ang] - 0.5 * sin(2*[ang]))/[pi],\n
})
})