seg1_x = 100:0.1:150;
seg1_y = 100:0.3:250;
seg1_z = getdk(seg1_y);

seg2_x = ones(1,501) * 150;
seg2_y = 250:0.1:300;
seg2_z = getdk(seg2_y);

seg3_x = 150:0.3:300;
seg3_y = 300:0.1:350;
seg3_z = getdk(seg3_y);

seg4_x = 400:0.1:450;
seg4_y = 400:0.1:450;
seg4_z = getdk(seg4_y);

seg5_x = 450:-0.1:400;
seg5_y = 450:0.1:500;
seg5_z = getdk(seg5_y);

figure;hold on;
plot3(seg1_x, seg1_y, seg1_z, 'b');

plot3(seg2_x, seg2_y, seg2_z);
plot3(seg3_x, seg3_y, seg3_z);
plot3(seg4_x, seg4_y, seg4_z);
plot3(seg5_x, seg5_y, seg5_z);
