--- a/sfepy/examples/diffusion/poisson_parametric_study.py
+++ b/sfepy/examples/diffusion/poisson_parametric_study.py
@@ -154,7 +154,7 @@
     return out
 
 def vary_omega1_size( problem ):
-    """Vary size of \Omega1. Saves also the regions into options['output_dir'].
+    r"""Vary size of \Omega1. Saves also the regions into options['output_dir'].
 
     Input:
       problem: Problem instance
@@ -205,7 +205,7 @@
         filename = join( output_dir,
                          ('log_%s.txt' % d_format) % ii )
         fd = open( filename, 'w' )
-        log_item = '$r(\Omega_1)$: %f\n' % diameter
+        log_item = r'$r(\Omega_1)$: %f\n' % diameter
         fd.write( log_item )
         fd.write( 'solution:\n' )
         nm.savetxt(fd, state())
--- a/sfepy/homogenization/coefs_phononic.py
+++ b/sfepy/homogenization/coefs_phononic.py
@@ -283,7 +283,7 @@
     return slogs, gaps, kinds
 
 def get_callback(mass, solver_kind, mtx_b=None, mode='trace'):
-    """
+    r"""
     Return callback to solve band gaps or dispersion eigenproblem P.
 
     Notes
@@ -326,7 +326,7 @@
     return eval(mode + '_callback')
 
 def find_zero(f0, f1, callback, freq_eps, zero_eps, mode):
-    """
+    r"""
     For f \in ]f0, f1[ find frequency f for which either the smallest (`mode` =
     0) or the largest (`mode` = 1) eigenvalue of problem P given by `callback`
     is zero.
--- a/sfepy/scripts/probe.py
+++ b/sfepy/scripts/probe.py
@@ -174,7 +174,7 @@
                 output('data ->', os.path.normpath(txt_filename))
 
 def integrate_along_line(x, y, is_radial=False):
-    """
+    r"""
     Integrate numerically (trapezoidal rule) a function :math:`y=y(x)`.
 
     If is_radial is True, multiply each :math:`y` by :math:`4 \pi x^2`.
--- a/sfepy/solvers/nls.py
+++ b/sfepy/solvers/nls.py
@@ -439,7 +439,7 @@
             make the iteration residuals decreasing while for 'error'
             the solution error estimates should decrease."""),
         ('ls_on', 'float', 0.99999, False,
-         """Start the backtracking line-search by reducing the step, if
+         r"""Start the backtracking line-search by reducing the step, if
             :math:`||d(x^i)|| / ||d(x^{i-1})||` is larger than `ls_on`,
             where :math:`d` is either :math:`f` or :math:`\Delta x`
             depending on `ls_mode`."""),