In 1893 Giuseppe Lauricella defined and studied four hypergeometric series FA, FB, FC, FD of three variables. They are (Lauricella 1893): 
  
for |x1| + |x2| + |x3| < 1 and 
  
for |x1| < 1, |x2| < 1, |x3| < 1 and 
  
for |x1|1/2 + |x2|1/2 + |x3|1/2 < 1 and 
  
for |x1| < 1, |x2| < 1, |x3| < 1. Here the Pochhammer symbol (q)i indicates the i-th rising factorial of q, i.e. 
  
where the second equality is true for all complex  except
 except  .
. 
These functions can be extended to other values of the variables x1, x2, x3 by means of analytic continuation. 
Lauricella also indicated the existence of ten other hypergeometric functions of three variables. These were named FE, FF, ..., FT and studied by Shanti Saran in 1954 (Saran 1954). There are therefore a total of 14 Lauricella–Saran hypergeometric functions. 
  Generalization to n variables
 These functions can be straightforwardly extended to n variables. One writes for example 
  
where |x1| + ... + |xn| < 1. These generalized series too are sometimes referred to as Lauricella functions. 
When n = 2, the Lauricella functions correspond to the Appell hypergeometric series of two variables: 
  
When n = 1, all four functions reduce to the Gauss hypergeometric function: 
  
Integral representation of FD
 In analogy with Appell's function F1, Lauricella's FD can be written as a one-dimensional Euler-type integral for any number n of variables: 
  
This representation can be easily verified by means of Taylor expansion of the integrand, followed by termwise integration. The representation implies that the incomplete elliptic integral Π is a special case of Lauricella's function FD with three variables: 
  
Finite-sum solutions of FD
 Case 1 :  ,
,  a positive integer
 a positive integer 
One can relate FD to the Carlson R function  via
 via 
 
 
with the iterative sum 
 and
 and  
 
where it can be exploited that the Carlson R function with  has an exact representation (see [1] for more information).
 has an exact representation (see [1] for more information). 
The vectors are defined as 
![{\displaystyle {\overline {b^{*}}}=[{\overline {b}},c-\sum _{i}b_{i}]}](./_assets_/80751fea2d6b7fc97fa4d37bc1ceffc8aa0d8cbe.svg) 
 
![{\displaystyle {\overline {z^{*}}}=[{\frac {1}{1-z_{1}}},\ldots ,{\frac {1}{1-z_{N-1}}},1]}](./_assets_/70acaa1030bfec1c2f47d68025a9ae24dae8cf8b.svg) 
 
where the length of  and
 and  is
 is  , while the vectors
, while the vectors  and
 and  have length
 have length  .
. 
Case 2:  ,
,  a positive integer
 a positive integer 
In this case there is also a known analytic form, but it is rather complicated to write down and involves several steps. See [2] for more information. 
 References
  - Appell, Paul; Kampé de Fériet, Joseph (1926). Fonctions hypergéométriques et hypersphériques; Polynômes d'Hermite (in French). Paris: Gauthier–Villars. JFM 52.0361.13. (see p. 114)
- Exton, Harold (1976). Multiple hypergeometric functions and applications. Mathematics and its applications. Chichester, UK: Halsted Press, Ellis Horwood Ltd. ISBN 0-470-15190-0. MR 0422713.
- Lauricella, Giuseppe (1893). "Sulle funzioni ipergeometriche a più variabili". Rendiconti del Circolo Matematico di Palermo (in Italian). 7 (S1): 111–158. doi:10.1007/BF03012437. JFM 25.0756.01. S2CID 122316343.
- Saran, Shanti (1954). "Hypergeometric Functions of Three Variables". Ganita. 5 (1): 77–91. ISSN 0046-5402. MR 0087777. Zbl 0058.29602. (corrigendum 1956 in Ganita 7, p. 65)
- Slater, Lucy Joan (1966). Generalized hypergeometric functions. Cambridge, UK: Cambridge University Press. ISBN 0-521-06483-X. MR 0201688. (there is a 2008 paperback with ISBN 978-0-521-09061-2)
- Srivastava, Hari M.; Karlsson, Per W. (1985). Multiple Gaussian hypergeometric series. Mathematics and its applications. Chichester, UK: Halsted Press, Ellis Horwood Ltd. ISBN 0-470-20100-2. MR 0834385. (there is another edition with ISBN 0-85312-602-X)