Return to singulr.dem CVS log

Up to [local] / OpenXM_contrib / gnuplot / demo

Annotation of OpenXM_contrib/gnuplot/demo/singulr.dem, Revision 1.1.1.2

1.1.1.2 ! maekawa     1: # $Id: singulr.dem,v 1.1.1.1.2.2 1999/08/19 14:35:07 lhecking Exp $
1.1       maekawa     2: #
                      3: # Demo that plots some surfaces with singularities.
                      4: # Author: Carsten Steger
                      5: #
                      6: # (x,y,x^2-y^2,2xy) is the graph of w=z^2 in 4-space.
                      7: # Therefore (x^2-y^2,2xy,x,y) is the graph of w=sqrt(z) in 4-space.
                      8: # Coordinates 1, 2, and 3 give the real part of either function,
                      9: # whereas coordinates 1, 2, and 4 give the imaginary part.
                     10: # The same holds for the cube function w=z^3. The graphs are given by
                     11: # (x,y,x^3-3xy^2,3x^2y-y^3) and (x^3-3xy^2,3x^2y-y^3,x,y).
                     12: # And so on...
                     13:
                     14:
                     15: set parametric
                     16: set hidden3d
                     17: set isosamples 21
                     18: set autoscale
                     19:
                     20:
                     21: set view 60,30
                     22: set urange [-3:3]
                     23: set vrange [-3:3]
                     24: set title "Real part of complex square root function"
                     25: splot u**2-v**2,2*u*v,u
                     26: pause -1 "Hit return to continue (1)"
                     27:
                     28:
                     29: set view 60,210
                     30: set title "Real part of complex square root function (different view)"
                     31: replot
                     32: pause -1 "Hit return to continue (2)"
                     33:
                     34:
                     35: set view 60,120
                     36: set urange [-3:3]
                     37: set vrange [-3:3]
                     38: set title "Imaginary part of complex square root function"
                     39: splot u**2-v**2,2*u*v,v
                     40: pause -1 "Hit return to continue (3)"
                     41:
                     42:
                     43: set view 60,300
                     44: set title "Imaginary part of complex square root function (different view)"
                     45: replot
                     46: pause -1 "Hit return to continue (4)"
                     47:
                     48:
                     49: set view 60,30
                     50: set urange [-3:3]
                     51: set vrange [-3:3]
                     52: set title "Real part of complex cube root function"
                     53: splot u**3-3*u*v**2,3*u**2*v-v**3,u
                     54: pause -1 "Hit return to continue (5)"
                     55:
                     56:
                     57: set view 60,210
                     58: set title "Real part of complex cube root function (different view)"
                     59: replot
                     60: pause -1 "Hit return to continue (6)"
                     61:
                     62:
                     63: set view 60,30
                     64: set urange [-3:3]
                     65: set vrange [-3:3]
                     66: set title "Imaginary part of complex cube root function"
                     67: splot u**3-3*u*v**2,3*u**2*v-v**3,v
                     68: pause -1 "Hit return to continue (7)"
                     69:
                     70:
                     71: set view 60,210
                     72: set title "Imaginary part of complex cube root function (different view)"
                     73: replot
                     74: pause -1 "Hit return to continue (8)"
                     75:
                     76:
                     77: set view 60,30
                     78: set isosamples 31
                     79: set urange [-1:1]
                     80: set vrange [-1:1]
                     81: set title "Real part of complex 4th root function"
                     82: splot u**4-6*u**2*v**2+v**4,4*u**3*v-4*u*v**3,u
                     83: pause -1 "Hit return to continue (9)"
                     84:
                     85:
                     86: set view 60,210
                     87: set title "Real part of complex 4th root function (different view)"
                     88: replot
                     89: pause -1 "Hit return to continue (10)"
                     90:
                     91:
                     92: set view 60,120
                     93: set urange [-1:1]
                     94: set vrange [-1:1]
                     95: set title "Imaginary part of complex 4th root function"
                     96: splot u**4-6*u**2*v**2+v**4,4*u**3*v-4*u*v**3,v
                     97: pause -1 "Hit return to continue (11)"
                     98:
                     99:
                    100: set view 60,300
                    101: set title "Imaginary part of complex 4th root function (different view)"
                    102: replot
                    103: pause -1 "Hit return to continue (12)"
                    104:
                    105:
                    106: set isosamples 21
                    107: set view 60,20
                    108: set urange [-3:3]
                    109: set vrange [-3:3]
                    110: set title "Enneper's surface"
                    111: splot u-u**3/3+u*v**2,v-v**3/3+v*u**2,u**2-v**2
                    112: pause -1 "Hit return to continue (13)"
                    113:
                    114:
                    115: set view 60,110
                    116: set title "Enneper's surface (different view)"
                    117: replot
                    118: pause -1 "Hit return to continue (14)"
                    119:
                    120:
                    121: set isosamples 31,11
                    122: set view 60,30
                    123: set title "Moebius strip"
                    124: set urange [0:2*pi]
                    125: set vrange [-0.25:0.25]
                    126: splot (2-v*sin(u/2))*sin(u),(2-v*sin(u/2))*cos(u),v*cos(u/2)
                    127: pause -1 "Hit return to continue (15)"
                    128:
                    129:
                    130: set view 60,210
                    131: set title "Moebius strip (view from opposite side)"
                    132: replot
                    133: pause -1 "Hit return to continue (16)"
                    134:
                    135: set nokey
                    136: set xrange [-10:10]
                    137: set yrange [-10:10]
                    138: set zrange [-3:3]
                    139: set urange [0:2*pi]
                    140: set vrange [0:2*pi]
                    141: set isosamples 39,60
                    142: set view 60,120
                    143: set title "Klein bottle"
                    144: splot  (2*sin(u)*cos(v/2)-sin(2*u)*sin(v/2)+8)*cos(v), \
                    145:        (2*sin(u)*cos(v/2)-sin(2*u)*sin(v/2)+8)*sin(v), \
                    146:         2*sin(u)*sin(v/2)+sin(2*u)*cos(v/2)
                    147: pause -1 "Hit return to continue (17)"
                    148:
                    149:
                    150: set urange [0:2*pi]
                    151: set vrange [0:4*pi/3]
                    152: set isosamples 39,40
                    153: set view 60,20
                    154: set title "Klein bottle with look at the 'inside'"
                    155: replot
                    156: pause -1 "Hit return to continue (18)"
                    157:
                    158: set data style lines
                    159: set xrange [-12:12]
                    160: set yrange [-12:12]
                    161: set zrange [*:*]
                    162: set nohidden3d
1.1.1.2 ! maekawa   163: set ticslevel 0.
        !           164: set view 50,15 #HBB: ,1,1.7
1.1       maekawa   165: set title "Klein bottle, glassblowers' version (look-through)"
                    166: splot "klein.dat"
                    167: pause -1 "Hit return to continue (19)"
                    168:
                    169:
                    170: set hidden3d
                    171: set view 70,305
                    172: set title "Klein bottle, glassblowers' version (solid)"
                    173: splot "klein.dat"
                    174: pause -1 "Hit return to continue (20)"
                    175:
                    176: reset