def METHOD_NAME(wb):
    uncertainties = wb.model.get_params_dict("free", "uncertainty")
    for param in uncertainties:
        # uncertainty based on the diagonals of the covariance matrix
        matrix_uncertainty = np.sqrt(
            wb.parameter_covariance_matrix.get_label_matrix([param]).matrix[0, 0]
        )
        assert np.isclose(uncertainties[param], matrix_uncertainty)