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]) 
