MODAL Receptance matrix calculation

I am triying create receptance matrix of beam exited by harmonic force. I wrote the following code;

for i=1:10

for k=1:10

for r=1:10

summe(i,k)=0;

w=0:0.01:2000;

summe(i,k)=summe(i,k)+(phi_c(i,r).*phi_c(k,r))./(wn_c(r).^2-w.^2+0.05i*wn_c(r).^2)

end

end

end

But there exist an error like "Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is

1-by-200001.

Error in (line 87)

summe(i,k)=summe(i,k)+(phi_c(i,r).*phi_c(k,r))./(wn_c(r).^2-w.^2+0.05i*wn_c(r).^2)"

What is the wrong?