i have an error:Assignment has more non-singleton rhs dimensions than non-singleton subscripts

Hi

i have this error in line 77 and 81 in matlab:

Assignment has more non-singleton rhs dimensions than non-singleton subscripts

please help me.

h = 1;
Lr = 1.5;
Lz = 1;
dr = 0.01429;
dz = 0.01429;
r = 0:dr:Lr;
z = 0:dz:Lz;
i = 0;
j = 0;
k_m = 1.80;
k_c = 1.65;
G_m = 0.76;
G_c = 1.24;
ro_m = 8.3;
ro_c = 3.7;
for z = 0:dz:Lz
    i=i+1;
    j=0;
    for r = 0:dr:Lr
        j=j+1;
        n=1;
        v_m(i,j) = ((z/h)^n);
        v_c(i,j) = (1-v_m(i,j));
        k(i,j) = k_m+(v_c(i,j)*(k_c-k_m)/(1+(1-v_c(i,j))*(3*(k_c-k_m)/3*k_m+4*G_m)));
        f1 = G_m*(9*k_m+8*G_m)/6*(k_m+2*G_m);
        G(i,j) = G_m+(v_c(i,j)*(G_c-G_m)/(1+(1-v_c(i,j))*((G_c-G_m)/G_m+f1)));
        E(i,j) = 9*k(i,j)*G(i,j)/(3*k(i,j)+inv(G(i,j)));
        nu(i,j) = (3*k(i,j)-2*G(i,j))/2*inv(3*k(i,j)+G(i,j));
        ro(i,j) = v_m(i,j)*ro_m+v_c(i,j)*ro_c;
        landa(i,j) = E(i,j)*nu(i,j)/((1+nu(i,j))*(1-2*nu(i,j)));
        mu(i,j) = E(i,j)/(2*(1+nu(i,j)));
        cl(i,j) = ((landa(i,j)+2*mu(i,j))/ro(i,j))^0.5;
        cs(i,j) = (mu(i,j)/ro(i,j))^0.5;
        q(i,j) = cs(i,j)/cl(i,j);
    end
end
z=0:dz:Lz;
r=0:dr:Lr;
dt=1e-8;
t=0:dt:50e-6;
u=zeros(length(z),length(r),length(t));
w=zeros(length(z),length(r),length(t));
for h=1:length(r)
 u(1,h,:)=0.001;%*sin(pi*t/5e-5);
 w(1,h,:)=0.001;
end
for k=2:length(t)-1
    for i=2:length(z)-1
        for j=2:length(r)-1

              u(i,j,2) = (dt)^2*((1/dr^2)*(u(i+1,j,1)-2*u(i,j,1)+u(i-1,j,1))+(1-q(i,j)^2)*(1/4*dr*dz)*(w(i+1,j+1,1)-w(i+1,j-1,1)...
                  -w(i-1,j+1,1)+w(i-1,j-1,1))+(1/i*dr)*(1/2*dr)*(u(i+1,j,1)-u(i-1,j,1))-(1/(i*dr)^2)*u(i,j,1)...
                  +(q(i,j)^2)*(1/dz^2)*(u(i,j+1,1)-2*u(i,j,1)+u(i,j-1,1)))+2*u(i,j,1)-u(i,j,1);

              w(i,j,2) = (dt)^2*((1/dr^2)*(w(i,j+1,1)-2*w(i,j,1)+w(i,j-1,1))+(q(i,j)^2)*(1/dr^2)*(w(i+1,j,1)-2*w(i,j,1)+w(i-1,j,1)...
                  +(1/i*dr)*(1/2*dr)*(w(i+1,j,1)-w(i-1,j,1)))+(1-q(i,j)^2)*(u(i+1,j+1,1)-u(i+1,j-1,1)-u(i-1,j+1,1)+u(i-1,j-1,1))...
                  +(1/2*dz)*(w(i,j+1,1)-w(i,j-1,1)))+2*w(i,j,1)-w(i,j,1);

          end
      end
  end

for k=2:length(t)-1
    for i=2:length(z)-1
        for j=2:length(r)-1

              u(i,j,k) = (dt)^2*((1/dr^2)*(u(i+1,j,k)-2*u(i,j,k)+u(i-1,j,k))+(1-q^2)*(1/4*dr*dz)*(w(i+1,j+1,k)-w(i+1,j-1,k)...
                  -w(i-1,j+1,k)+w(i-1,j-1,k))+(1/i*dr)*(1/2*dr)*(u(i+1,j,k)-u(i-1,j,k))-(1/(i*dr)^2)*u(i,j,k)...
                  +(q^2)*(1/dz^2)*(u(i,j+1,k)-2*u(i,j,k)+u(i,j-1,k)))+2*u(i,j,k)-u(i,j,k-1);

              w(i,j,k) = (dt)^2*((1/dr^2)*(w(i,j+1,k)-2*w(i,j,k)+w(i,j-1,k))+(q^2)*(1/dr^2)*(w(i+1,j,k)-2*w(i,j,k)+w(i-1,j,k)...
                  +(1/i*dr)*(1/2*dr)*(w(i+1,j,k)-w(i-1,j,k)))+(1-q^2)*(u(i+1,j+1,k)-u(i+1,j-1,k)-u(i-1,j+1,k)+u(i-1,j-1,k))...
                  +(1/2*dz)*(w(i,j+1,k)-w(i,j-1,k)))+2*w(i,j,k)-w(i,j,k-1);

          end
      end
  end