#ifndef PSI_FUNCTIONS_H
 #define PSI_FUNCTIONS_H
 
 namespace psi {
 
 // Digamma function
 double Psi_0(double x);
 
 // Trigamma function
 double Psi_1(double x);
 
 } // namespace psi
 
 #endif // PSI_FUNCTIONS_H