matlabで配列の抽出についての質問です。

前提・実現したいこと

matlabで火山ガスが拡散プルーム方式で風向きが変わったときの排出されたときの測定シミュレーションを行っています。ある直線式(このプログラムでいうe3(j)=(b/a).*(x-x1)+y1)を通る部分だけのプルームを求めたいです。

発生している問題・エラーメッセージ

エラーメッセージ

該当のソースコード

alpha=-45;%%噴煙の角度(風向き)
a=-1*tand(alpha);%%傾き
b=-1;
x1=1000;%%ライダーのx座標
y1=0;%%ライダーのy座標

for x=0:10:wx
for y=-1000:1:wy
X=x.*cosd(alpha)-y.*sind(alpha);
Y=x.*sind(alpha)+y.*cosd(alpha);

      e0(j)=x;
      e1(j)=y;
      e2(j)=a.*x;%%噴煙の真ん中を通る線
      e3(j)=(b/a).*(x-x1)+y1;%%噴煙への垂線
      
      sigumay=(cy/sqrt(2)).*((X).^(1-(n/2)));%%y軸方向の拡散幅
      sigumaz=(cz/sqrt(2)).*((X).^(1-(n/2)));%%ｚ軸方向の拡散幅
      
      fy=exp(-(Y).*(Y)./(2.*sigumay.*sigumay));%%y軸方向の濃度分布           
      fz=exp(-(wz-He).^2./(2.*sigumaz.*sigumaz))+exp(-(wz+He).^2./(2.*sigumaz.*sigumaz));%%z軸方向の濃度分布
      
      pulm(j)=Q.*fy.*fz./(2*pi.*sigumay.*sigumaz.*u);%%プルーム式
     
      
      j=j+1;

end
i=i+1;
end

plot3(e0,e1,pulm)
grid on
hold on
xlim([0 1000])
ylim([-1000 1000])
zlim([0 9.0e+24])
xlabel('x軸','FontSize',10)
ylabel('y軸','FontSize',10)
zlabel('濃度分布[個/㎥]','FontSize',10)