pyrbit.mle_utils.confidence_ellipse

confidence_ellipse(inferred_parameters, estimated_covariance_matrix, confidence_levels=[0.68, 0.95], ax=None, colors=['#B0E0E6', '#87CEEB'], plot_kwargs={'color': 'red', 'label': 'ML estimate', 'marker': 'D'})

estimated_covariance_matrix = numpy.linalg.inv(J)