Исправить ошибки в программе построения графиков распределения

m=1;
X=-4:0.05:4;
Y=-4:0.05:4
[x,t]=meshgrid(X,[0.1 0.2 0.5 1]);
[x,t]=meshgrid(Y,[0.1 0.2 0.5 1]);
Dx=1;
Dy=0.2;
Lz=2;
a='t';
for ii=1:2
    figure
     plot(X,m./(Lz.*4.*pi.*t.*sqrt(Dx*Dy*t)).*exp(-x.^2./(4*Dx*t)-y.^2/(4*Dy*t)))
    grid on
   legend([a,'=0.1'],[a,'=0.2'],[a,'=0.5'],[a,'=1'])
    xlabel('x')
    ylabel('c(x,t)')
    D=t;
    t=0.5;
    a='D';
end

m=1;
X=-4:0.05:4;
Y=-4:0.05:4;
[x,t]=meshgrid(X,[0.1 0.2 0.5 1]);
[y,t]=meshgrid(Y,[0.1 0.2 0.5 1]);
Dx=1;
Dy=0.2;
Lz=2;
a='t';
for ii=1:2
    figure
     plot(X,m./(Lz.*4.*pi.*t.*sqrt(Dx*Dy*t)).*exp(-x.^2./(4*Dx*t)-y.^2./(4*Dy*t)))
    grid on
   legend([a,'=0.1'],[a,'=0.2'],[a,'=0.5'],[a,'=1'])
    xlabel('x')
    ylabel('c(x,t)')
    D=t;
    t=0.5;
    a='D';
end

m=1;
X=-4:0.05:4;
[x,t]=meshgrid(X,[0.1 0.2 0.5 1]);
y=0;
Dx=1;
Dy=0.2;
Lz=2;
a='c(x,0,t)';
b='x';
for ii=1:2
    figure
    plot(X,m./(Lz*4*pi*t.*sqrt(Dx*Dy)).*exp(-x.^2./(4*Dx*t)-y.^2./(4*Dy*t)))
    grid on
    legend('t=0.1','t=0.2','t=0.5','t=1')
    xlabel(b)
    ylabel(a)
    y=x;
    x=0;
    a='c(0,y,t)';
    b='y';
end