1- About the "digits" you saved me from a big mistake! :D
I have forgotten that "Digits" should be used!
2- About the definition of "f(r)", It is not the main problem but you are right! :)
3- About the FR procedure, it may need some improvments, but it is not the main problem. I don't know what are autocompile and hfloat! and I haven't used them.
4-About maple 18, I should try it! :)
5- Finally, the main problem: On my system when it reaches to the evalf at the beginnig of the Numerical Result section, it stops for a long time (more than 1000 seconds for N=100 !), so it seems that it has problems in numerically evaluating the integrals which are in A[i] (coefficients of SS).
6- Question: Does it stop for a noticeable time at the begining of Numerical Result section on your system?