Descrição2D 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 fevereiro de 2007
Origem	Own drawing, plotted in MuPAD.
Autor	Inductiveload
Permissão (Reutilizar este ficheiro)	Public domainPublic domainfalsefalse Eu, titular dos direitos de autor desta obra, dedico-a ao domínio público, com aplicação em todo o mundo. Nalguns países isto pode não ser legalmente possível; se assim for: Concedo a todos o direito de usar esta obra para qualquer fim, sem quaisquer condições, a menos que tais condições sejam impostas por lei.
Outras versões	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)