Ecuacion general del modelo de escape instantaneo. C(x, y, z, M, u, He, t, est) donde est es "A", "B", "C", "D", "E" o "F"

FuitaInst.hpprgm

#cas
FuitaInst(x, y, z, M, u, He, t, est):=
BEGIN
 LOCAL a := 0;
  LOCAL b := 0;
  LOCAL c := 0;
  LOCAL d := 0;
  CASE
    IF est == "A" THEN
      a := 0.527;
      b := 0.865;
      c := 0.28;
      d := 0.9;
    END;
    IF est == "B" THEN
      a := 0.371;
      b := 0.866;
      c := 0.23;
      d := 0.85;
    END;
    IF est == "C" THEN
      a := 0.209;
      b := 0.897;
      c := 0.22;
      d := 0.8;
    END;
    IF est == "D" THEN
      a := 0.128;
      b := 0.905;
      c := 0.2;
      d := 0.76;
    END;
    IF est == "E" THEN
      a := 0.098;
      b := 0.902;
      c := 0.15;
      d := 0.73;
    END;
    IF est == "F" THEN
      a := 0.065;
      b := 0.902;
      c := 0.12;
      d := 0.67;
    END;
  END;
  return (M/((2*π)^(3/2)*0.13*x*0.5*a*x^b*c*x^d))*EXP(-0.5*((x-u*t)/(0.13*x))^2)*EXP(-0.5*(y/(0.5*a*x^b))^2)*(EXP(-0.5*((z-He)/(c*x^d))^2)+EXP(-0.5*((z+He)/(c*x^d))^2));
END;
#end