function [equations] = system_of_equation2array(system_of_equation)
equations_divided_in_cells = children(system_of_equation);
number_of_equations = size(equations_divided_in_cells,1);
equations = sym(zeros(number_of_equations,1));
for i = 1 : number_of_equations
lhs_of_equation_i = equations_divided_in_cells{i}{1};
rhs_of_equation_i = equations_divided_in_cells{i}{2};
equations(i) = lhs_of_equation_i == rhs_of_equation_i;
end
end
Usage:
syms x(t) y(t) t a
dynamical_system = [
diff(x,t) == a*x
diff(y,t) == -y
];
% BEFORE:
>> dynamical_system(1)
ans =
subs(diff(x(t), t), t, 1) == a*x(1)
subs(diff(y(t), t), t, 1) == -y(1)
% AFTER:
dynamical_system = system_of_equation2array(dynamical_system);
>> dynamical_system(1)
ans =
diff(x(t), t) == a*x(t)