我有以下代码,可以绘制氢原子2pz轨道的x-y切片:
我该如何让这个运作起来?
%probability densities
pd_psi_210 = @(r,theta,phi) exp(-r).*(r.^2).*(cos(theta).^2)/(32*pi);
%configuring the range
[x y z] = meshgrid(-10:.1:10,-10:.1:10,-2:.1:2);
[THETA,PHI,R] = cart2sph(x,y,z);
%create array of probability density magnitudes
psi_210_vals = pd_psi_210(R,THETA,PHI);
%plotting
imagesc(psi_210_vals(:,:,1)); %x-y plane
我想绘制轨道的三维等高线图。我试过这个方法(但好像不是我想要的结果):
isosurface(psi_210_vals(:,:,:)); %3D contour
我该如何让这个运作起来?