yes, i meant that with adaptive = false. Nevertheless, i get an output of
NaN + NaN*im
when evaluating
psi(a0+50.0*b0*im,x0,T0)
even after adding adaptive = false in the psi function definition, whereas the above-mentioned MATLAB’s package cfaffine.m function
cfaffine(-a0*i+b0*50,x0,T0,K0,K1,H0,H1,R0,R1,[],[],[],1)
happily computes the finite number of
0.3534 + 0.2271i.