figure
hold on
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(pts{1}(conn(i,idx)),pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end
axis equal

if(sym_planes >= 1)
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(pts{1}(conn(i,idx)),-pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end  
end

if(sym_planes == 2)
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(-pts{1}(conn(i,idx)),pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
  plot3(-pts{1}(conn(i,idx)),-pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end  
end

clear idx


g = 9.81;
d = 5.0;
rho = 1025;
S = pi*d^2/4;
vol = pi*d^3/12;
kd = k*d;

Z = F{7,3}./(-w.^2.*(real(F{3,3})+rho*g)+sqrt(-1)*w.*imag(F{3,3})+rho*g*S);


figure
plot(kd,real(F{3,3})/(rho*vol),kd,imag(F{3,3})./(rho*vol*w),kd,abs(Z));
xlabel('w^2*d/g');
legend('a22/(rho*V)','b22/(rho*V*w)');


%figure
%plot(kd,real(F{7,3})/(rho*vol*g),kd,imag(F{7,3})/(rho*vol*g));