Блог

У меня есть следующий код, который работает для построения среза x-y орбитали водорода 2pz:

 %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
  

Я хотел бы построить 3D-контурный график орбиты. Я пробовал это (и, похоже, это не дает мне того, что я хотел):

 isosurface(psi_210_vals(:,:,:)); = contour
  

Как я могу заставить это работать?

Вам просто нужно указать базовую сетку и желаемый уровень. Например:

 >> isosurface(-10:.1:10, -10:.1:10, -2:.1:2, psi_210_vals, 0.001);
>> axis equal