octave:726> pwd ans = /home/aaron/Documents/3Ddrawings/BirdBathBowl octave:727> clock ans = Columns 1 through 3: 2019.00000000000000 9.00000000000000 2.00000000000000 Columns 4 through 6: 12.00000000000000 5.00000000000000 4.94767900000000 octave:728> gleen.p=p; octave:729> gleen.m=m; octave:730> disp(horzcat(p,m')); -9.949999999999999 -0.020000000000000 0.500000000000000 -5.000000000000000 2.000000000000000 45.000000000000000 0.000000000000000 10.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -75.000000000000000 6.000000000000000 -10.000000000000000 -86.000000000000000 octave:731> v=sampleSpline(p,m); size(v,1) 1 -9.94999999999999929 -0.02000000000000000 0.00872686779075879 -5.00000000000000000 2.00000000000000000 0.99999999999999989 0.00872686779075879 0.04177549265876909 0.00785900424973544 2 -5.000000000000000 2.000000000000000 1.000000000000000 0.000000000000000 10.000000000000000 0.000000000000000 0.9999999999999999 0.5600000000000002 -0.0880000000000000 3 0.000000000000000 10.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -0.453589838486225 0.010717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 6.000000000000000 -10.000000000000000 -14.300666256711942 -3.73205080756888 -8.23523212815030 1.96728293571918 ans = 86 octave:732> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:733> plot(v(:,1),v(:,2),'.-') octave:734> p p = -9.949999999999999 -0.020000000000000 -5.000000000000000 2.000000000000000 0.000000000000000 10.000000000000000 5.000000000000000 0.000000000000000 6.000000000000000 -10.000000000000000 octave:735> p=[-15,-0.02;-5,5;0,30;5,0;7,-10] p = -15.000000000000000 -0.020000000000000 -5.000000000000000 5.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 7.000000000000000 -10.000000000000000 octave:736> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -5.00000000000000e+00 5.00000000000000e+00 1.00000000000000e+00 8.72686779075879e-03 4.88546264418483e-02 4.72686779075873e-05 2 -5.000000000000000 5.000000000000000 1.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 1.000000000000000 2.600000000000001 -0.360000000000000 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 7.000000000000000 -10.000000000000000 -14.300666256711942 -3.73205080756888 3.38238393592485 -2.00817926607020 ans = 100 octave:737> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:738> plot(v(:,1),v(:,2),'.-') octave:739> m m = Columns 1 through 3: 0.500000000000000 45.000000000000000 0.000000000000000 Columns 4 and 5: -75.000000000000000 -86.000000000000000 octave:740> m(5)=-80; octave:741> p p = -15.000000000000000 -0.020000000000000 -5.000000000000000 5.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 7.000000000000000 -10.000000000000000 octave:742> p(2,:)=[-5 10]; octave:743> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -5.00000000000000e+00 1.00000000000000e+01 1.00000000000000e+00 0.00872686779075879 0.19885462644184831 -0.00995273132209241 2 -5.000000000000000 10.000000000000000 1.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 1.000000000000000 2.000000000000000 -0.280000000000000 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 7.000000000000000 -10.000000000000000 -5.671281819617707 -3.732050807568878 -0.932308282622269 0.149166843203354 ans = 85 octave:744> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:745> plot(v(:,1),v(:,2),'.-') octave:746> p p = -15.000000000000000 -0.020000000000000 -5.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 7.000000000000000 -10.000000000000000 octave:747> p(5,:)=[6.5,-10]; octave:748> p(2,:)=[-8 10]; octave:749> octave:749> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 1.00000000000000e+00 0.00872686779075879 0.46811885410059950 -0.03783939336870171 2 -8.000000000000000 10.000000000000000 1.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 0.9999999999999999 0.6875000000000000 -0.0625000000000000 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 6.500000000000000 -10.000000000000000 -5.671281819617707 -3.73205080756888 -4.57641104349636 1.74666698050967 ans = 93 octave:750> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:751> plot(v(:,1),v(:,2),'.-') octave:752> m m = Columns 1 through 3: 0.500000000000000 45.000000000000000 0.000000000000000 Columns 4 and 5: -75.000000000000000 -80.000000000000000 octave:753> m(5)=-78; octave:754> m(2)=60; octave:755> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 1.73205080756888e+00 0.00872686779075879 0.36354016730504568 -0.02289958096933688 2 -8.000000000000000 10.000000000000000 1.732050807568877 0.000000000000000 30.000000000000000 0.000000000000000 1.7320508075688767 0.5044872981077808 -0.0510617061317363 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 6.500000000000000 -10.000000000000000 -4.704630109478456 -3.73205080756888 -5.22084551692253 2.17628996279378 ans = 79 octave:756> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:757> plot(v(:,1),v(:,2),'.-') octave:758> m(5)=-76; octave:759> m(2)=70; octave:760> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 6.500000000000000 -10.000000000000000 -4.010780933535846 -3.73205080756888 -5.68341163421760 2.48466737432383 ans = 84 octave:761> plot(v(:,1),v(:,2),'.-') octave:762> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:763> plot(v(:,1),v(:,2),'.-') octave:764> m m = Columns 1 through 3: 0.500000000000000 70.000000000000000 0.000000000000000 Columns 4 and 5: -75.000000000000000 -76.000000000000000 octave:765> m(5)=-75; octave:766> p p = -15.000000000000000 -0.020000000000000 -8.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 6.500000000000000 -10.000000000000000 octave:767> p(5,1)=5+(10/cosd(90+m(4)))*sind(90+m(4)) p = -15.000000000000000 -0.020000000000000 -8.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 7.679491924311227 -10.000000000000000 octave:768> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -3.732050807568878 0.000000000000000 -2.853589838486225 0.330717967697245 4 5.000000000000000 0.000000000000000 -3.732050807568878 7.679491924311227 -10.000000000000000 -3.732050807568878 -3.73205080756888e+00 7.42243772059697e-16 -1.84672764595224e-16 ans = 76 octave:769> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:770> plot(v(:,1),v(:,2),'.-') octave:771> p=p(1:4,:); octave:772> m=m(1:4,:); error: A(I,J): row index out of bounds; value 4 out of bound 1 octave:772> m=m(1:4); octave:773> m(4)=76; octave:774> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:775> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 4.010780933535846 0.000000000000000 -4.402156186707169 0.640431237341434 ans = 81 octave:776> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:777> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:778> plot(v(:,1),v(:,2),'.-') octave:779> m(4)=-76; octave:780> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -4.010780933535846 0.000000000000000 -2.797843813292832 0.319568762658566 ans = 73 octave:781> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:782> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:783> plot(v(:,1),v(:,2),'.-') octave:784> m(4)=-78; octave:785> disp(horzcat(p,m')); -15.000000000000000 -0.020000000000000 0.500000000000000 -8.000000000000000 10.000000000000000 70.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -78.000000000000000 octave:786> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -4.704630109478456 0.000000000000000 -2.659073978104309 0.291814795620862 ans = 74 octave:787> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:788> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:789> plot(v(:,1),v(:,2),'.-') octave:790> m(4)=-80; octave:791> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -8.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.21847922274993925 -0.00217658889003596 2 -8.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2506306451363446 -0.0351956653210215 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 72 octave:792> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:793> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:794> plot(v(:,1),v(:,2),'.-') octave:795> p(2,:)=[-6 10]; octave:796> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -6.00000000000000e+00 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.06389653832931785 0.00653750148862610 2 -6.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.747477419454622 0.750840860181793 -0.108866367978112 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 84 octave:797> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:798> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:799> plot(v(:,1),v(:,2),'.-') octave:800> for k=1:size(v,1),fprintf('[%.2f,%.2f],\n',v(k,:));end [-15.00,-0.02], [-14.72,-0.01], [-14.44,0.01], [-14.16,0.04], [-13.88,0.08], [-13.59,0.14], [-13.31,0.21], [-13.03,0.29], [-12.75,0.40], [-12.61,0.46], [-12.47,0.52], [-12.33,0.58], [-12.19,0.66], [-12.05,0.73], [-11.91,0.81], [-11.77,0.90], [-11.62,0.99], [-11.48,1.08], [-11.34,1.19], [-11.20,1.29], [-11.06,1.40], [-10.92,1.52], [-10.78,1.64], [-10.64,1.77], [-10.50,1.91], [-10.36,2.05], [-10.22,2.20], [-10.08,2.35], [-9.94,2.51], [-9.80,2.68], [-9.66,2.85], [-9.52,3.03], [-9.38,3.21], [-9.09,3.61], [-8.81,4.03], [-8.53,4.48], [-8.25,4.96], [-7.97,5.47], [-7.69,6.02], [-7.41,6.59], [-7.12,7.20], [-6.84,7.85], [-6.56,8.53], [-6.00,10.00], [-5.81,10.54], [-5.62,11.13], [-5.25,12.44], [-4.50,15.44], [-3.00,22.06], [-2.25,25.12], [-1.88,26.47], [-1.69,27.08], [-1.50,27.65], [-1.31,28.16], [-1.22,28.40], [-1.12,28.63], [-1.03,28.83], [-0.94,29.03], [-0.80,29.29], [-0.70,29.44], [-0.61,29.58], [-0.49,29.72], [-0.40,29.81], [-0.20,29.95], [-0.01,30.00], [0.00,30.00], [0.16,29.94], [0.31,29.77], [0.39,29.64], [0.47,29.48], [0.55,29.30], [0.62,29.10], [0.70,28.87], [0.78,28.62], [0.86,28.34], [0.94,28.04], [1.09,27.38], [1.25,26.64], [1.56,24.95], [1.88,23.00], [2.50,18.54], [3.75,8.68], [4.38,4.00], [5.00,-0.00], [6.76,-10.00], octave:801> disp(horzcat(p,m')); -15.000000000000000 -0.020000000000000 0.500000000000000 -6.000000000000000 10.000000000000000 70.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -80.000000000000000 octave:802> %for k=1:size(v,1),fprintf('[%.2f,%.2f],\n',v(k,:));end octave:802> system tcsh ans = 0 octave:803> p p = -15.000000000000000 -0.020000000000000 -6.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:804> p(2,:)=[-5,20]; octave:805> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -5.00000000000000e+00 2.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.32410688449638614 -0.01247795712754620 2 -5.000000000000000 20.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.1010090322181513 -0.0501009032218151 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 67 octave:806> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; octave:807> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:808> plot(v(:,1),v(:,2),'.-') octave:809> p p = -15.000000000000000 -0.020000000000000 -5.000000000000000 20.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:810> p=[-15,-0.02;-10,5;-5,20;0,30;5,0] p = -15.000000000000000 -0.020000000000000 -10.000000000000000 5.000000000000000 -5.000000000000000 20.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:811> m m = Columns 1 through 3: 0.500000000000000 70.000000000000000 0.000000000000000 Column 4: -80.000000000000000 octave:812> m=[.5,20,70,0,-80] m = Columns 1 through 3: 0.500000000000000 20.000000000000000 70.000000000000000 Columns 4 and 5: 0.000000000000000 -80.000000000000000 octave:813> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 5.00000000000000e+00 3.63970234266202e-01 0.00872686779075879 0.52611520603045614 -0.06541211591772154 2 -10.000000000000000 5.000000000000000 0.363970234266202 -5.000000000000000 20.000000000000000 2.747477419454622 0.363970234266202 1.104916422402595 -0.115542093851167 3 -5.000000000000000 20.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.1010090322181513 -0.0501009032218151 4 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 85 octave:814> v(end+1,:)=[p(end,1)+(10/cosd(90+m(4)))*sind(90+m(4)),-10]; warning: division by zero octave:815> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 5.00000000000000e+00 3.63970234266202e-01 0.00872686779075879 0.52611520603045614 -0.06541211591772154 2 -10.000000000000000 5.000000000000000 0.363970234266202 -5.000000000000000 20.000000000000000 2.747477419454622 0.363970234266202 1.104916422402595 -0.115542093851167 3 -5.000000000000000 20.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.1010090322181513 -0.0501009032218151 4 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 85 octave:816> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:817> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:818> plot(v(:,1),v(:,2),'.-') octave:819> p p = -15.000000000000000 -0.020000000000000 -10.000000000000000 5.000000000000000 -5.000000000000000 20.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:820> p=[-15,-0.02;-10,5;0,30;5,0] p = -15.000000000000000 -0.020000000000000 -10.000000000000000 5.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:821> m=[.5,50,0,-80] m = Columns 1 through 3: 0.500000000000000 50.000000000000000 0.000000000000000 Column 4: -80.000000000000000 octave:822> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 5.00000000000000e+00 1.19175359259421e+00 0.00872686779075879 0.36055853436485458 -0.03230078158460124 2 -10.000000000000000 5.000000000000000 1.191753592594210 0.000000000000000 30.000000000000000 0.000000000000000 1.1917535925942100 0.5116492814811581 -0.0380824640740579 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 72 octave:823> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:824> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:825> plot(v(:,1),v(:,2),'.-') octave:826> disp(horzcat(p,m')); -15.000000000000000 -0.020000000000000 0.500000000000000 -10.000000000000000 5.000000000000000 50.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -80.000000000000000 octave:827> m m = Columns 1 through 3: 0.500000000000000 50.000000000000000 0.000000000000000 Column 4: -80.000000000000000 octave:828> m(2)=60; octave:829> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 5.00000000000000e+00 1.73205080756888e+00 0.00872686779075879 0.25249909136992121 -0.01068889298561457 2 -10.000000000000000 5.000000000000000 1.732050807568877 0.000000000000000 30.000000000000000 0.000000000000000 1.7320508075688767 0.4035898384862247 -0.0326794919243112 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 75 octave:830> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:831> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:832> plot(v(:,1),v(:,2),'.-') octave:833> m(2)=70; octave:834> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 5.00000000000000e+00 2.74747741945462e+00 0.00872686779075879 0.04941376899277217 0.02992817148981523 2 -10.000000000000000 5.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.2005045161090757 -0.0225252258054538 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 76 octave:835> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:836> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:837> plot(v(:,1),v(:,2),'.-') octave:838> p p = -15.000000000000000 -0.020000000000000 -10.000000000000000 5.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:839> p(2,:)=[-10,8]; octave:840> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 8.00000000000000e+00 2.74747741945462e+00 0.00872686779075879 0.40941376899277226 -0.01807182851018478 2 -10.000000000000000 8.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.1105045161090756 -0.0165252258054538 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 68 octave:841> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:842> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:843> plot(v(:,1),v(:,2),'.-') octave:844> p(2,:)=[-10,10]; octave:845> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.64941376899277237 -0.05007182851018478 2 -10.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.0505045161090756 -0.0125252258054538 3 0.000000000000000 30.000000000000000 0.000000000000000 5.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -2.465743636076459 0.253148727215292 ans = 68 octave:846> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:847> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:848> plot(v(:,1),v(:,2),'.-') octave:849> p p = -15.000000000000000 -0.020000000000000 -10.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 5.000000000000000 0.000000000000000 octave:850> p(4,1)=7 p = -15.000000000000000 -0.020000000000000 -10.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 7.000000000000000 0.000000000000000 octave:851> v=sampleSpline(p,m); size(v,1) 1 -1.50000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.64941376899277237 -0.05007182851018478 2 -10.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.0505045161090756 -0.0125252258054538 3 0.000000000000000 30.000000000000000 0.000000000000000 7.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -1.026551576789307 0.059186668404303 ans = 69 octave:852> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:853> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:854> plot(v(:,1),v(:,2),'.-') octave:855> disp(horzcat(p,m')); -15.000000000000000 -0.020000000000000 0.500000000000000 -10.000000000000000 10.000000000000000 70.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 7.000000000000000 0.000000000000000 -80.000000000000000 octave:856> for k=1:size(v,1),fprintf('[%.2f,%.2f],\n',v(k,:));end [-15.00,-0.02], [-14.84,-0.00], [-14.69,0.04], [-14.53,0.12], [-14.38,0.23], [-14.22,0.36], [-14.06,0.52], [-13.91,0.70], [-13.75,0.91], [-13.59,1.14], [-13.44,1.39], [-13.28,1.66], [-13.12,1.95], [-12.97,2.26], [-12.81,2.58], [-12.50,3.28], [-12.19,4.03], [-11.88,4.82], [-11.25,6.50], [-10.00,10.00], [-7.50,16.99], [-6.25,20.35], [-5.62,21.94], [-5.00,23.43], [-4.38,24.82], [-4.06,25.47], [-3.75,26.09], [-3.44,26.67], [-3.12,27.21], [-2.97,27.46], [-2.81,27.71], [-2.66,27.94], [-2.50,28.16], [-2.34,28.37], [-2.19,28.57], [-2.03,28.76], [-1.88,28.94], [-1.80,29.02], [-1.64,29.18], [-1.48,29.32], [-1.33,29.46], [-1.17,29.57], [-1.02,29.68], [-0.86,29.77], [-0.66,29.86], [-0.51,29.92], [-0.35,29.96], [-0.20,29.99], [-0.04,30.00], [0.00,30.00], [0.22,29.95], [0.44,29.81], [0.55,29.70], [0.66,29.57], [0.77,29.42], [0.88,29.25], [0.98,29.06], [1.09,28.85], [1.20,28.62], [1.31,28.37], [1.42,28.09], [1.53,27.81], [1.75,27.17], [1.97,26.47], [2.19,25.71], [2.62,24.00], [3.06,22.07], [3.50,19.96], [7.00,0.00], [8.76,-10.00], octave:857> p p = -15.000000000000000 -0.020000000000000 -10.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 7.000000000000000 0.000000000000000 octave:858> p(1,1)=-19 p = -19.000000000000000 -0.020000000000000 -10.000000000000000 10.000000000000000 0.000000000000000 30.000000000000000 7.000000000000000 0.000000000000000 octave:859> v=sampleSpline(p,m); size(v,1) 1 -1.90000000000000e+01 -2.00000000000000e-02 8.72686779075879e-03 -1.00000000000000e+01 1.00000000000000e+01 2.74747741945462e+00 0.00872686779075879 0.06389653832931785 0.00653750148862610 2 -10.000000000000000 10.000000000000000 2.747477419454622 0.000000000000000 30.000000000000000 0.000000000000000 2.7474774194546216 0.0505045161090756 -0.0125252258054538 3 0.000000000000000 30.000000000000000 0.000000000000000 7.000000000000000 0.000000000000000 -5.671281819617707 0.000000000000000 -1.026551576789307 0.059186668404303 ans = 93 octave:860> v(end+1,:)=[p(end,1)+(10/cosd(90+m(end)))*sind(90+m(end)),-10]; octave:861> hold off;plot(p(:,1),p(:,2),'r');grid on;hold on;axis equal octave:862> plot(v(:,1),v(:,2),'.-') octave:863> disp(horzcat(p,m')); -19.000000000000000 -0.020000000000000 0.500000000000000 -10.000000000000000 10.000000000000000 70.000000000000000 0.000000000000000 30.000000000000000 0.000000000000000 7.000000000000000 0.000000000000000 -80.000000000000000 octave:864> for k=1:size(v,1),fprintf('[%.2f,%.2f],\n',v(k,:));end [-19.00,-0.02], [-18.72,-0.01], [-18.44,0.01], [-18.16,0.04], [-17.88,0.08], [-17.59,0.14], [-17.31,0.21], [-17.03,0.29], [-16.75,0.40], [-16.61,0.46], [-16.47,0.52], [-16.33,0.58], [-16.19,0.66], [-16.05,0.73], [-15.91,0.81], [-15.77,0.90], [-15.62,0.99], [-15.48,1.08], [-15.34,1.19], [-15.20,1.29], [-15.06,1.40], [-14.92,1.52], [-14.78,1.64], [-14.64,1.77], [-14.50,1.91], [-14.36,2.05], [-14.22,2.20], [-14.08,2.35], [-13.94,2.51], [-13.80,2.68], [-13.66,2.85], [-13.52,3.03], [-13.38,3.21], [-13.09,3.61], [-12.81,4.03], [-12.53,4.48], [-12.25,4.96], [-11.97,5.47], [-11.69,6.02], [-11.41,6.59], [-11.12,7.20], [-10.84,7.85], [-10.56,8.53], [-10.00,10.00], [-7.50,16.99], [-6.25,20.35], [-5.62,21.94], [-5.00,23.43], [-4.38,24.82], [-4.06,25.47], [-3.75,26.09], [-3.44,26.67], [-3.12,27.21], [-2.97,27.46], [-2.81,27.71], [-2.66,27.94], [-2.50,28.16], [-2.34,28.37], [-2.19,28.57], [-2.03,28.76], [-1.88,28.94], [-1.80,29.02], [-1.64,29.18], [-1.48,29.32], [-1.33,29.46], [-1.17,29.57], [-1.02,29.68], [-0.86,29.77], [-0.66,29.86], [-0.51,29.92], [-0.35,29.96], [-0.20,29.99], [-0.04,30.00], [0.00,30.00], [0.22,29.95], [0.44,29.81], [0.55,29.70], [0.66,29.57], [0.77,29.42], [0.88,29.25], [0.98,29.06], [1.09,28.85], [1.20,28.62], [1.31,28.37], [1.42,28.09], [1.53,27.81], [1.75,27.17], [1.97,26.47], [2.19,25.71], [2.62,24.00], [3.06,22.07], [3.50,19.96], [7.00,0.00], [8.76,-10.00], octave:865> system tcsh ans = 0 octave:866> diary off