Procedimiento numérico
Número de moléculas N1 que tienen una energía menor que Ec
donde xc=Ec/kT
erf(x) se denomina integral de los errores
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 proporcion(double Ec){ double xc=Ec*11604.49/T; double y=2*(Math.sqrt(Math.PI)*(1.0-erfcc(Math.sqrt(xc)))/2- |