|||
结果图
MATLAB程序
function [q]=electron_density()
q=zeros(51,51);
for i=1:50
for j=1:50
for m=1:10
for n=1:10
q(i,j)=q(i,j)+(11-m)*(11-n)/(11*11);
q(i+1,j)=q(i+1,j)+m*(11-n)/(11*11);
q(i,j+1)=q(i,j+1)+(11-m)*n/(11*11);
q(i+1,j+1)=q(i+1,j+1)+m*n/(11*11);
end
end
end
end
raw=zeros(1,51);column=zeros(52,1)
q_raw=[q;raw];q_column=[q_raw column]
figure(1);pcolor(q_column);colorbar;
end
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-21 11:45
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社