Lx=1
Ly=1
@interact
def _(Nx=(1..9), Ny=(1..9)):
psi(x,y)=sqrt(4/(Lx*Ly))*sin(Nx*pi/Lx*x)*sin(Ny*pi/Ly*y)
Pnorm=plot3d(psi(x,y),(x,0,Lx),(y,0,Ly),adaptive=True, color='red')
Pswap=plot3d(psi(y,x),(x,0,Ly),(y,0,Lx),adaptive=True, color='blue')
Pasym=plot3d(psi(x,y)+psi(y,x),(x,0,Ly),(y,0,Lx),adaptive=True, color='green')
html('<font color="red">$\psi(x,y)\;=\;%s$'%latex(psi(x,y)))
html('</font><font color="blue">$\psi(y,x)\;=\;%s$'%latex(psi(y,x)))
html('</font><font color="green">$(\psi(x,y)-\psi(y,x))\;=\;%s$'%latex(psi(x,y)-psi(y,x)))
Pnorm.save('2dpsinorm1.png')
Pswap.save('2dpsiswap1.png')
Pasym.save('2dpsiasym.png')