En matemática, la fórmula de Riemann–Siegel es una fórmula asintótica para el error que se comete en la ecuación funcional aproximada de la función zeta de Riemann, una aproximación de la función zeta mediante las suma de dos series de Dirichlet finitas. Esta fue encontrada por Siegel (1932) en unos manuscritos no publicados de Bernhard Riemann  alrededor de los años 1850s. Siegel convirtió esta en la fórmula integral de Riemann–Siegel, una expresión para la función zeta en la que intervienen integrales de contorno. Es a menudo usada para calcular valores de la fórmula de Riemann–Siegel, a veces en combinación con el algoritmo de Odlyzko–Schönhage, lo cual aumenta la velocidad de los cálculos considerablemente.
Si M y N son dos números enteros no negativos, entonces la función zeta es igual a 

  
    
      
        ζ
        (
        s
        )
        =
        
          ∑
          
            n
            =
            1
          
          
            N
          
        
        
          
            1
            
              n
              
                s
              
            
          
        
        +
        γ
        (
        1
        −
        s
        )
        
          ∑
          
            n
            =
            1
          
          
            M
          
        
        
          
            1
            
              n
              
                1
                −
                s
              
            
          
        
        +
        R
        (
        s
        )
      
    
    {\displaystyle \zeta (s)=\sum _{n=1}^{N}{\frac {1}{n^{s}}}+\gamma (1-s)\sum _{n=1}^{M}{\frac {1}{n^{1-s}}}+R(s)}
  

donde:

  
    
      
        
          γ
          (
          s
          )
          =
          
            π
            
              1
              
                /
              
              2
              −
              s
            
          
          Γ
          (
          s
          
            /
          
          2
          )
          
            /
          
          Γ
          (
          (
          1
          −
          s
          )
          
            /
          
          2
          )
        
      
    
    {\displaystyle \displaystyle \gamma (s)=\pi ^{1/2-s}\Gamma (s/2)/\Gamma ((1-s)/2)}
      es un factor que aparece en la ecuación funcional ζ(s) = γ(s) ζ(1 − s), y
  
    
      
        R
        (
        s
        )
        =
        
          
            
              −
              Γ
              (
              1
              −
              s
              )
            
            
              2
              π
              i
            
          
        
        ∫
        
          
            
              (
              −
              x
              
                )
                
                  s
                  −
                  1
                
              
              
                e
                
                  −
                  N
                  x
                
              
              d
              x
            
            
              
                e
                
                  x
                
              
              −
              1
            
          
        
      
    
    {\displaystyle R(s)={\frac {-\Gamma (1-s)}{2\pi i}}\int {\frac {(-x)^{s-1}e^{-Nx}dx}{e^{x}-1}}}
       es una integral de contorno, cuyo contorno comienza y termina en +∞ y circunferencias de singularidades de un valor absoluto a lo sumo 2πM. La ecuación funcional aproximada da una estimación del tamaño del término error.Siegel (1932) y Edwards (1974) derivaron la fórmula de Riemann–Siegel formula para este fin, mediante la aplicación del método del descenso más rápido a esta integral, lo cual da una expansión asintótica para el término error R(s) como una serie de potencias negativas de Im(s). 
En aplicaciones,  s está usualmente en la línea crítica, y los enteros positivos M y N son escogidos de tal manera que estén cercanos a (2π Im(s))1/2. Gabcke (1979) encontró unos buenos límites para el término de error de la fórmula Riemann–Siegel.