diff --git a/cvc5_pythonic_api/cvc5_pythonic.py b/cvc5_pythonic_api/cvc5_pythonic.py
index bae370f..d5e1e43 100644
--- a/cvc5_pythonic_api/cvc5_pythonic.py
+++ b/cvc5_pythonic_api/cvc5_pythonic.py
@@ -6325,15 +6325,14 @@ def statistics(self):
         return self.solver.getStatistics()
 
     def unsat_core(self):
-        """Return a subset (as a list of Bool expressions) of the assumptions provided to the last check().
+        """Return a subset (as a list of Bool expressions) of the
+        assertions and assumptions provided to the last check().
 
-        These are the unsat ("failed") assumptions.
-
-        To enable this, set the option "produce-unsat-assumptions" to true.
+        To enable this, set the option "produce-unsat-cores" to true.
 
         >>> a,b,c = Bools('a b c')
         >>> s = Solver()
-        >>> s.set('produce-unsat-assumptions','true')
+        >>> s.set('produce-unsat-cores','true')
         >>> s.add(Or(a,b),Or(b,c),Not(c))
         >>> s.check(a,b,c)
         unsat
@@ -6347,8 +6346,8 @@ def unsat_core(self):
         >>> s.check(a,b)
         sat
         """
-        core = self.solver.getUnsatAssumptions()
-        return [BoolRef(c) for c in core]
+        core = self.solver.getUnsatCore()
+        return [_to_expr_ref(c, Context(self)) for c in core]
 
 
 def SolverFor(logic, ctx=None, logFile=None):
diff --git a/test/pgm_outputs/unsat_assumptions.py.out b/test/pgm_outputs/unsat_cores.py.out
similarity index 100%
rename from test/pgm_outputs/unsat_assumptions.py.out
rename to test/pgm_outputs/unsat_cores.py.out
diff --git a/test/pgms/unsat_assumptions.py b/test/pgms/unsat_cores.py
similarity index 57%
rename from test/pgms/unsat_assumptions.py
rename to test/pgms/unsat_cores.py
index 3ff4e99..2cf9288 100644
--- a/test/pgms/unsat_assumptions.py
+++ b/test/pgms/unsat_cores.py
@@ -3,24 +3,21 @@
 
 def reset_solver(s):
     s.reset()
-    s.set('produce-unsat-assumptions','true')
+    s.set('produce-unsat-cores','true')
 
-def validate_unsat_assumptions(assumptions, core):
-    # checks that the produced unsat assumptions (core) match the assumptions (assumptions) sent to the check function
-    return sum([c in assumptions for c in core]) == len(core) 
+def validate_unsat_core(input_formulas, core):
+    # checks that the produced unsat core match the input formulas sent to the check function
+    return sum([c in input_formulas for c in core]) == len(core)
 
 
-def check_unsat_assumptions(assertions, core):
-    # This function checks wether, given assertions,  the produced unsat assumptions (core) also lead to unsat result
+def check_unsat_core(core):
+    # This function checks whether the unsat core is unsatisfiable
     slvr = Solver()
-    slvr.set('produce-unsat-assumptions','true')
-    for a in assertions:
-        slvr.add(a)
     return s.check(*core) == unsat
 
 
 # To make make sure the unsat_core function works there should be at least one nontrivial solution - a solution that doesn't contain all the assumptions sent in the check function.
-nontrivial_counter = 0  
+nontrivial_counter = 0
 
 p1, p2, p3 = Bools('p1 p2 p3')
 x, y = Ints('x y')
@@ -37,10 +34,9 @@ def check_unsat_assumptions(assertions, core):
 
 core = s.unsat_core()
 
-
-assert validate_unsat_assumptions(assumptions,core)
-assert check_unsat_assumptions(assertions,core)
-if len(core) < len(assumptions):
+assert validate_unsat_core(assertions + assumptions,core)
+assert check_unsat_core(core)
+if len(core) < len(assumptions) + len(assertions):
     nontrivial_counter += 1
 
 # example 2 - booleans
@@ -60,11 +56,11 @@ def check_unsat_assumptions(assertions, core):
 assumptions = [a,b,c]
 result = s.check(*assumptions)
 
-unsat_core = s.unsat_core()
+core = s.unsat_core()
 
-assert validate_unsat_assumptions(assumptions,unsat_core) 
-assert check_unsat_assumptions(assertions,assumptions)
-if len(unsat_core) < len(assumptions):
+assert validate_unsat_core(assertions + assumptions, core)
+assert check_unsat_core(core)
+if len(core) < len(assumptions) + len(assertions):
     nontrivial_counter += 1
 
 # example 3 - booleans
@@ -83,11 +79,11 @@ def check_unsat_assumptions(assertions, core):
 assumptions = [a,b,c,d]
 result = s.check(*assumptions)
 
-unsat_core = s.unsat_core()
+core = s.unsat_core()
 
-assert validate_unsat_assumptions(assumptions,unsat_core)
-assert check_unsat_assumptions(assertions,assumptions)
-if len(unsat_core) < len(assumptions):
+assert validate_unsat_core(assumptions + assertions, core)
+assert check_unsat_core(core)
+if len(core) < len(assumptions) + len(assertions):
     nontrivial_counter += 1
 
 # example 4 - reals
@@ -108,13 +104,13 @@ def check_unsat_assumptions(assertions, core):
 assumptions = [x > 0, y > 0, z > 0]
 result = s.check(*assumptions)
 
-unsat_core = s.unsat_core()
+core = s.unsat_core()
 
-assert validate_unsat_assumptions(assumptions,unsat_core)
-assert check_unsat_assumptions(assertions,assumptions)
-if len(unsat_core) < len(assumptions):
+assert validate_unsat_core(assumptions + assertions, core)
+assert check_unsat_core(core)
+if len(core) < len(assumptions) + len(assertions):
     nontrivial_counter += 1
-    
+
 
 # example 5 - strings
 
@@ -135,16 +131,14 @@ def check_unsat_assumptions(assertions, core):
 
 result = s.check( Length(s2) < 2)
 
-unsat_core = s.unsat_core()
+core = s.unsat_core()
 
-assert validate_unsat_assumptions([Length(s2) < 2], unsat_core)
-assert check_unsat_assumptions(assertions,[ Length(s2) < 2 ])
-if len(unsat_core) < len([ Length(s2) < 2 ]):
+assert validate_unsat_core([Length(s2) < 2] + assertions, core)
+assert check_unsat_core(core)
+if len(core) < len([ Length(s2) < 2 ]) + len(assertions):
     nontrivial_counter += 1
 
 # check that there is at least one nontrivial unsat core
 assert nontrivial_counter >= 1
 
 print('success')
-
-