diff --git a/ext/TrixiNLsolveExt.jl b/ext/TrixiNLsolveExt.jl
index 5177c72dcba..3db0e057f12 100644
--- a/ext/TrixiNLsolveExt.jl
+++ b/ext/TrixiNLsolveExt.jl
@@ -73,15 +73,15 @@ end
 # non-linear equations that arise from the relation of the stability polynomial to the Butcher tableau.
 # For details, see Proposition 3.2, Equation (3.3) from 
 # Hairer, Wanner: Solving Ordinary Differential Equations 2
-function Trixi.solve_a_butcher_coeffs_unknown!(a_unknown, num_stages, monomial_coeffs,
+function Trixi.solve_a_butcher_coeffs_unknown!(a_unknown, num_stages, num_stage_evals, monomial_coeffs,
                                                c_s2, c;
-                                               verbose, max_iter = 100000) #TODO: Add an argument of num_stage_evals here
+                                               verbose, max_iter = 100000)
 
     # Define the objective_function
     function objective_function!(c_eq, x)
         return PairedExplicitRK3_butcher_tableau_objective_function!(c_eq, x,
                                                                      num_stages,
-                                                                     num_stages,
+                                                                     num_stage_evals,
                                                                      monomial_coeffs,
                                                                      c_s2)
     end
diff --git a/src/time_integration/paired_explicit_runge_kutta/methods_PERK3.jl b/src/time_integration/paired_explicit_runge_kutta/methods_PERK3.jl
index 17a841c7ac3..0e7dfcc69b9 100644
--- a/src/time_integration/paired_explicit_runge_kutta/methods_PERK3.jl
+++ b/src/time_integration/paired_explicit_runge_kutta/methods_PERK3.jl
@@ -56,7 +56,7 @@ function compute_PairedExplicitRK3_butcher_tableau(num_stages, tspan,
         # Solve the nonlinear system of equations from monomial coefficient and
         # Butcher array abscissae c to find Butcher matrix A
         # This function is extended in TrixiNLsolveExt.jl
-        a_unknown = solve_a_butcher_coeffs_unknown!(a_unknown, num_stages,
+        a_unknown = solve_a_butcher_coeffs_unknown!(a_unknown, num_stages, num_stages,
                                                     monomial_coeffs, cS2, c;
                                                     verbose)
     end