Fixed two typos.

Author: JonasBreuling

README.md

@@ -104,7 +104,7 @@ More examples are given in the [examples](examples/) directory, which includes
 
 ## Work-precision
 
-In order to investigate the work precision of the implemented solvers, we use different DAE examples with differentiation index 1, 2 and 3 as well as IDE example.
+In order to investigate the work precision of the implemented solvers, we use different DAE examples with differentiation index 1, 2 and 3 as well as IDE examples.
 
 ### Index 1 DAE - Brenan
 

scipy_dae/integrate/_dae/radau.py

@@ -682,8 +682,8 @@ def _step_impl(self):
                 else:
                     error_embedded = -self.solve_lu(LU_real, F)
             else:
-                # compute implicit embedded method with `newton_iter_embedded`` iterations
-                yp_hat_new0 = -(y / h + + b_hat1_implicit * yp + self.b_hat_implicit @ Yp) / MU_REAL 
+                # compute implicit embedded method with `newton_iter_embedded` iterations
+                yp_hat_new0 = -(y / h + b_hat1_implicit * yp + self.b_hat_implicit @ Yp) / MU_REAL 
                 y_hat_new = y_new.copy() # initial guess
                 for _ in range(self.newton_iter_embedded):
                     yp_hat_new = yp_hat_new0 + y_hat_new / (h * MU_REAL)