I am trying to implement the following loop inside of matlab

The difficulty I am having is implementing the sum sign such that I get the correct values. The first row of my X matrix is correct but the second row is not and I cant figure out why.

J = [200 100 300];
k = [3e6 2e6 4e6];

n = 3;
points = 300;

for j = 1:points
    w(j)=j;
end

x = zeros(n+1,points);
x(1,:)=1;
for j = 1:points
    sumJx = 0;
    for i = 2:n+1
        sumJx = sum(J(1:i-1)*x(1:i-1,j));
        x(i,j) = x(i-1) - (w(j)^2 / k(i-1)) * sumJx;   
    end
end
f_omega = x(n+1,:);
plot(w,f_omega);
ylim([-40 3])