x=0:0.1
i/2;y=0:0.1
i/2;
[x,y]=meshgrid(x,y);
p1=0.74;
p2=1.2;
p3=2;
p4=5;
p5=1.2;
p6=1.2;
p7=1.2;
l=700;
a=100;
x1k=70;
x2k=140+0.5*l*sin(x);
x3k=202+l*sin(x);
x4k=265+l*sin(x);
x5k=90+0.5*l*sin(x);
x17=-(p7+4*(p1*x1k+p2*x2k+p3*x3k+p4*x4k+p5*x5k)).*sin(y)./(2*a*(sqrt(3)*cos(y)-sin(y)));
x16=zeros(size(x,1),size(x,1));
for i=1:size(x,1)
for j=1:size(x,1)
x16(i,j)=-(p7+4*(p1*x1k+p2*x2k(i,j)+p3*x3k(i,j)+p4*x4k(i,j)+p5*x5k(i,j)))*sin(y(i,j))/(2*a*(sqrt(3)*cos(y(i,j))-sin(y(i,j))));
end
end
figure(1)
surf(x,y,x16);
figure(2)
surf(x,y,x17);