function msd=msd(t,D,H,g,X,Y)



msd=(1/2).*D.*exp(1).^((-1/2).*D.^(-1).*t.^(-1).*(g.*H+X).^2).*t.*(( ...
  2.*pi).^(1/2).*(D.*t).^(1/2).*(g.*H+((-1)+exp(1).^((1/2).*D.^(-1) ...
  .*g.*H.*t.^(-1).*(g.*H+2.*X))).*X)+exp(1).^((1/2).*D.^(-1).*t.^( ...
  -1).*(g.*H+X).^2).*pi.*(D.*t+X.^2).*(erf(2.^(-1/2).*(D.*t).^(-1/2) ...
  .*X)+(-1).*erf(2.^(-1/2).*(D.*t).^(-1/2).*(g.*H+X)))).*(erf(2.^( ...
  -1/2).*(D.*t).^(-1/2).*Y)+(-1).*erf(2.^(-1/2).*(D.*t).^(-1/2).*(H+ ...
  Y)));