Yes, solve
takes your second input, a_tilde_sym(iRange)
to be a set of equations rather than the variable(s) to be solved for. There's a good reason for this. Your first two equations (bounday_conditions
(sic)) are in terms of the scalar variables a_tilde2
and a_tilde3
, not the vector [a_tilde2;a_tilde3]
. You can solve this by explicitly indicating your variables:
s = solve(bounday_conditions, a_tilde_sym(2), a_tilde_sym(3))
or
[a_tilde2,a_tilde3] = solve(bounday_conditions, a_tilde_sym(2), a_tilde_sym(3))