#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