机器人腿部落点 maltab问题，求大神解决

二楼是我写的程序
for a=1:10
H(a)=rand(1)*0.1+0.15;
d4(a)=400*H(a);
l4(a)=2*d4(a);
l2(a)=400-3*d4(a);
for i=1:1000
theta1(i)=60*rand(1)-30;
theta2(i)=90*rand(1)-60;
theta4(i)=60*rand(1)+20;
theta3(i)=90*rand(1)+40;
x(i)=l4(a)*cosd(theta1(i))*sind(theta2(i)+theta3(i))*cosd(theta4(i))-l4(a)*sind(theta1(i))*sind(theta4(i))+cosd(theta1(i))*cosd(theta2(i)+theta3(i))*d4(a)+l2(a)*cosd(theta1(i)+cosd(theta2(i)));
y(i)=l4(a)*sind(theta1(i))*sind(theta2(i)+theta3(i))*cosd(theta4(i))+l4*cosd(theta1(i))*sind(theta4(i))+d4(a)*sind(theta1(i))*cosd(theta2(i)+theta3(i))+l2(a)*sind(theta1(i)+cosd(theta2(i)));
z(i)=-l4(a)*cosd(theta2(i)+theta3(i))*cosd(theta4(i))+sind(theta2(i)+theta3(i))*d4(a)+l2(a)*sind(theta2(i));
end
ymax(a)=max(y);
ymin(a)=min(y);
n(a)= (ymax(a)-ymin(a))/200;
s(a)=0;
for j=1:1:200
xmin(j)=200;
xmax(j)=0;
for i=1:1000
if(y(i)>ymax(a)-n(a)*j)&(y(i)<=ymax(a)-n(a)*(j-1))&(x(i)>xmax(j))
xmax(j)=x(i);
end
if (y(i)>ymax(a)-n(a)*j)&(y(i)<=ymax(a)-n(a)*(j-1))&(x(i)<xmin(j))
xmin(j)=x(i);
end
end
s(a)=(xmax(j)-xmin(j))+s(a);
end
sum(a)=s(a)*n(a);
end
plot(H, sum,'rx')

没人吗，我自己顶一下啊，着急啊

顶下啊