fixed bug where old IVP Solver Toolbox syntax was used for matrix-vec…

…tor IVP conversions
tamaskis · Dec 8, 2022 · 73cdce0 · 73cdce0
1 parent 015bf38
commit 73cdce0
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/EXAMPLES.mlx b/EXAMPLES.mlx
diff --git a/Riccati_Differential_Equation.pdf b/Riccati_Differential_Equation.pdf
diff --git a/solve_riccati_ode.m b/solve_riccati_ode.m
@@ -7,7 +7,7 @@
 %   [t,P] = solve_riccati_ode(A,B,Q,R,S,PT,tspan)
 %
 % Copyright © 2021 Tamas Kis
-% Last Update: 2022-08-28
+% Last Update: 2022-12-08
 % Website: https://tamaskis.github.io
 % Contact: tamas.a.kis@outlook.com
 %
@@ -92,7 +92,7 @@
     dPdt = @(t,P) -(A.'*P+P*A-(P*B+S)/R*(B.'*P+S.')+Q);
 
     % converts matrix-valued ODE to vector-valued ODE
-    dydt = mat2vec_ode(dPdt);
+    dydt = mat2vec_ODE(dPdt,n);
 
     % converts matrix initial condition to vector initial condition
     yT = mat2vec_IC(PT);
@@ -102,7 +102,7 @@
 
     % transforms solution matrix for vector-valued ODE into solution array
     % for matrix-valued ODE
-    P = vec2mat_sol(y);
+    P = vec2mat_sol(y,n);
 
     % reorders t so that time is increasing
     t = flipud(t);