clear
%
% written by HP Huang
%
x1 = zeros(51); y1 = zeros(51); x2 = zeros(2); y2 = zeros(2);
x3 = zeros(2); y3 = zeros(2); x4 = zeros(2); y4 = zeros(2);
x5 = zeros(2); y5 = zeros(2);
%
phi1 = pi; phi2 = atan(2/sqrt(3));
dphi = (phi1-phi2)/(length(x1)-1);
for i = 1:length(x1)
phi = phi1-(i-1)*dphi;
x1(i) = 8*cos(phi);
y1(i) = 4*sin(phi);
phi*180/pi;
end
%
x2(1) = x1(length(x1)); y2(1) = y1(length(x1));
x2(2) = 8; y2(2) = 8*(tan(pi/6));
%
x3(1) = x2(2); y3(1) = y2(2);
x3(2) = x3(1)+1*cos(pi/6); y3(2) = y3(1)-1*sin(pi/6);
%
x4(1) = x3(2); y4(1) = y3(2);
x4(2) = x4(1); y4(2) = 0;
%
x5(1) = x4(2); y5(1) = y4(2);
x5(2) = x1(1); y5(2) = y1(1);
%
x1_rot=x1*cos(-pi/6)-y1*sin(-pi/6);
y1_rot=y1*cos(-pi/6)+x1*sin(-pi/6);
x2_rot=x2*cos(-pi/6)-y2*sin(-pi/6);
y2_rot=y2*cos(-pi/6)+x2*sin(-pi/6);
x3_rot=x3*cos(-pi/6)-y3*sin(-pi/6);
y3_rot=y3*cos(-pi/6)+x3*sin(-pi/6);
x4_rot=x4*cos(-pi/6)-y4*sin(-pi/6);
y4_rot=y4*cos(-pi/6)+x4*sin(-pi/6);
x5_rot=x5*cos(-pi/6)-y5*sin(-pi/6);
y5_rot=y5*cos(-pi/6)+x5*sin(-pi/6);
z0 = 0;
fileID = fopen('half_fish1.txt','w');
%
i1 = 1;
for i2 = 1:length(x1)
fprintf(fileID,'%3i %3i %9.6f %9.6f
%9.6f\n',i1,i2,x1_rot(i2),y1_rot(i2),z0);
end
%
i1 = 2;
for i2 = 1:length(x2)
fprintf(fileID,'%3i %3i %9.6f %9.6f
%9.6f\n',i1,i2,x2_rot(i2),y2_rot(i2),z0);
end
%
i1 = 3;
for i2 = 1:length(x3)
fprintf(fileID,'%3i %3i %9.6f %9.6f
%9.6f\n',i1,i2,x3_rot(i2),y3_rot(i2),z0);
end
%
i1 = 4;
for i2 = 1:length(x4)
fprintf(fileID,'%3i %3i %9.6f %9.6f
%9.6f\n',i1,i2,x4_rot(i2),y4_rot(i2),z0);
end
%
i1 = 5;
for i2 = 1:length(x5)
fprintf(fileID,'%3i %3i %9.6f %9.6f
%9.6f\n',i1,i2,x5_rot(i2),y5_rot(i2),z0);
end
%
fclose(fileID);
%
plot(x1_rot,y1_rot,'k-',x2_rot,y2_rot,'r-',x3_rot,y3_rot,'g',x4_rot,y4_rot,'b-',x5_rot,y5_rot,'m-',...
'LineWidth',2)
axis([-12 12 -12 12])