Формула диффузии одномерного случая

F=[2.5]
q=[1]
n=[2.14]
m=[0.81]
A=[200]
M=[47.27]
H=[8]
T=[14]
V1=[1.2]
Cm=(A*M*F*m*n*q)./H.^2*(V1*T).^1/3
D=[1.2]; 
t=[1:24];
x1=[100];
x2=[50]
y1=sqrt(1./(4*pi*D*t.*Cm))*exp(-(x1-25).^2./(4*D*t));
y2=sqrt(1./(4*pi*D*t.*Cm))*exp(-(x2-1.)^2./(4*D*t));
plot(x,y1);
hold on;
plot(x,y2);
title('Формула диффузии одномерного случая');
xlabel('x'); ylabel('c');

F = 2.5; q = 1;
n = 2.14; m = 0.81;
A = 200; M = 47.27;
H = 8; T = 14;
V1 = 1.2;
Cm = (A*M*F*m*n*q)/H^2*(V1*T)^1/3;
D = 1.2; t = 1:24;
x1 = 100; x2 = 50;
y1 = sqrt(1./(4*pi*D.*t*Cm)).*exp(-(x1-25)^2./(4*D.*t));
y2 = sqrt(1./(4*pi*D.*t*Cm)).*exp(-(x2-1)^2./(4*D.*t));
plot(t,y1);
hold on;
plot(t,y2);
title('Формула диффузии одномерного случая');
xlabel('t'); ylabel('c(t)');

clear ,clc
 
F = 2.5; q = 1;
n = 2.14; m = 0.81;
A = 200; M = 47.27;
H = 8; T = 14;
V1 = 1.2;
Cm = (A*M*F*m*n*q)/(H^2*(V1*T)^(1/3));
D = 1.2; 
x = 0:0.1:20;
t1 = 2; t2 = 2.2;
x01 = 6; x02 = 12;
y1 = sqrt(1./(4*pi*D*t1*Cm)).*exp(-(x-x01).^2./(4*D*t1));
y2 = sqrt(1./(4*pi*D*t2*Cm)).*exp(-(x-x02).^2./(4*D*t2));
plot(x,y1,'b',x,y2,'r');
str = {['t=',num2str(t1)],['t=',num2str(t2)]};
legend(str)
title('Формула диффузии одномерного случая');
xlabel('x'); ylabel('c(x)');