# set terminal canvas solid butt size 600,400 fsize 10 lw 1 fontscale 1 name "surface1_17" jsdir "." # set output 'surface1.17.js' set dummy u,v set label 1 "increasing v" at 6, 0, -1 left norotate back nopoint offset character 0, 0, 0 set label 2 "increasing u" at 0, -5, -1.5 left norotate back nopoint offset character 0, 0, 0 set label 3 "floor(u)%3=0" at 5, 6.5, -1 left norotate back nopoint offset character 0, 0, 0 set label 4 "floor(u)%3=1" at 5, 6.5, 0.100248 left norotate back nopoint offset character 0, 0, 0 set arrow 1 from 5, -5, -1.2 to 5, 5, -1.2 head back nofilled linetype -1 linewidth 1.000 set arrow 2 from -5, -5, -1.2 to 5, -5, -1.2 head back nofilled linetype -1 linewidth 1.000 set arrow 3 from 5, 6, -1 to 5, 5, -1 head back nofilled linetype -1 linewidth 1.000 set arrow 4 from 5, 6, 0.100248 to 5, 5, 0.100248 head back nofilled linetype -1 linewidth 1.000 set parametric set view 70, 20, 1, 1 set samples 51, 51 set isosamples 30, 33 set hidden3d back offset 1 trianglepattern 3 undefined 1 altdiagonal bentover set ztics -1.00000,0.25,1.00000 norangelimit set title "\"fence plot\" using single parametric surface with undefined points" set xlabel "X axis" set xlabel offset character -3, -2, 0 font "" textcolor lt -1 norotate set xrange [ -5.00000 : 5.00000 ] noreverse nowriteback set ylabel "Y axis" set ylabel offset character 3, -2, 0 font "" textcolor lt -1 rotate by -270 set yrange [ -5.00000 : 5.00000 ] noreverse nowriteback set zlabel "Z axis" set zlabel offset character -5, 0, 0 font "" textcolor lt -1 norotate set zrange [ -1.00000 : 1.00000 ] noreverse nowriteback sinc(u,v) = sin(sqrt(u**2+v**2)) / sqrt(u**2+v**2) GPFUN_sinc = "sinc(u,v) = sin(sqrt(u**2+v**2)) / sqrt(u**2+v**2)" xx = 6.08888888888889 dx = 1.11 x0 = -5 x1 = -3.89111111111111 x2 = -2.78222222222222 x3 = -1.67333333333333 x4 = -0.564444444444444 x5 = 0.544444444444445 x6 = 1.65333333333333 x7 = 2.76222222222222 x8 = 3.87111111111111 x9 = 4.98 xmin = -4.99 xmax = 5 n = 10 zbase = -1 splot [u=.5:3*n-.5][v=-4.99:4.99] xmin+floor(u/3)*dx, v, ((floor(u)%3)==0) ? zbase : (((floor(u)%3)==1) ? sinc(xmin+u/3.*dx,v) : 1/0) notitle