2011-11-07 17 views
4

私は軌道2PZ水素のxyスライスをプロットするために働く次のコードを持っている:MATLAB - 3D等高線図(水素2P軌道)

%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(:,:,:)); %3D contour 

これをどうやって動作させることができますか?

答えて

6

基本グリッドと必要なレベルを指定するだけで済みます。たとえば、次のように

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

enter image description here

関連する問題