例子：网格50*50，网格距为1，每个网格里有均匀分布的100个粒子，粒子电量为1，使用粒子云分室法计算格点的电荷量。

结果图

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