clear all
clc
format short e
a=2.46e-10;
y=2.7*(1.602e-19);
K=1.38e-23;
T=300;
q=(1.055e-34)*((4*pi)/(3*a));
s=input('Enter a number:');
m=input('Enter a number:');
syms p
Ep=y*sqrt(1+4*cos((pi*s)/m)^2+(4*cos((pi*s)/m)*cos(a*p)));
v=diff(Ep,p);
Ams=(1/pi)*int(v*sin(m*p),p,-pi,pi); 
F0=1/(1+exp(Ep/(K*T)));
Bms=int(F0*cos(m*a*p),p,-q,q);
Cm=Ams*Bms

Ep=-y*sqrt(1+4*cos((pi*s)/m)^2+(4*cos((pi*s)/m)*cos(a*p)));
v=diff(Ep,p);
Ams=(1/pi)*int(v*sin(m*p),p,-pi,pi);
F0=1/(1+exp(Ep/(K*T)));
Bms=int(F0*cos(m*a*p),p,-q,q);
Cm=Ams*Bms