The estimation algorithm is simulated field algorithm. theta 0.8653265542218813 inverse_phi 0.07843172509982103 psi 0.36919426955263107 running_time in unit of a second 21.029