Simulation of sewage discharge

clear
length = 200;
width = 100;
u=2;E=0.36;M=30; %定义计算边界条件
t=100; %时间限制
x=1:1:length;
y=1:1:width;
x0=0;%污染源头位置
y0=50;

figure                                 % 创建图窗
myVideo=VideoWriter('VideoFile_test2'); % 指定视频文件名

myVideo.FrameRate=5;                  %调整视频帧率
open(myVideo)                          % open开始写入
for t=1:1:120
    for i=1:length
        for j=1:width
            C=@(t1)(2.*M./(4.*pi.*E.*(t-t1)).*exp((-(i-u.*(t-t1)).^2-(j-y0).^2)./(4.*E.*(t-t1))));  % 均匀流中的连续源扩散解析解
            c(i,j)=integral(C,0,t);
        end
    end
    surfc(y,x,c);%产生由X、Y、c指定的有色参数化曲面，即三维有色图。
    xlabel('X/m');zlabel('Y/m');%加X、Y标题X/m、Y/m
    shading flat;%去掉图上的黑色网格线
    title('污染物连续点源浓度随空间的变化');%加图标题
    colorbar;%给图形添加色标
    drawnow                           % 即时显示当前图窗
    Fr=getframe(gcf);                 % 抓取图窗
    writeVideo(myVideo,Fr);           % 写入文件
end
close(myVideo)                        % close文件，结束