Maybe I got your question wrong and missing something complex in the equations. But my comment would be using two for loops. Outer one for index i=1:n and inner one for j=1:n. This might not be an elegant way of doing it or maybe have some mistales but do comment if you have doubt. I can try to improve.
for i=1:n
%compute S(i), assuming you have some formula for that
for j=1:n
%compute your S(i,j) if you have some formula for that
end
%compute you c(i)
%compute your d(i)
%I assume you can compute vj term outside these two loops code in a separate j=1:n for loop.
end
