Fix deprecated ODEProblem API calls in documentation

Updated documentation to use the new ODEProblem API that was introduced
to replace the deprecated 4-argument signature. Changes include:

- Changed ODEProblem(sys, u0, tspan, p) to ODEProblem(sys, [u0; p], tspan)
- Changed ODEProblem(sys, [], tspan, []) to ODEProblem(sys, [], tspan)

This fixes the deprecation warnings in the "Copy-Pastable Simplified Example"
and other tutorial sections.

Files modified:
- docs/src/tutorials/ode_modeling.md: Line 39
- docs/src/tutorials/programmatically_generating.md: Line 52
- docs/src/basics/Events.md: Lines 421, 442
- docs/src/basics/Composition.md: Line 60

Fixes #3996

🤖 Generated with [Claude Code](https://claude.com/claude-code)

Co-Authored-By: Claude <[email protected]>

Author: ChrisRackauckas

docs/src/basics/Composition.md

@@ -57,7 +57,7 @@ p = [decay1.a => 0.1
      decay2.a => 0.2]
 
 using OrdinaryDiffEq
-prob = ODEProblem(simplified_sys, x0, (0.0, 100.0), p)
+prob = ODEProblem(simplified_sys, [x0; p], (0.0, 100.0))
 sol = solve(prob, Tsit5())
 sol[decay2.f]
 ```

docs/src/basics/Events.md

@@ -418,7 +418,7 @@ ev = ModelingToolkit.SymbolicDiscreteCallback(
 @mtkcompile sys = System(
     D(x) ~ c * cos(x), t, [x], [c]; discrete_events = [ev])
 
-prob = ODEProblem(sys, [x => 0.0], (0.0, 2pi), [c => 1.0])
+prob = ODEProblem(sys, [x => 0.0, c => 1.0], (0.0, 2pi))
 sol = solve(prob, Tsit5())
 sol[c]
 ```
@@ -439,7 +439,7 @@ will be saved. If we repeat the above example with `c` not a `discrete_parameter
 @mtkcompile sys = System(
     D(x) ~ c * cos(x), t, [x], [c]; discrete_events = [1.0 => [c ~ Pre(c) + 1]])
 
-prob = ODEProblem(sys, [x => 0.0], (0.0, 2pi), [c => 1.0])
+prob = ODEProblem(sys, [x => 0.0, c => 1.0], (0.0, 2pi))
 sol = solve(prob, Tsit5())
 sol.ps[c] # sol[c] will error, since `c` is not a timeseries value
 ```

docs/src/tutorials/ode_modeling.md

@@ -36,7 +36,7 @@ end
 
 using OrdinaryDiffEq
 @mtkcompile fol = FOL()
-prob = ODEProblem(fol, [], (0.0, 10.0), [])
+prob = ODEProblem(fol, [], (0.0, 10.0))
 sol = solve(prob)
 
 using Plots

docs/src/tutorials/programmatically_generating.md

@@ -49,7 +49,7 @@ eqs = [D(x) ~ (h - x) / τ] # create an array of equations
 # Perform the standard transformations and mark the model complete
 # Note: Complete models cannot be subsystems of other models!
 fol = mtkcompile(model)
-prob = ODEProblem(fol, [], (0.0, 10.0), [])
+prob = ODEProblem(fol, [], (0.0, 10.0))
 using OrdinaryDiffEq
 sol = solve(prob)