I haven't anywhere seen symsum doing the indexing in a structured data. Do you have to use symsum only?
I believe something on the following lines could be helpful, please do verify if the following equation does what you are trying to do.
sum(S(i,1:n) .* (x(1:n) - beta*y(1:n)) - gamma * (x(i) - beta*(y(i))))