Remove non "-->" arrow versions (except when they are given as example)

Author: TorkelE

docs/src/catalyst_applications/advanced_simulations.md

@@ -375,10 +375,10 @@ Let us consider the Brusselator model:
 using Catalyst, DifferentialEquations, Plots
 
 brusselator = @reaction_network begin
-    A, ∅ → X
-    1, 2X + Y → 3X
-    B, X → Y
-    1, X → ∅
+    A, ∅ --> X
+    1, 2X + Y --> 3X
+    B, X --> Y
+    1, X --> ∅
 end
 u0 = [:X => 1.0, :Y => 0.0]
 tspan = (0.0, 50.0)

docs/src/catalyst_functionality/dsl_description.md

@@ -211,7 +211,7 @@ end
 Pair 3:
 ```@example tut2
 rn1 = @reaction_network begin
-  (1.0,2.0,3.0), (S1,S2,S3) → (P1,P2,P3)
+  (1.0,2.0,3.0), (S1,S2,S3) --> (P1,P2,P3)
 end
 ```
 ```@example tut2
@@ -278,7 +278,7 @@ using SpecialFunctions
 rn = @reaction_network begin
   2.0*X^2, 0 --> X + Y
   t*gamma(Y), X --> ∅
-  pi*X/Y, Y → ∅
+  pi*X/Y, Y --> ∅
 end
 ```
 where here `t` always denotes Catalyst's time variable. Please note that many

docs/src/catalyst_functionality/example_networks/basic_CRN_examples.md

@@ -99,10 +99,10 @@ using Catalyst, DifferentialEquations, Plots
 
 rs = @reaction_network begin
     @parameters A B
-    A, ∅ → X
-    1, 2X + Y → 3X
-    B, X → Y
-    1, X → ∅
+    A, ∅ --> X
+    1, 2X + Y --> 3X
+    B, X --> Y
+    1, X --> ∅
 end
 tspan = (0.0,50.0)
 u0 = [:X => 1.0, :Y => 1.0]

docs/src/faqs.md

@@ -288,7 +288,7 @@ ModelingToolkit. e.g., this is should work
 using Catalyst
 myHill(x) = 2*x^3/(x^3+1.5^3)
 rn = @reaction_network begin
-  myHill(X), ∅ → X
+  myHill(X), ∅ --> X
 end
 ```
 In some cases, it may be necessary or desirable to register functions with

docs/src/inverse_problems/parameter_estimation.md

@@ -21,10 +21,10 @@ Next, we declare our model. For our example, we will use the Brusselator, a
 simple oscillator.
 ```@example pe1
 brusselator = @reaction_network begin
-    A, ∅ → X
-    1, 2X + Y → 3X
-    B, X → Y
-    1, X → ∅
+    A, ∅ --> X
+    1, 2X + Y --> 3X
+    B, X --> Y
+    1, X --> ∅
 end
 p_real = [:A => 1., :B => 2.]
 ```

docs/src/inverse_problems/petab_ode_param_fitting.md

@@ -26,7 +26,7 @@ true_sol = solve(oprob_true, Tsit5())
 data_sol = solve(oprob_true, Tsit5(); saveat=1.0)
 data_ts = data_sol.t[2:end]
 data_vals = (0.8 .+ 0.4*rand(10)) .* data_sol[:P][2:end]
-
+→
 # Plots the true solutions and the (synthetic data) measurements.
 using Plots
 plot(true_sol; idxs=:P, label="True solution", lw=8)