fix: dependencies errors with scipy

src/optymus/methods/_second_order.py

@@ -100,7 +100,7 @@ def penalized_obj(x):
         if jnp.linalg.norm(g) < tol:
             break
 
-        d = jax.scipy.linalg.solve(M, -g)
+        d = jnp.linalg.solve(M, -g)
 
         r = line_search(f=penalized_obj, x=x, d=d, learning_rate=learning_rate)
         x = r['xopt']