Descripció2D Wavefunction (1,2) Surface Plot.png	The quantum wavefunction of a particle in a 2D box of dimensions Lx and Ly. The wavenumbers are: nx=1 ny=2
Data	13 de febrer de 2007
Font	Own drawing, plotted in MuPAD.
Autor	Inductiveload
Permís (Com reutilitzar aquest fitxer)	Public domainPublic domainfalsefalse Jo, el titular del copyright d'aquesta obra, l'allibero al domini públic. Això s'aplica a tot el món. En alguns països això pot no ser legalment possible, en tal cas: Jo faig concessió a tothom del dret d'usar aquesta obra per a qualsevol propòsit, sense cap condició llevat d'aquelles requerides per la llei.
Altres versions	See 2D Quantum Wavefunctions for more plots.

This code does not require any packages to be loaded.

nx = 1
ny = 2
Plot3D[
  Sin[nx Pi x]Sin[ny Pi y],
  {x, 0, 1}, 
  {y, 0, 1},
  ImageSize -> 2500,
  BoxRatios -> {2, 2, 2},
  AxesLabel -> {x/Subscript[L, x], y/Subscript[L, y], Subscript[ψ, "1,1"]},
  PlotRange -> {-1, 1},
  TextStyle -> {FontSize -> 96}
  ]

nx := 1:
ny := 2: 
f := sin(nx*x*PI) * sin(ny*y*PI):

  conts := 16:
  projectionlevel:=-1.5:

colourfunc := zip(RGB::Red, RGB::Yellow,
                  (a, b) -> (f(x,y)+1)/(2)*a
                           -(f(x,y)-1)/(2)*b):

funcplot := plot::Function3d(f(x,y),
                             x = 0..1,
                             y = 0..1,
                             Mesh = [24, 24],
                             Submesh=[8,8],
                             LineColor = RGB::Black.[0.99],
                             LineWidth = 0.2，
                             FillColorFunction = colourfunc,                             
                             AxesTitleFont = ["BKM-cmr17", 20],
                             TicksLabelFont = ["Courier New", 12],                         
                             ViewingBoxZRange = -1.5..1
                            ):

contours := plot::modify(funcplot,
                         ZContours = [Automatic, conts],
                         LineWidth = 0.2,
                         LineColor = RGB::Gray90.[0.5],
                         XLinesVisible = FALSE,
                         YLinesVisible = FALSE,
                         Filled = FALSE
                        ):
                        
projcontours := plot::Transform3d([0, 0, projectionlevel],
                                  [1, 0, 0, 0, 1, 0, 0, 0, 0], 
                                   plot::modify(funcplot, ZContours = [Automatic, conts],
                                                LineWidth = 0.2,
                                                LineColorType = Dichromatic,
                                                LineColor = RGB::Red.[0.99],
                                                LineColor2 = RGB::Yellow.[0.99],
                                                XLinesVisible = FALSE,
                                                YLinesVisible = FALSE,
                                                Filled = FALSE
                                               )
                                  ):
                             
camera := plot::Camera([17, -24, 45],
                       [0.5, 0.4, -0.2],
                       0.038
                      ):
                             
ploteverything := plot::Canvas(funcplot, contours, projcontours,
                               Width = 200*unit::mm, 
                               Height = 150*unit::mm,
                               AxesTitles = ["x/L(x)", "y/L(y)", "z"],                                                              
                               camera
                              ): 
                          
plot(ploteverything)