Is there a smooth function on an interval in $\mathbb R$, not analytic on any subinterval, whose Taylor series at every point has positive radius of convergence? The Fabius function might be an example, but this is questionable since the n'th derivative has maximum $2^{\sigma(n)}$, where $\sigma(n)=\frac{n(n+1)}{2}$, which is not quite good enough using a crude estimate for radius of convergence if there are points where many derivatives are close to the maximum.

We define $$\Psi(x)= \sum_{k\ge 0} 2^{k}\psi_{\sigma_k}(xx_k),\quad \psi_{\sigma}(y)=\exp{{\vert x\vert}^{\frac{1}{s1}}}, $$ where $(x_k)_{k\ge 0}$ is dense in $\mathbb R^d$ and $(\sigma_k)_{k\ge 0}$ is decreasing and valued in $[s_1,s_0]\subset(1,+\infty)$. That function is good explicit substitute to Fabius function since it is smooth and nowhere analytic: even better, it is multidimensional and its analytic wavefrontset is all the cotangent space (minus the zero section). To prove this use Gevrey classes. It seems likely that the radius of convergence of the Taylor series is positive on a dense subset of $\mathbb R^d$. Anyhow it is a good candidate. Bazin. 

