#5671 psi(complex float)

[Barton Willis]

The polygamma functions with a complex float argument return a nounform:

(%i3) psi[0](2.0 + 0.1*%i);

(%o3) psi[0](0.1*%i+2.0)

But Maxima does have code for this case:

(%i3)  bfpsi(0,2.0,15);
(%o3)                         4.22784335098467b-1
(%i4)  bfpsi(0,2.0+%i,15);
(%o4)            5.76674047468581b-1 %i + 5.94650320622477b-1

It seems it returns a big float even when the input is a binary64.