Statistical: Plot Partition function & other thermodynamic variables (U, dU, Cv, F, S) - Part-2

In this video, Partition function for the systems with finite energy levels and finite no. of non interacting particles is calculated and plotted using Bose Einstein and Fermi-Dirac distribution. Using partition function, other thermodynamic variables like internal energy, energy fluctuation, Helmholtz free energy, specific heat and entropy are calculated and plotted with temperature for different no. of particles in SCILAB software.

[@flashbang8253]

Hello sir !
Is the formula correct for Bose-Einstein statistics ?
for j = 1:length(T)
b = -e/(Kb*T(j));
lz = 0
for i=1:nl
lz = lz + log(1-exp(b*((i-1)*dE-u)))
end
logz(j) = lz
z(j) = exp(lz)
end

[@arkachaudhuri4028]

Can u please write the same program in Python or kindly tell how to write it in Python.

[@spdas5942]

Sir, more examples on the three partition functions ,please .

[@pranjalibhattacharjee4321]

Sir for the Grand Cannonical partition function of BE statistics, can you tell me what values should I take? I've been trying using various values but nothings seems to work out

[@practicalhope5408]

Just change the formula for partition function. If still u r not getting tell me what values do you need

[@pranjalibhattacharjee4321]

@@practicalhope5408 I had changed the partition function. But the desired graph didn't come

[@practicalhope5408]

@@pranjalibhattacharjee4321 u can send your program and curves at [e-mail]
I will check whenever I am free

[@pranjalibhattacharjee4321]

@@practicalhope5408 okay!

[@gauri5601]

SIR PLEASE HELP ME. Write a Scilab code to draw a two-dimensional random walk from the dice rolls...plz make a coding of this question..i too much work on it but unfortunetly ,i always found error in it plzz sir help me

[@practicalhope5408]

//positive random walk depending upon the value of dice
clc;clear;clf;
x(1)=0;y(1)=0;s=1:100;
for i=1:100
x(i+1)=x(i)+round(5*rand()+1);
y(i+1)=y(i)+round(5*rand()+1);
d(i)=sqrt(x(i+1)^2+y(i+1)^2);
end
subplot(1,2,1)
plot(x,y);
subplot(1,2,2)
plot(s,d');

[@gauri5601]

@@practicalhope5408 thanks sir 🙏