octave:7> p = [0,2; 1,2; 2,.5; 4,-.1] p = 0.00000 2.00000 1.00000 2.00000 2.00000 0.50000 4.00000 -0.10000 octave:8> slope=[0 0 -70 -10]; octave:9> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 0.50000 -2.74748 0.00000 -1.75252 0.25252 3 2.00000 0.50000 -2.74748 4.00000 -0.10000 -0.17633 -2.74748 2.38564 -0.58095 octave:10> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:11> p = [0,2; 1,2; 2,1; 4,-.1] p = 0.00000 2.00000 1.00000 2.00000 2.00000 1.00000 4.00000 -0.10000 octave:12> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 1.00000 -2.74748 0.00000 -0.25252 -0.74748 3 2.00000 1.00000 -2.74748 4.00000 -0.10000 -0.17633 -2.74748 2.01064 -0.45595 octave:13> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:14> p = [0,2; 1,2; 2,1; 2.5,.3;4,-.1] p = 0.00000 2.00000 1.00000 2.00000 2.00000 1.00000 2.50000 0.30000 4.00000 -0.10000 octave:15> slope=[0 0 -70 -20 -10]; octave:16> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 1.00000 -2.74748 0.00000 -0.25252 -0.74748 3 2.00000 1.00000 -2.74748 2.50000 0.30000 -0.36397 -2.7475 3.3179 -1.2458 4 2.50000 0.30000 -0.36397 4.00000 -0.10000 -0.17633 -0.3639702 0.0695116 -0.0030951 octave:17> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:18> slope=[0 0 -80 -20 -10]; octave:19> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 1.00000 -5.67128 0.00000 2.67128 -3.67128 3 2.00000 1.00000 -5.67128 2.50000 0.30000 -0.36397 -5.6713 15.0131 -12.9410 4 2.50000 0.30000 -0.36397 4.00000 -0.10000 -0.17633 -0.3639702 0.0695116 -0.0030951 octave:20> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:21> slope=[0 0 -70 -20 -10]; octave:22> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 1.00000 -2.74748 0.00000 -0.25252 -0.74748 3 2.00000 1.00000 -2.74748 2.50000 0.30000 -0.36397 -2.7475 3.3179 -1.2458 4 2.50000 0.30000 -0.36397 4.00000 -0.10000 -0.17633 -0.3639702 0.0695116 -0.0030951 octave:23> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:24> p = [0,2; 1,2; 2,1; 2.6,.2;4,-.1] p = 0.00000 2.00000 1.00000 2.00000 2.00000 1.00000 2.60000 0.20000 4.00000 -0.10000 octave:25> v=sampleSpline(p,slope,2,.1); 1 0 2 0 1 2 0 0 0 0 2 1.00000 2.00000 0.00000 2.00000 1.00000 -2.74748 0.00000 -0.25252 -0.74748 3 2.00000 1.00000 -2.74748 2.60000 0.20000 -0.36397 -2.7475 3.0982 -1.2355 4 2.60000 0.20000 -0.36397 4.00000 -0.10000 -0.17633 -0.363970 0.186722 -0.057003 octave:26> hold off;plot(v(:,1),v(:,2));hold on;plot(p(:,1),p(:,2),'or');grid on octave:27> ls *.m octave:28> whos v Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== v 48x2 768 double Total is 96 elements using 768 bytes octave:29> v v = 0.00000 2.00000 0.50000 2.00000 1.00000 2.00000 1.06250 1.99883 1.12500 1.99459 1.18750 1.98620 1.25000 1.97254 1.31250 1.95253 1.37500 1.92507 1.43750 1.88907 1.50000 1.84343 1.56250 1.78707 1.62500 1.71887 1.65625 1.67999 1.68750 1.63775 1.71875 1.59200 1.75000 1.54261 1.78125 1.48945 1.81250 1.43237 1.84375 1.37123 1.87500 1.30591 1.90625 1.23626 1.93750 1.16215 1.96875 1.08345 2.00000 1.00000 2.03750 0.90126 2.07500 0.81085 2.09375 0.76864 2.13125 0.68997 2.16875 0.61865 2.20625 0.55429 2.24375 0.49649 2.28125 0.44486 2.31875 0.39901 2.35625 0.35856 2.39375 0.32310 2.43125 0.29225 2.50625 0.24282 2.58125 0.20714 2.60000 0.20000 2.68750 0.16954 2.77500 0.14172 2.86250 0.11629 2.95000 0.09304 3.03750 0.07173 3.12500 0.05213 3.30000 0.01716 4.00000 -0.10000 octave:30> ls *.scad octave:31> ls -l *,.scad parse error: syntax error >>> ls -l *,.scad ^ octave:31> ls -l *.scad octave:32> system('less parabola.scad') ans = 0 octave:33> for k=1:length(f),fprintf(1,' [%.2f,%.2f],\n',v(k,:));end error: 'f' undefined near line 1 column 16 octave:33> for k=1:length(v),fprintf(1,' [%.2f,%.2f],\n',v(k,:));end [0.00,2.00], [0.50,2.00], [1.00,2.00], [1.06,2.00], [1.12,1.99], [1.19,1.99], [1.25,1.97], [1.31,1.95], [1.38,1.93], [1.44,1.89], [1.50,1.84], [1.56,1.79], [1.62,1.72], [1.66,1.68], [1.69,1.64], [1.72,1.59], [1.75,1.54], [1.78,1.49], [1.81,1.43], [1.84,1.37], [1.88,1.31], [1.91,1.24], [1.94,1.16], [1.97,1.08], [2.00,1.00], [2.04,0.90], [2.08,0.81], [2.09,0.77], [2.13,0.69], [2.17,0.62], [2.21,0.55], [2.24,0.50], [2.28,0.44], [2.32,0.40], [2.36,0.36], [2.39,0.32], [2.43,0.29], [2.51,0.24], [2.58,0.21], [2.60,0.20], [2.69,0.17], [2.78,0.14], [2.86,0.12], [2.95,0.09], [3.04,0.07], [3.12,0.05], [3.30,0.02], [4.00,-0.10], octave:34> system('head parabola.scad') ans = 0 octave:35> for k=length(v):-1:1,fprintf(1,' [%.2f,%.2f],\n',v(k,:));end [4.00,-0.10], [3.30,0.02], [3.12,0.05], [3.04,0.07], [2.95,0.09], [2.86,0.12], [2.78,0.14], [2.69,0.17], [2.60,0.20], [2.58,0.21], [2.51,0.24], [2.43,0.29], [2.39,0.32], [2.36,0.36], [2.32,0.40], [2.28,0.44], [2.24,0.50], [2.21,0.55], [2.17,0.62], [2.13,0.69], [2.09,0.77], [2.08,0.81], [2.04,0.90], [2.00,1.00], [1.97,1.08], [1.94,1.16], [1.91,1.24], [1.88,1.31], [1.84,1.37], [1.81,1.43], [1.78,1.49], [1.75,1.54], [1.72,1.59], [1.69,1.64], [1.66,1.68], [1.62,1.72], [1.56,1.79], [1.50,1.84], [1.44,1.89], [1.38,1.93], [1.31,1.95], [1.25,1.97], [1.19,1.99], [1.12,1.99], [1.06,2.00], [1.00,2.00], [0.50,2.00], [0.00,2.00], octave:36> for k=length(v):-1:1,fprintf(1,' [%.2f,%.2f],\n',v(k,:) .* [-1,1]);end [-4.00,-0.10], [-3.30,0.02], [-3.12,0.05], [-3.04,0.07], [-2.95,0.09], [-2.86,0.12], [-2.78,0.14], [-2.69,0.17], [-2.60,0.20], [-2.58,0.21], [-2.51,0.24], [-2.43,0.29], [-2.39,0.32], [-2.36,0.36], [-2.32,0.40], [-2.28,0.44], [-2.24,0.50], [-2.21,0.55], [-2.17,0.62], [-2.13,0.69], [-2.09,0.77], [-2.08,0.81], [-2.04,0.90], [-2.00,1.00], [-1.97,1.08], [-1.94,1.16], [-1.91,1.24], [-1.88,1.31], [-1.84,1.37], [-1.81,1.43], [-1.78,1.49], [-1.75,1.54], [-1.72,1.59], [-1.69,1.64], [-1.66,1.68], [-1.62,1.72], [-1.56,1.79], [-1.50,1.84], [-1.44,1.89], [-1.38,1.93], [-1.31,1.95], [-1.25,1.97], [-1.19,1.99], [-1.12,1.99], [-1.06,2.00], [-1.00,2.00], [-0.50,2.00], [-0.00,2.00], octave:37> system('tail parabola.scad') ans = 0 octave:38> for k=1:length(v),fprintf(1,' [%.2f,%.2f],\n',v(k,:)-[1 0]);end [-1.00,2.00], [-0.50,2.00], [0.00,2.00], [0.06,2.00], [0.12,1.99], [0.19,1.99], [0.25,1.97], [0.31,1.95], [0.38,1.93], [0.44,1.89], [0.50,1.84], [0.56,1.79], [0.62,1.72], [0.66,1.68], [0.69,1.64], [0.72,1.59], [0.75,1.54], [0.78,1.49], [0.81,1.43], [0.84,1.37], [0.88,1.31], [0.91,1.24], [0.94,1.16], [0.97,1.08], [1.00,1.00], [1.04,0.90], [1.08,0.81], [1.09,0.77], [1.13,0.69], [1.17,0.62], [1.21,0.55], [1.24,0.50], [1.28,0.44], [1.32,0.40], [1.36,0.36], [1.39,0.32], [1.43,0.29], [1.51,0.24], [1.58,0.21], [1.60,0.20], [1.69,0.17], [1.78,0.14], [1.86,0.12], [1.95,0.09], [2.04,0.07], [2.12,0.05], [2.30,0.02], [3.00,-0.10], octave:39> for k=1:length(v),fprintf(1,' [%.2f,%.2f],\n',v(k,:));end [0.00,2.00], [0.50,2.00], [1.00,2.00], [1.06,2.00], [1.12,1.99], [1.19,1.99], [1.25,1.97], [1.31,1.95], [1.38,1.93], [1.44,1.89], [1.50,1.84], [1.56,1.79], [1.62,1.72], [1.66,1.68], [1.69,1.64], [1.72,1.59], [1.75,1.54], [1.78,1.49], [1.81,1.43], [1.84,1.37], [1.88,1.31], [1.91,1.24], [1.94,1.16], [1.97,1.08], [2.00,1.00], [2.04,0.90], [2.08,0.81], [2.09,0.77], [2.13,0.69], [2.17,0.62], [2.21,0.55], [2.24,0.50], [2.28,0.44], [2.32,0.40], [2.36,0.36], [2.39,0.32], [2.43,0.29], [2.51,0.24], [2.58,0.21], [2.60,0.20], [2.69,0.17], [2.78,0.14], [2.86,0.12], [2.95,0.09], [3.04,0.07], [3.12,0.05], [3.30,0.02], [4.00,-0.10], octave:40> for k=length(v):-1:1,fprintf(1,' [%.2f,%.2f],\n',v(k,:) .* [-1,1]);end [-4.00,-0.10], [-3.30,0.02], [-3.12,0.05], [-3.04,0.07], [-2.95,0.09], [-2.86,0.12], [-2.78,0.14], [-2.69,0.17], [-2.60,0.20], [-2.58,0.21], [-2.51,0.24], [-2.43,0.29], [-2.39,0.32], [-2.36,0.36], [-2.32,0.40], [-2.28,0.44], [-2.24,0.50], [-2.21,0.55], [-2.17,0.62], [-2.13,0.69], [-2.09,0.77], [-2.08,0.81], [-2.04,0.90], [-2.00,1.00], [-1.97,1.08], [-1.94,1.16], [-1.91,1.24], [-1.88,1.31], [-1.84,1.37], [-1.81,1.43], [-1.78,1.49], [-1.75,1.54], [-1.72,1.59], [-1.69,1.64], [-1.66,1.68], [-1.62,1.72], [-1.56,1.79], [-1.50,1.84], [-1.44,1.89], [-1.38,1.93], [-1.31,1.95], [-1.25,1.97], [-1.19,1.99], [-1.12,1.99], [-1.06,2.00], [-1.00,2.00], [-0.50,2.00], [-0.00,2.00], octave:41> exit