Les dejamos un programa escrito en Matlab/Octave p2Bose.m para calcular las propiedades del gas de Bose (N0, U,Cv,etc) del problema 2 de la guia en función de número de partículas y temperatura. Aquellos que sepan usar algun otro lenguaje o programa los invitamos a usarlo.
A modo de ejemplo les dejamos un gráfico producido por el programa p2Bose.m para el gas.
Follow
¿Cómo sabés que la solución de la función cuadrática para z(N,T,V) es la que tiene el discriminante restando?
Lo que pasa que la otra solucion de la cuadrática da un z que es mayor que uno. Eso lo ves facil por ejemplo en el limite de T=0 donde una solución de la cuadratica da N/(N+1) como debe ser y la otra 2, el limite de T infinito tambien se ve facil y la otra solucion de la cuadratica da un z-> 1, cuando deberia ser cero por el limite clasico de alta temperatura.
Gracias Pablo, muy útil todo. Intenté usar el mismo código para adaptarlo al problema tres, pero algunas cosas me fallan, después te voy a consultar. Una cosa que me resultó útil para leer tu código fue convertir las expresiones a simbólicas y usar el comando 'pretty' para verlas en escritura en usual. Una cosa que te quiero preguntar luego sin que se me pase es si hay una forma facil de graficar las fn(z) y las gn(z).