I know $\sum\_{k=0}^{n} \sin(k)$ is bounded by a constant
 and $\sum\_{k=0}^{n} \sin(k^2)$ is [not bounded](https://mathoverflow.net/questions/201250/is-sum-k-1n-sink2-bounded-by-a-constant-m) by a constant. 
Then, what about $\sum\_{k=0}^{n} (|\sin(k)|-2/\pi)$?


From numerical calculation, $\max\_{n=0...10^8}(\sum\_{k=0}^{n} (|\sin(k)|-2/\pi))=0.0900478$ which is much small compared to $\max\_{n=0...10^8}\sum\_{k=0}^{n} \sin(k^2)=1882.86$.


So, I suppose $\sum\_{k=0}^{n} (|\sin(k)|-2/\pi)$ can be bounded by a constant, but I don't know how to prove it.



Is $\sum_{k=0}^n (|\sin(k)|-2/\pi) $ bounded by a constant $M$?
As a partial answer, we can show that if $S\_N=\sum\_{k=0}^N (|\sin(k)|-2/\pi)$, then:


$S\_N=O(1)-\frac{4}{\pi}\sum\_{k=2}^{N}\sum\_{m=1}^{[k \log^2 k]}\frac{\cos 2mk}{4m^2-1}=O(\sum\_{m=1}^{[N\log^2 N]}{\frac{\min(N\log^2N, ||\frac{m}{\pi}||^{-1})}{m^2}})$, 


where as usual $||x||$ represents the distance to the closest integer. 


In particular if we know that $||\frac{m}{\pi}||$ is not too "small" (e.g it is $>>m^{-1+\epsilon}$) for enough $m$ we may be able to show boundness or at least some better estimate than the uniform distribution $o(N)$ noted by @Wojowu above. Conversely, if there are infinitely many $m$ with $||\frac{m}{\pi}|| < C\frac{1}{m^2}$, the estimate doesn't work at least as boundness goes, though a cleverer estimate of the double trigonometric sum above could work.


Proof: We use the Fourier series $|\sin k|=\frac{8}{\pi}\sum\_{m=1}^{\infty}\frac{\sin^2 mk}{4m^2-1}$ which converges absolutely and we cut it at $[k\log^2 k]$ for $k \ge 2$ as the terms for $k=0,1$ in $S\_N$ are absorbed in the $O(1)$ above, so we won't bother with them. 


Then the tail of the Fourier series for $|\sin k|$ is obviously bounded by $C\sum\_{m > k\log^2 k}\frac{1}{m^2}=O(\frac{1}{k\log^2 k})$, hence summing on $k \ge 2$ we get the tail $O(1)$.


For the rest, using $2\sin^2(mk)=1-\cos(2mk)$ and $2\frac{1}{4m^2-1}=\frac{1}{2m-1}-\frac{1}{2m+1}$, we get that $\frac{8}{\pi}\sum\_{m=1}^{[k \log^2 k]}\frac{\sin^2 mk}{4m^2-1}=\frac{2}{\pi}+O(\frac{1}{k\log^2 k})-\frac{4}{\pi}\sum\_{m=1}^{[k \log^2 k]}\frac{\cos 2mk}{4m^2-1}$


Putting the above together we get the first statement above, namely that:


$S\_N=O(1)-\frac{4}{\pi}\sum\_{k=2}^{N}\sum\_{m=1}^{[k \log^2 k]}\frac{\cos 2mk}{4m^2-1}$


Rearranging the terms with the same $1 \le m \le [N\log^2 N]$ and using the well known (and easy to prove) estimate:


$\sum\_{j=m}^{n}\cos(2\pi jx)=O(\min(n-m+1, ||x||^{-1}))$, we get our claimed estimate result:


$\sum\_{k=2}^{N}\sum\_{m=1}^{[k \log^2 k]}\frac{\cos 2mk}{4m^2-1}=O(\sum\_{m=1}^{[N\log^2 N]}{\frac{\min(N\log^2N, ||\frac{m}{\pi}||^{-1})}{m^2}})$


since for fixed $m$, the sum in $k$ has at most $N \log^2N$ terms and the corresponding $x$ is of course $\frac{m}{\pi}$


