Imagen:2D Wavefunction (1,2) Surface Plot.png

Description	The quantum wavefunction of a particle in a 2D box of dimensions Lx and Ly. The wavenumbers are: nx=1 ny=2
Source	Own drawing, plotted in MuPAD.
Date	13/02/2007
Author	Inductiveload
Permission (Reusing this image)	I, the copyright holder of this work, hereby release it into the public domain. This applies worldwide. In case this is not legally possible: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law. Afrikaans | Alemannisch | Aragonés | العربية | Asturianu | Български | Català | Česky | Cymraeg | Dansk | Deutsch | Eʋegbe | Ελληνικά | English | Español | Esperanto | Euskara | Estremeñu | فارسی | Français | Galego | 한국어 | हिन्दी | Hrvatski | Ido | Bahasa Indonesia | Íslenska | Italiano | עברית | Kurdî / كوردی | Latina | Lietuvių | Latviešu | Magyar | Македонски | Bahasa Melayu | Nederlands | ‪Norsk (bokmål)‬ | ‪Norsk (nynorsk)‬ | 日本語 | Polski | Português | Ripoarisch | Română | Русский | Shqip | Slovenčina | Slovenščina | Српски / Srpski | Svenska | ไทย | Tagalog | Türkçe | Українська | Tiếng Việt | Walon | ‪中文(简体)‬ | ‪中文(繁體)‬ | zh-yue-hant | +/-
Other versions	See 2D Quantum Wavefunctions for more plots.

[edit] Mathematica Code

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}
  ]

[edit] MuPAD Code

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)

Las siguientes páginas enlazan a esta imagen:

• Función de ondas
• Partícula en una caja