cintrue.m

来自「这是在网上下的一个东东」· M 代码 · 共 36 行

%
% this routine will plot the true source
% history for the Skaggs and Kabala paper.
%
% function cin=cintrue(t)
%
% t = time array at which true solution is sought
% cin = true concentration vector
%
function cin=cintrue(t)

nt=size(t,2);

nterms=3;
off=zeros(nterms,1);
m=off;
d=off;
off(1)=130.;
off(2)=150.;
off(3)=190.;
m(1)=1.;
m(2)=0.3;
m(3)=0.5;
d(1)=5;
d(2)=10;
d(3)=7;

cin=zeros(1,nt);

for i=1:nterms
	cin=cin+m(i)*exp(-(t-off(i)).^2/(2*d(i)^2));
end % for i

return