Notre objectif principal est de modéliser numériquement la répartition des bosons et fermions au sein d'un puits de potentiel. Pour le mener à bien, il est nécessaire d'avoir accès à la probabilité d'une certaine configuration des particules dans le piège harmonique. Qui dit probabilités, dit fonctions d'ondes, puisque la probabilité d'une configuration est donnée par le module au carré de la fonction d'onde correspondante. Il devient donc indispensable de les connaître selon le type et le nombre des particules étudiées.
Pour deux fermions indépendants, la fonction d'onde s'écrit :
respectant l'antisymétrie et le principe d'exlusion de Pauli imposés par sa famille. C'est une fonction à deux dimensions que l'on peut tracer :
La probabilité en lien avec une certaine configuration s'écrit :
On peut tracer cette probabilité dans ce cas :
On la généralise pour N fermions indépendants de la manière suivante :
En ce qui concerne deux bosons sans interaction, la fonction d'onde s'exprime comme :
Ici, la symétrie bosonique est bien respectée. De la même manière, N bosons indépendants auront pour fonction d'onde :