Hey, as I've said, I came back to redeem myself. I totally neglected the nature of \mu and thus blundered spectacularly.
The solution to initial question "What should be the minimal friction coefficient at angle θ so that the ball won't slip?" I will mark with (*) for clarity. Now I will focus on understanding of the issue at hand.
Okay, so I will set up the whole problem as it was given in post #1, \theta begging from horizontal, initial velocity and rotation being equal to 0.
From initial system of equations: <br />
(1)\left\{\begin{matrix}<br />
mR\dot{\theta}^2=N-mg\, sin\theta\\ <br />
mR\ddot{\theta}=mg\,cos\theta-f\;\\<br />
<br />
f=\mu\,N\end{matrix}\right.<br />
One very important thing to notice is that \mu is a function of \theta,<br />
\mu=\left\{\begin{matrix}<br />
\mu_s(\theta)\\<br />
<br />
\mu_{k}=const.\end{matrix}\right.<br />
In interval between [0,\theta_{0}), the ball is slipping and \mu=\mu_{k}, from which you can find (with a bit of work and wolframalpha ) the following diff. equation:\ddot{\theta}+\mu_{k}\,\dot{\theta}^2=\frac{g}{R}(cos\theta-\mu_{k}\,sin\theta)\;,\dot{\theta}(0)=0\;,\theta(0)=0
and it's "solution" \dot{\theta}^2=\frac{2g}{R(4\mu_{k}^2+1)}(-3\mu_{k}e^{-2\mu_{k}\theta}+sin\theta\,(1-2\mu_{k}^2)+3\mu_{k}\,cos\theta)\;(2) From angle \theta_{0} ball is rolling and thus energy is conserved:\frac{1}{2}mR^2\dot{\theta}^2+\frac{1}{5}mR^2\dot\theta^2+mgR(1-sin\theta)=mgR(1-sin\theta_{0})+\frac{1}{2}mR^2\omega_{0}^2+\frac{1}{5}mR^2\omega_{0}^2 \Rightarrow \dot{\theta}^2=\frac{10g}{7R}(sin\theta-sin\theta_{0})+\omega_{0}^2\;(3)
From equation (1) and (3) you can find \mu_{s} for curiosity sake, notice how it depends on angle:\mu_{s}(\theta)=\frac{2\,cos\theta}{17\,sin\theta-10\,sin\theta_{0}+\frac{7R}{g}\omega_{0}^2}\;(*)
This really depressed me when I saw velocity dependence (which was against my assumption in previous post).
To find \omega_{0}^2 you simply need to plug in \theta_{0} in equation (2).
Friction "coefficient" really is the culprit of the problem. An idea came to my mind to take the derivate of (2) and (3) equate them and evaluate them at \theta_{0}. However, finding \mu_{k} explicitly is not possible. Maybe you could use Taylor expansion, but playing around with
https://www.desmos.com/calculator ensured me that such \mu_{k} really does exist.
Cheers!