Enfriamiento de un gas

Procedimiento numérico

La temperatura T₁del gas después de eliminar las moléculas más energéticas del gas, es decir, aqullas cuya energía es superior a la de corte E_c es

$T_{1} = \frac{T_{0} (3 \sqrt{π} erf (\sqrt{x_{c}}) - (4 x_{c} \sqrt{x_{c}} + 6 \sqrt{x_{c}}) \exp (- x_{c}))}{3 (\sqrt{π} erf (\sqrt{x_{c}}) - 2 \sqrt{x_{c}} \exp (- x_{c}))}$

El número de moléculas que permenecen en el recipiente es

$N_{1} = N_{0} (erf (\sqrt{x_{c}}) - \frac{2}{\sqrt{π}} \sqrt{x_{c}} \exp (- x_{c}))$

donde x_c=E_c/kT

erf(x) se denomina integral de los errores

$erf (x) = \frac{2}{\sqrt{π}} \int_{0}^{x} \exp (- y^{2}) d y$

double erfcc(double x) {
	double t,z,ans;
	z=Math.abs(x);
	t=1.0/(1.0+0.5*z);
	ans=t*Math.exp(-z*z-1.26551223+t*(1.00002368+t*(0.37409196+t*(0.09678418+
		t*(-0.18628806+t*(0.27886807+t*(-1.13520398+t*(1.48851587+
		t*(-0.82215223+t*0.17087277)))))))));
	return x >= 0.0 ? ans : 2.0-ans;
}
  double erf(double x){
    return (1-erfcc(x));
  }
   void setNuevo(double temperatura){
     T0=temperatura;
     N0=1.0;
  }
  void setElimina(double Ec){
     this.Ec=Ec;
  }
  void setEquilibrio(){
    double xc=Ec*11604.49/T0;
    double T1=T0*(3*Math.sqrt(Math.PI)*erf(Math.sqrt(xc))-(4*xc*Math.sqrt(xc)+6*Math.sqrt(xc))*Math.exp(-xc))
/(3*(Math.sqrt(Math.PI)*erf(Math.sqrt(xc))-2*Math.sqrt(xc)*Math.exp(-xc)));
    double N1=N0*(erf(Math.sqrt(xc))-2*Math.sqrt(xc)*Math.exp(-xc)/Math.sqrt(Math.PI));
    T0=T1;
    N0=N1; //para la siguiente iteracción
  }
  double funcion(double x){ //distribución
     double y=2*N0*Math.sqrt(x)*Math.exp(-11604.49*x/T0)*1250083.9/
(Math.sqrt(Math.PI)*Math.pow(T0, 1.5));
     return y;
  }