F =
Hi Sunday,
The code runs without error (after defining h) here on Answers using 2024b. You're seeing a different result in 2021a?
syms a b c d e f g h % have to define h
F = [0 0 0 a; 0 0 0 b; 0 0 0 c; e f g h]
syms q r s t u
V = [q 0 0 0; r s 0 0; 0 0 t 0; 0 0 0 u]
Vi = inv(V) % calculate inverse of V
K = F*Vi % compute the next generation matrix
eigenvalues= eig(K) % Find the eigenvalues of \(K\)
abs_eigenvalues = abs(eigenvalues) % Compute the absolute values of the eigenvalues.
spectral_radius = max(abs_eigenvalues) %Find the maximum absolute eigenvalue. Error It's show fatal error
