  {We provide a probabilistic approach in order to investigate the smoothness
of the solution to the Poisson and Dirichlet problems in $L$-shaped domains. In
particular, we obtain (probabilistic) integral representations for the
solution. We also recover Grisvard's classic result on the angle-dependent
breakdown of the regularity of the solution measured in a Besov scale.
