The  partition function is geven as'

Zvib(Beta)"=integrate(exp^Beta*alpha^2*h^2/4*units*m *(B-2A)+Beta*alpha^2*h^2*B^/8*units*m*(n+c)^2+Beta*alpha^2*h^2*rho^2*drho,rho=(c+n) from to lambda+c

how can plot a graph of V(r) against E0,B0 and EO*B0;

given that V(r)=V0/b(a-exp**-eta*r/1-exp**-eta*r)-V1exp^-eta*r/(1-exp^-eta*r)+V2exp^-2eta*r/(1-exp^-eta*r)^2;

E0=q*rho(electron)*velocity(vs)/mu;

B0=h/2^2*omega^2/16