Wolfram Cloud Document

L = 1;
nX = 1; nY = 1; nZ = 1;
psi[x_, y_, z_] := Sin[nX Pi x/L] Sin[nY Pi y/L] Sin[nZ Pi z/L]

ContourPlot3D[
psi[x, y, z] == 0.5,
{x, 0, L}, {y, 0, L}, {z, 0, L},
Mesh -> None, ContourStyle -> Opacity[0.5, Blue]
]

Out[]=

(*Parámetrosdelpozo*)L=1;(*Longituddelpozo*)nMax=4;(*Númeromáximodeestadosavisualizar*)m=1;(*Masadelapartícula*)ħ=1;(*ConstantedePlanckreducida*)(*Definirlafuncióndeondaestacionaria*)psi[n_,x_]:=Sqrt[2/L]Sin[nPix/L](*Definirladensidaddeprobabilidad*)probDensity[n_,x_]:=Abs[psi[n,x]]^2(*Graficarlasfuncionesdeonda*)Plot[Evaluate[Table[psi[n,x],{n,1,nMax}]],{x,0,L},PlotLegends"Expressions",PlotRangeAll,AxesLabel{"x","ψ_n(x)"},PlotLabel"Funciones de Onda Estacionarias en 1D",FillingAxis](*Graficarlasdensidadesdeprobabilidad*)Plot[Evaluate[Table[probDensity[n,x],{n,1,nMax}]],{x,0,L},PlotLegends"Expressions",PlotRangeAll,AxesLabel{"x","|ψ_n(x)|^2"},PlotLabel"Densidades de Probabilidad en 1D",FillingAxis]

Out[]=

	2 sin(πx)
	2 sin(2πx)
	2 sin(3πx)
	2 sin(4πx)

Out[]=

	2 2 sin(πx)
	2 2 sin(2πx)
	2 2 sin(3πx)
	2 2 sin(4πx)

In[]:=

(*Parámetrosadicionales*)E[n_]:=(n^2Pi^2ħ^2)/(2mL^2)(*Nivelesdeenergía*)(*Funcióndeondadependientedeltiempo*)Psi[n_,x_,t_]:=psi[n,x]Exp[-IE[n]t/ħ](*Graficarlaparterealdelafuncióndeondaendiferentestiempos*)Animate[Plot[Re[Psi[n,x,t]],{x,0,L},PlotRange{-Sqrt[2/L],Sqrt[2/L]},AxesLabel{"x","Re[ψ_n(x,t)]"},PlotLabelStringForm["Función de Onda Real para n=`1`, t=`2`",n,t]],{t,0,4mL^2/(Pi^2ħ^2),0.1}]

SetDelayed

:Tag  in [n_] is Protected.

Out[]=

$Failed

Animate

:Animate argument t,0,

4m

2

L

2

π

2

ħ

,0.1 does not have the correct form for a variable specification.

Out[]=

AnimatePlotRe[Psi[n,x,t]],{x,0,L},PlotRange-

2

L

,

2

L

,AxesLabel{x,Re[ψ_n(x,t)]},PlotLabelFunción de Onda Real para n=n, t=t,t,0,

4m

2

L

2

π

2

ħ

,0.1

In[]:=

(*Parámetrosdelpozo*)L=1;(*Longituddelpozoenambasdirecciones*)nMax=2;(*Númeromáximodeestadosencadadirección*)m=1;(*Masadelapartícula*)ħ=1;(*ConstantedePlanckreducida*)(*Definirlafuncióndeondaestacionariaen2D*)psi2D[nx_,ny_,x_,y_]:=Sqrt[4/L^2]Sin[nxPix/L]Sin[nyPiy/L](*Definirladensidaddeprobabilidaden2D*)probDensity2D[nx_,ny_,x_,y_]:=Abs[psi2D[nx,ny,x,y]]^2(*Crearunamalladecoordenadas*)resolution=50;x=Linspace[0,L,resolution];y=Linspace[0,L,resolution];data=Table[{x[[i]],y[[j]],probDensity2D[1,1,x[[i]],y[[j]]]},{i,1,resolution},{j,1,resolution}];data=Flatten[data,1];(*GraficarladensidaddeprobabilidadutilizandoPlot3D*)Plot3D[psi2D[1,1,x,y],{x,0,L},{y,0,L},PlotRangeAll,AxesLabel{"x","y","ψ_{1,1}(x,y)"},PlotLabel"Función de Onda Estacionaria 2D para n_x=n_y=1",ColorFunction"Rainbow",MeshNone](*Graficarladensidaddeprobabilidad*)DensityPlot[probDensity2D[1,1,x,y],{x,0,L},{y,0,L},PlotRangeAll,ColorFunction"Rainbow",AxesLabel{"x","y","|ψ_{1,1}(x,y)|^2"},PlotLabel"Densidad de Probabilidad 2D para n_x=n_y=1"](*Graficarnivelesdeenergía*)E2D[nx_,ny_]:=(Pi^2ħ^2)/(2mL^2)(nx^2+ny^2)(*Mostrarlosnivelesdeenergía*)Table[E2D[nx,ny],{nx,1,nMax},{ny,1,nMax}]