initial commit / add all books

37 files changed, 352 insertions, 0 deletions

diff --git a/821/CH4/EX4.1/4_1.sce b/821/CH4/EX4.1/4_1.sce
new file mode 100755
index 000000000..841c0521c
--- /dev/null
+++ b/821/CH4/EX4.1/4_1.sce
@@ -0,0 +1,8 @@
+P=760;//pressure of 14g of nitrogen in mm of Hg//

+V=22.4;//Volume occupied by 14g of Nitrogen in litres//

+P1=380;//changed pressure of 14g of nitrogen in mm of Hg//

+V1=(P*V)/P1;//changed volume of 14g of nitrogen at 380mm pressure in litres//

+printf('Volume of 14g of Nitrogen at 380 mm Hg pressure=V1=%flitres',V1);

+V2=5.6;//changed volume of 14g of nitrogen in litres//

+P2=(P*V)/V2;//Pressure of 14g of nitrogen of volume 5.6litres//

+printf('\nPressure of 14g of Nitrogen of volume 5.6litres=P2=%fmm Hg',P2);

diff --git a/821/CH4/EX4.10/4_10.sce b/821/CH4/EX4.10/4_10.sce
new file mode 100755
index 000000000..d5ad3eadc
--- /dev/null
+++ b/821/CH4/EX4.10/4_10.sce
@@ -0,0 +1,7 @@
+printf('Density is proportional to molecular weight');

+printf('\nRelative density of O2 with respect to H2 is 32/2=16.');

+t1=12;//time for certain volume of gas to stream through hole in mins//

+p2=16;//Relative density of O2 w.r.t H2//

+t2=15;//time for same volume(as gas) of oxygen to stream through hole in mins//

+p1=((t1/t2)^2)*p2;//according to graham's law//

+printf('\nRealtive density of gas to oxygen=p1=%f',p1);//here the answer given in textbook is wrongly printed the actual answer is the one we got here through execution// 

diff --git a/821/CH4/EX4.11/4_11.sce b/821/CH4/EX4.11/4_11.sce
new file mode 100755
index 000000000..61973a02f
--- /dev/null
+++ b/821/CH4/EX4.11/4_11.sce
@@ -0,0 +1,10 @@
+T=320;//47C in kelvin//

+R=8.31*10^7;//Universal gas constant in erg per degree per mole//

+M=32;//molecular weight of O2 in gram per mole//

+C2=(3*R*T)/M;//mean square velocity of Oxygen in (cm/sec)^2//

+Crms=sqrt(C2);//Root mean square velocity of Oxygen in cm/sec//

+printf('Root mean square velocity of Oxygen=Crms=%fcm/sec',Crms);

+Cm=sqrt(8*R*T/(%pi*M));//mean velocity of Oxygen in cm/sec//

+printf('\nMean velocity of Oxygen=Cm=%fcm/sec',Cm);

+Cmpv=sqrt(2*R*T/M);//mean probable velocity in cm/sec//

+printf('\nMean probable velocity of Oxygen=Cmpv=%fcm/sec',Cmpv);//in textbook in Cmpv value it is misprinted as 10^9 but it is actually 10^4//

diff --git a/821/CH4/EX4.12/4_12.sce b/821/CH4/EX4.12/4_12.sce
new file mode 100755
index 000000000..1d2bc1061
--- /dev/null
+++ b/821/CH4/EX4.12/4_12.sce
@@ -0,0 +1,9 @@
+P=10^6;//pressure of gas in dyn per cm^2//

+p=0.00333;//density of gas in gram per cm^3//

+C2=3*P/p;//mean square velocity of gas in (cm/sec)^2//

+Crms=sqrt(C2);//Root mean square velocity of Oxygen in cm/sec//

+printf('Root mean square velocity of Gas=Crms=%fcm/sec',Crms);

+Cm=Crms/1.085;//mean velocity of Gas in cm/sec//

+printf('\nMean velocity of Oxygen=Cm=%fcm/sec',Cm);

+Cmpv=Cm/1.128;//mean probable velocity in cm/sec//

+printf('\nMean probable velocity of Oxygen=Cmpv=%fcm/sec',Cmpv);

diff --git a/821/CH4/EX4.13/4_13.sce b/821/CH4/EX4.13/4_13.sce
new file mode 100755
index 000000000..fd4cad0a4
--- /dev/null
+++ b/821/CH4/EX4.13/4_13.sce
@@ -0,0 +1,5 @@
+printf('Mean free path is inversely proportional to N which is proportional to pressure');

+L1=5.8*10^-6;//mean free path of N2 at NTP in cm//

+K=58;//preesure is raised 58 times so factor is 58//

+L2=L1/K;//mean free path of N2 at 58atm pressure in cm//

+printf('\nMean free path of N2 at 58atm pressure=1*10^-7cm');

diff --git a/821/CH4/EX4.14/4_14.sce b/821/CH4/EX4.14/4_14.sce
new file mode 100755
index 000000000..56085b18f
--- /dev/null
+++ b/821/CH4/EX4.14/4_14.sce
@@ -0,0 +1,9 @@
+T=273;//temperature 0C in kelvin//

+R=8.31*10^7;//Universal gas constant in erg per degree per mole//

+M=28;//molecular weight of N2 in gram per mole//

+printf('Since 22.4Litres of Nitrogen gas at 0C and 1atm pressure will contain 6.023*10^23Molecules.');

+N=2.69*10^19;//no. of molecules in molecules per cm^3//

+G=3.78*10^-8;//collision diameter in cm//

+K=sqrt(%pi*R*T/M);

+Z11=2*N^2*G^2*K;//number of collisions per second of Nitrogen at 0C and 1atm//

+printf('\nNumber of molecular collisions per second of Nitrogen at 0C and 1atm=%f=10.43*10^28molecular collisions per sec per cm^3',Z11);

diff --git a/821/CH4/EX4.15/4_15.sce b/821/CH4/EX4.15/4_15.sce
new file mode 100755
index 000000000..8f8a03f63
--- /dev/null
+++ b/821/CH4/EX4.15/4_15.sce
@@ -0,0 +1,18 @@
+T=273;//temperature 0C in kelvin//

+R=8.31*10^7;//Universal gas constant in erg per degree per mole//

+M=28;//molecular weight of N2 in gram per mole//

+printf('Since 22.4Litres of Nitrogen gas at 0C and 1atm pressure will contain 6.023*10^23Molecules.');

+N=2.69*10^19;//no. of molecules in molecules per cm^3//

+Cm=sqrt(8*R*T/(%pi*M));//mean velocity of Nitrogen in cm/sec//

+printf('\nMean velocity of Nitrogen=Cm=%fcm/sec',Cm);

+V=22400;//volume of nitrogen in cm^3//

+p=M/V;//Density of nitrogen in gram per cm^3//

+printf('\nDensity of Nitrogen=p=%f=1.25*10^-3gram per cm^3',p);

+n=10.99*10^-5;//Viscosity of N2 in poise//

+L=(3*n)/(Crms*p);//mean free path of nitrogen in cm//

+printf('\nMean free path of Nitrogen=L=5.81*10^-6cm');

+G=sqrt(1/(1.414*%pi*L*N));//Collision diameter of Nitrogen in cm//

+printf('\nCollission diameter of Nitrogen=G=3.80*10^-8cm');

+K=sqrt(%pi*R*T/M);

+Z11=2*N^2*G^2*K;//number of collisions per second of Nitrogen at 0C and 1atm//

+printf('\nNumber of molecular collisions per second of Nitrogen at NTP=%f=10.52*10^28molecular collisions per sec per cm^3',Z11);

diff --git a/821/CH4/EX4.16/4_16.sce b/821/CH4/EX4.16/4_16.sce
new file mode 100755
index 000000000..6fdd04dfd
--- /dev/null
+++ b/821/CH4/EX4.16/4_16.sce
@@ -0,0 +1,8 @@
+T1=27;//initial temperature in C//

+T2=177;//final temperature in C//

+printf('NH3 being a nonlinear molecule has 3 translational,3 rotational and 6 vibrational degrees of freedom.');

+printf('\nCv=3*0.5*R+3*0.5*R+6*R=9*R');

+Cv=18;//Molar heat capacity in cal per deg per mol//

+dU=Cv*(T2-T1);//Change in internal energy of a mole in cal per mole//

+printf('\nChange in internal energy of a mole of NH3=dU=%fcal per mole',dU);

+printf('\nThe actual increase in energy may not be 2700cal per mol\nbecause at the given temperature,none or only some of the vibrational degrees of freedom may be contributing to the total energy.');//answer should come as 2700 not 2400,it was just misprinted//

diff --git a/821/CH4/EX4.17/4_17.sce b/821/CH4/EX4.17/4_17.sce
new file mode 100755
index 000000000..629b244e6
--- /dev/null
+++ b/821/CH4/EX4.17/4_17.sce
@@ -0,0 +1,8 @@
+printf('CH4 being a nonlinear molecule has 3 translational,3 rotational and 9 vibrational degrees of freedom.');

+printf('\nCv=3*0.5*R+3*0.5*R+9*R=12*R');

+Cv=24;//Molar heat capacity in cal per degree per mol//

+printf('\nActual Cv is 16cal per degree per mol');

+R=2;

+Cvr=16-R;//Real molar heat capapcity in cal per degree per mole//

+printf('\nReal molar heat capapcity of a mole of CH4=Cvr=%fcal per degree per mole',Cvr);

+printf('\nThe Real molar heat capacity at constant volume is 10cal less than the theoretical value.\nSince each vibrational degree of freedom can contribute 2cal,this means that 5 vibrational degrees\nare not contributing at the given temperature.');

diff --git a/821/CH4/EX4.18/4_18.sce b/821/CH4/EX4.18/4_18.sce
new file mode 100755
index 000000000..95fdf525e
--- /dev/null
+++ b/821/CH4/EX4.18/4_18.sce
@@ -0,0 +1,10 @@
+a=3.60;//Van der Waals constant for CO2 in L^2 atm per mol^2//

+b=4.30*10^-2;//Van der Waals constant for CO2 in litre per mol//

+N=300;//No. of moles of CO2 the cylinder contains//

+V=100;//volume of cylinder in litres//

+P=150;//Maximum pressure the cylinder can withstand in atm//

+R=0.082;

+K1=P+(N^2*a/V^2);

+K2=V-(N*b);

+T=(K1*K2)/(N*R);//The lowest temperature at which cylinder can explode in Kelvin//

+printf('The lowest temperature at which cylinder can explode=T=%fK',T);

diff --git a/821/CH4/EX4.19/4_19.sce b/821/CH4/EX4.19/4_19.sce
new file mode 100755
index 000000000..3e4441505
--- /dev/null
+++ b/821/CH4/EX4.19/4_19.sce
@@ -0,0 +1,15 @@
+tA=280;//time of flow for liquid A in seconds//

+tB=200;//time of flow for liquid B in seconds//

+pA=1;//density of liquid A in gram per cm^3//

+pB=1.1;//density of liquid B in gram per cm^3//

+h=10;//height of liquid responsible for the flow in cm//

+g=980;//gravity constant in dyns//

+V=1;//volume of liquid in ml//

+L=10;//length of the capillary in cm//

+r=0.1;//radius of the capillary in cm//

+PA=h*pA*g;//Pressure of liquid A//

+PB=h*pB*g;//Pressure of liquid B//

+nA=(%pi*PA*tA*r^4)/(8*L*V);//Viscosity of Liquid A in centipoise//

+printf('\nViscosity of Liquid A=nA=%fcentipoise',nA);

+nB=(%pi*PB*tB*r^4)/(8*L*V);//Viscosity of Liquid B in centipoise//

+printf('\nViscosity of Liquid B=nB=%fcentipoise',nB);

diff --git a/821/CH4/EX4.2/4_2.sce b/821/CH4/EX4.2/4_2.sce
new file mode 100755
index 000000000..12181dacb
--- /dev/null
+++ b/821/CH4/EX4.2/4_2.sce
@@ -0,0 +1,5 @@
+V=5.6;//Volume occupied by 8g of Oxygen at 0 C in litres//

+T=273;//Temperature at which 8g of Oxygen occupies 5.6litres in Kelvin//

+V1=11.2;//Changed volume of 8g of Oxygen in litres//

+T1=(V1*T)/V;//Temperature at which 8g of Oxygen occupies 11.6litres in kelvin// 

+printf('Temperature of 8g of Oxygen Occupying 11.6litres=T1=%fKelvin=273degrees',T1);

diff --git a/821/CH4/EX4.20/4_20.sce b/821/CH4/EX4.20/4_20.sce
new file mode 100755
index 000000000..2297adf5e
--- /dev/null
+++ b/821/CH4/EX4.20/4_20.sce
@@ -0,0 +1,6 @@
+ST=520;//surface tension of mercury in dyn per cm//

+d=13.6;//density of mercury in gram per cm^3//

+g=980;//gravity constant in dyns//

+r=0.1;//radius of the capillary tube in cm//

+h=(ST*2)/(g*r*d);//Depression observed in capillary tube in cm//

+printf('Depression observed in capillary tube=h=%fcm',h);//the answer is given wrong in textbook,it should actually be double the one given//

diff --git a/821/CH4/EX4.21/4_21.sce b/821/CH4/EX4.21/4_21.sce
new file mode 100755
index 000000000..d06773e5a
--- /dev/null
+++ b/821/CH4/EX4.21/4_21.sce
@@ -0,0 +1,6 @@
+T=293;//temperature 20C in kelvin//

+R=8.31*10^7;//Universal gas constant in erg per degree per mole//

+dr=165;//Change of surface tension in erg per cm^2 per mol litre//

+C=0.06;//concentration we are considering in M//

+t=(C*dr)/(R*T);//surface excess concentration of phenol in mol per cm^2//

+printf('surface excess concentration of phenol=t=4.08*10^-10mol per cm^2');

diff --git a/821/CH4/EX4.22/4_22.sce b/821/CH4/EX4.22/4_22.sce
new file mode 100755
index 000000000..b2553a036
--- /dev/null
+++ b/821/CH4/EX4.22/4_22.sce
@@ -0,0 +1,10 @@
+M=256;//Molecular weight of acid in grams//

+w=2.56*10^-5;//weight of palmitic acid in grams//

+N=w/M;//No. of molecules of acid//

+A=123;//Total area occupied in cm^2//

+AM=A/N;//Area per molecule in cm^2//

+printf('Area per molecule=20.4*10^-16cm^2.\nThis is the area of crosssection A for a vertical position.');

+V=w/0.81;//Volume of acid in cm^3//

+VM=V/N;//Volume of one molecule of acid in cm^3//

+L=VM/AM;//Length of molecule in cm//

+printf('\nLength of a molecule=25.7*10^-8cm');

diff --git a/821/CH4/EX4.23/4_23.sce b/821/CH4/EX4.23/4_23.sce
new file mode 100755
index 000000000..51eef4328
--- /dev/null
+++ b/821/CH4/EX4.23/4_23.sce
@@ -0,0 +1,12 @@
+T=298;//Temperature in kelvin//

+ST=0.2;//Lowered surface tension in dyn per cm//

+K=1.37*10^-16;//Boltzman's constant//

+AM=(K*T)/ST;//Area per molecule on the surface in cm^2//

+A=300;//Total area occupied in cm^2//

+N=A/AM;//no. of molecules//

+W=3*10^-7;//weight of insoluble substance in grams//

+w=W/N;//weight of one molecule in grams//

+N1=6.023*10^23;

+MW=w*N1;//Molecular weight of substance in grams//

+printf('Molecular weight of substance=MW=%f=123',MW);

+

diff --git a/821/CH4/EX4.24/4_24.sce b/821/CH4/EX4.24/4_24.sce
new file mode 100755
index 000000000..ba75ff75c
--- /dev/null
+++ b/821/CH4/EX4.24/4_24.sce
@@ -0,0 +1,8 @@
+MW=461;//molecular weight of lead iodide in  grams//

+CR1=3000;//initial count rate in count per minute//

+CR2=900;//final count rate in count per minute//

+CR=CR1-CR2;

+printf('The count for lead iodide as absorbed=CR=%fcount per minute',CR);

+printf('\nThe ratio of weights of lead iodide and radio lead iodide in solution is equal to that of the same ratio on surface');

+printf('\nWeight of lead iodide in solution=0.0014grams');

+printf('\nWeight of radio lead iodide is proportional to the count.\nWeight of lead iodide on the surface=0.0014*21/9.\nMolecular weight of lead iodide=461.\nArea of the surface per gram=4266cm^2 g^-1');

diff --git a/821/CH4/EX4.25/4_25.sce b/821/CH4/EX4.25/4_25.sce
new file mode 100755
index 000000000..9bdd8994e
--- /dev/null
+++ b/821/CH4/EX4.25/4_25.sce
@@ -0,0 +1,8 @@
+printf('From the linear plot of the langmuir isotherm the intercept Xm=1/0.19=5.26milligrams');

+Xm=5.26*10^-3;

+printf('\nThis is the weight of N2 that forms a unimolecular layer of 2g charcoal.');

+MW=28;//molecular weight of N2//

+N=Xm*6.023*10^23/MW;//No. of molecules of N2//

+TA=N*16*10^-16;//Total area in cm^2//

+A=TA/2;//Area per gram in cm^2//

+printf('\nArea of N2 per gram=%fcm^2',A);

diff --git a/821/CH4/EX4.26/4_26.sce b/821/CH4/EX4.26/4_26.sce
new file mode 100755
index 000000000..c2c61ceae
--- /dev/null
+++ b/821/CH4/EX4.26/4_26.sce
@@ -0,0 +1,8 @@
+printf('From the linear plot of the langmuir isotherm the intercept=0.35*10^-3 and slope=9.47*10^-2');

+printf('\nVolume is the inverse of summation of intercept and slope and that is 10.52cc');

+Vm=10.52;//volume in cc//

+m=Vm/22400;//No. of moles of N2//

+N=m*6.023*10^23;//No. of molecules of N2//

+TA=N*16*10^-16;//Total area in cm^2//

+A=TA/17.5;//Area per gram in cm^2//

+printf('\nArea of N2 per gram=%fcm^2',A);

diff --git a/821/CH4/EX4.28/4_28.sce b/821/CH4/EX4.28/4_28.sce
new file mode 100755
index 000000000..a7c4bff41
--- /dev/null
+++ b/821/CH4/EX4.28/4_28.sce
@@ -0,0 +1,10 @@
+d=13.56;//density of mercury in gr per cm^3//

+V=1/d;//Volume of mercury in cm^3//

+D=0.07*10^-4;//Diameter of the globule in cm//

+r=D/2;//radius of globule in cm//

+printf('One globule of mercury occpies a volume of 1.33*3.14*r^3cm^3\nSurface area of one globule=4*3.14*r^2');

+Vg=1.33*3.14*r^3;//volume of one globule in cm^3//

+y=0.0738/Vg;//No of globules in 0.0738 volume//

+printf('\nNo. of globules in 0.0738cm^3=y=%f',y);

+SA=y*4*3.14*r^2;//Surface area of y globules in cm^2//

+printf('\nSurface area of y globules=SA=%fcm^2',SA);

diff --git a/821/CH4/EX4.29/4_29.sce b/821/CH4/EX4.29/4_29.sce
new file mode 100755
index 000000000..a3e55ea3d
--- /dev/null
+++ b/821/CH4/EX4.29/4_29.sce
@@ -0,0 +1,10 @@
+T=300;//temperature in kelvin//

+R=8.31*10^7;//Universal gas constant in ergs//

+r=2*10^-5;//average radius of particles in cm//

+h=4.15*10^-3;//Vertical seperation in cm//

+p=1.21;//density of latex in gram per cm^3//

+g=980;//gravity constant in dynes//

+printf('Since the dispersion medium is water its density is p1=1.\nwhen the system equilibrated average number of colloidal particles seen in the field is halved so N0/N=2.\nthe required expression is derived based on kinetic energy of particles');

+p1=1;//density of water//

+L=log(2)*0.75*(R*T)/(%pi*g*h*(p-p1)*r^3);

+printf('\nValue of Avagadro number is L=%f=6.037*10^23molecules per mol.',L);

diff --git a/821/CH4/EX4.3/4_3.sce b/821/CH4/EX4.3/4_3.sce
new file mode 100755
index 000000000..81e455d02
--- /dev/null
+++ b/821/CH4/EX4.3/4_3.sce
@@ -0,0 +1,7 @@
+P=380;//pressure of 11g of CO2 at 273K in mm of Hg//

+T=273;//Initial temperature of 11g of CO2 in kelvin//

+V=11.2;//Volume occupied by 11g of CO2 in litres at 273K//

+P1=760;//changed pressure of 11g of CO2 at 546K in mm of Hg//

+T1=546;//Final temperature of 11g of CO2 in kelvin//

+V1=(P*V*T1)/(T*P1);//changed volume of 11g of CO2 at 760mm pressure in litres//

+printf('Volume of 11g of CO2 at 760 mm Hg pressure at 546K=V1=%flitres',V1);

diff --git a/821/CH4/EX4.30/4_30.sce b/821/CH4/EX4.30/4_30.sce
new file mode 100755
index 000000000..7a90f8287
--- /dev/null
+++ b/821/CH4/EX4.30/4_30.sce
@@ -0,0 +1,12 @@
+Vf=117.857*10^-9;//Volume falling within field of view in ml//

+AN=9.66;//Average no. of particles per count//

+N=AN/Vf;//No. of particles per count//

+printf('No .of particles per ml=N=%f',N);

+W=1.5*10^-7;//Weight of iron oxide per ml//

+w=W/N;//weight of one particle in grams//

+printf('\nWeight of one particle=w=1.83*10^-15grams');

+d=5.2;//density of iron oxide in g per cm^3//

+V=w/d;//Volume of one particle in cm^3//

+printf('\nVolume of one particle=V=3.52*10^-16cm^3');

+R=2.032*10^-5;//radius of particle in cm//

+printf('Radius of particle=R=2.032*10^-5cm\nDiameter of particle=4.064*10^-5cm');

diff --git a/821/CH4/EX4.31/4_31.sce b/821/CH4/EX4.31/4_31.sce
new file mode 100755
index 000000000..3a03ae3c8
--- /dev/null
+++ b/821/CH4/EX4.31/4_31.sce
@@ -0,0 +1,8 @@
+D=78;//dielectric constant//

+V=0.009;//viscosity of suspension in dyn sec per cm^2//

+P=3.2;//potential gradient in volt per cm//

+t=38.2;

+d=0.033;//standard distance between electrodes in cm//

+u=90000;

+EP=(4*%pi*V*u*d)/(D*t*P);//Electrokinetic potential of the catalyst in volts//

+printf('Electrokinetic potential of the catalyst=EP=%fvolts',EP);

diff --git a/821/CH4/EX4.32/4_32.sce b/821/CH4/EX4.32/4_32.sce
new file mode 100755
index 000000000..6e399c1a6
--- /dev/null
+++ b/821/CH4/EX4.32/4_32.sce
@@ -0,0 +1,12 @@
+D=80.36;//dielectric constant//

+V=0.01;//viscosity of suspension in dyn sec per cm^2//

+P=10;//potential gradient in volt per cm//

+v=3*10^-3;//observed velocity in cm per sec//

+u=90000;

+EP=(4*%pi*V*u*v)/(D*P);//Zeta potential of the catalyst in volts//

+printf('Zeta potential of the catalyst=EP=%fvolts',EP);

+printf('\nThe effective thickness of the double layer can be taken to be 1cm');

+e=9*10^-4;

+r=0.5*10^-4;

+N=4*%pi*e*r^2;//total no. of charges carried by a particle//

+printf('\nTotal no. of charges carried by a particle=N=2.83*10^-11');

diff --git a/821/CH4/EX4.33/4_33.sce b/821/CH4/EX4.33/4_33.sce
new file mode 100755
index 000000000..81689875d
--- /dev/null
+++ b/821/CH4/EX4.33/4_33.sce
@@ -0,0 +1,7 @@
+D=80.36;//dielectric constant//

+V=0.01;//viscosity of suspension in dyn sec per cm^2//

+P=100;//potential gradient in volt per cm//

+u=90000;

+EP=0.03;//Electrokinetic potential of the catalyst in volts//

+FW=(P*EP*D*3600)/(4*%pi*V*u);//amount of flow of water in cm per hour//

+printf('Amount of flow of water through diaphragm=FW=%fcm per hour',FW);

diff --git a/821/CH4/EX4.34/4_34.sce b/821/CH4/EX4.34/4_34.sce
new file mode 100755
index 000000000..4a1d7b14a
--- /dev/null
+++ b/821/CH4/EX4.34/4_34.sce
@@ -0,0 +1,8 @@
+D=80.36;//dielectric constant//

+V=0.01;//viscosity of suspension in dyn sec per cm^2//

+P=40;//potential gradient in volt per cm//

+r=0.05;//radius of capillary in cm//

+u=90000;

+ZP=0.05;//Zeta potential of the catalyst in volts//

+FW=(P*ZP*D*3600*r^2)/(4*V*u);//amount of flow of water in cc per hour//

+printf('Amount of flow of water through diaphragm=FW=%fcc per hour',FW);

diff --git a/821/CH4/EX4.35/4_35.sce b/821/CH4/EX4.35/4_35.sce
new file mode 100755
index 000000000..4a44386c8
--- /dev/null
+++ b/821/CH4/EX4.35/4_35.sce
@@ -0,0 +1,14 @@
+C1=5*10^-4;

+C2=2.5*10^-4;

+x=(C1^2)/(C2+(2*C1));//The amount of Na+ transported at equilibrium in M//

+printf('The amount of Na+ transported at equilibrium=x=2*10^-4M');

+NaR=C2+x;//Na+ on RHS//

+printf('\nNa+ on RHS=4.5*10^-4M.\nCl- on RHS=x=2*10^-4M.\nNa+ on LHS=Cl- on LHS=C1-x=3*10^-4M');

+printf('\nE is in volts and F in coulombs.for homogeneity of units either R should be expressed in coulombs or F in calories.F is 96500C or 23060cal.');

+F=23060;

+n=1;

+T=300;//temperature in kelvin//

+R=2;

+K=1.5;

+E=(R*T*log(K))/(n*F);//Emf across membrane in volts//

+printf('\nEmf across membrane=E=%fvolts',E);

diff --git a/821/CH4/EX4.36/4_36.sce b/821/CH4/EX4.36/4_36.sce
new file mode 100755
index 000000000..00bb8b323
--- /dev/null
+++ b/821/CH4/EX4.36/4_36.sce
@@ -0,0 +1,8 @@
+C1=0.04;

+C2=0.02;

+x=(C1^2)/(C2+(2*C1));//The amount of Na+ transported at equilibrium in M//

+printf('The amount of Na+ transported at equilibrium=x=%fM',x);

+NaR=C2+x;//Na+ on RHS//

+printf('\nNa+ on RHS=NaR=%fM',NaR);

+NaL=C1-x;//Na+ on LHS//

+printf('\nNa+ on LHS=NaL=%fM',NaL);

diff --git a/821/CH4/EX4.37/4_37.sce b/821/CH4/EX4.37/4_37.sce
new file mode 100755
index 000000000..3cacfab9c
--- /dev/null
+++ b/821/CH4/EX4.37/4_37.sce
@@ -0,0 +1,13 @@
+MW=58.45;//Molar weight of NaCl in grams//

+d=2.17;//density of NaCl in g/cc//

+MV=MW/d;//Molar volume in cc//

+printf('Molar Volume of NaCl=MV=%fcc',MV);

+printf('\nThis must contain 6.023*10^23 NaCl units.\nIn one unit cell of NaCl there are 8 corner Na+ ions and 6 on the face centres.\nThe number of Na+ ions in one unit cell is 8*1/8+6*1/2=4.\nThere are 12Cl- ions on the edges,and one in the centre.\nThe number of Cl- ions are 12*1/4+1=4 ');

+printf('\nThe volume of the unit cell containing 4 NaCl units=179*10^-24cc');

+a=5.63*10^-8;//unit cell length in cm//

+printf('\nUnit cell length of NaCl crystal=a=5.63*10^-8cm');

+Id=a/2;//Interionic distance in cm//

+printf('\nInterionic distance in the crystal=Id=2.815*10^-8cm');

+Q=5.9;//glancing angle in degrees//

+L=2*Id*sin(Q*%pi/180)*10^8;//wavelenth in angstrums//

+printf('\nWavelength of the Xrays used=L=%fAngstrums',L);

diff --git a/821/CH4/EX4.39/4_39.sce b/821/CH4/EX4.39/4_39.sce
new file mode 100755
index 000000000..c48971985
--- /dev/null
+++ b/821/CH4/EX4.39/4_39.sce
@@ -0,0 +1,12 @@
+N=2;//No. of atoms or unit cells in BCC structure//

+L=6.023*10^23;//Avagadro number//

+a=4.291*10^-8;//Unit cell edge length in cm//

+Na=23;//weight of Na//

+p=(N*Na)/(L*a^3);//density of Na at room temperature in g/cm^3//

+printf('Density of Na at room temperature=p=%fg/cm^3',p);

+P=p*1.0398;//density of Na at -195degrees temperature in g/cm^3//

+printf('\nDensity of Na at -195degrees temperature=P=%fg/cm^3',P);

+a1=5.35*10^-8;//unit cell edge length in cm//

+N1=(P*L*a1^3)/(Na);//No. of unitcells at -195degrees//

+printf('\nNumber of unitcells at -195degrees=N1=%f',N1);

+printf('\nAt -195degrees temperature Na have 4 unitcells which means it assumes an FCC structure.');

diff --git a/821/CH4/EX4.4/4_4.sce b/821/CH4/EX4.4/4_4.sce
new file mode 100755
index 000000000..5119d8837
--- /dev/null
+++ b/821/CH4/EX4.4/4_4.sce
@@ -0,0 +1,6 @@
+P=570/760;//pressure of 10L of H2 at 300K in Atmospheres//

+T=300;//Temperature of 10L of H2 in kelvin//

+V=10;//Volume occupied by H2 at 300K in litres//

+R=0.082;//Value of R in litre-atmospheres//

+n=(P*V)/(R*T);//Number of moles of hydrogen present//

+printf('Number of moles of Hydrogen present=n=%fmol',n);

diff --git a/821/CH4/EX4.5/4_5.sce b/821/CH4/EX4.5/4_5.sce
new file mode 100755
index 000000000..224dc4db9
--- /dev/null
+++ b/821/CH4/EX4.5/4_5.sce
@@ -0,0 +1,8 @@
+TP=746;//Total pressure of gas at 298K in mm of Hg//

+PP=24;//Partial pressure of Water vapour at 298K in mm of Hg//

+PPG=TP-PP;//partial pressure of gas at 298K in mm of Hg//

+printf('partial pressure of dry gas at 298k=%fmm of Hg',PPG);

+V=200;//Volume occupied by gas at 298K in millilitres//

+P1=760;

+V1=(PPG*V)/P1;//Volume of dry gas at a pressure of 760mm Hg//

+printf('\nVolume of dry gas at pressure 760mm of Hg=V1=%fml',V1);

diff --git a/821/CH4/EX4.6/4_6.sce b/821/CH4/EX4.6/4_6.sce
new file mode 100755
index 000000000..0a48c2694
--- /dev/null
+++ b/821/CH4/EX4.6/4_6.sce
@@ -0,0 +1,16 @@
+mA=0.11;//amount of gas A in grams//

+MWA=44;//Molecular weight of gas A in grams//

+nA=mA/MWA;//No. of moles of gas A//

+printf('No. of moles of gas A=nA=%fmol',nA);

+mB=0.17;//amount of gas B in grams//

+MWB=34;//Molecular weight of gas B in grams//

+nB=mB/MWB;//No. of moles of gas B//

+printf('\nNo. of moles of gas B=nB=%fmol',nB);

+MFA=nA/(nA+nB);//mole fraction of gas A//

+printf('\nmole fraction of gas A=MFA=%f',MFA);

+MFB=nB/(nA+nB);//mole fraction of gas B//

+printf('\nmole fraction of gas B=MFB=%f',MFB);

+PPA=MFA*759;//Partial pressure of gas A//

+printf('\nPartial pressure of gas A=PPA=%fmm of Hg',PPA);

+PPB=MFB*759;//Partial pressure of gas B//

+printf('\nPartial pressure of gas B=PPB=%fmm of Hg',PPB);

diff --git a/821/CH4/EX4.7/4_7.sce b/821/CH4/EX4.7/4_7.sce
new file mode 100755
index 000000000..c6d93deef
--- /dev/null
+++ b/821/CH4/EX4.7/4_7.sce
@@ -0,0 +1,11 @@
+P1=0.3;//pressure of gas A in atm//

+VA=50;//volume of gas A in ml//

+V=500;//Volume of vessel in ml//

+PA=P1*VA/V;//Pressure of gas A in vessel//

+printf('Pressure of gas A in vessel=PA=%fatm',PA);

+P2=0.4;//pressure of gas B in atm//

+VB=250;//volume of gas B in ml//

+PB=P2*VB/V;//Pressure of gas A in vessel//

+printf('\nPressure of gas B in vessel=PB=%fatm',PB);

+P=PA+PB;//Total pressure in vessel in atm//

+printf('\nTotal pressure in vessel=P=%fatm',P);

diff --git a/821/CH4/EX4.8/4_8.sce b/821/CH4/EX4.8/4_8.sce
new file mode 100755
index 000000000..20893de3c
--- /dev/null
+++ b/821/CH4/EX4.8/4_8.sce
@@ -0,0 +1,10 @@
+printf('The volume of 1 mole of H2 at 1 atm and 273.15K=22.415Litres');

+printf('\nThe volume at 0.5atm is obtained by considering the constancy of PV');

+P1=1;

+V1=22.415;

+P2=0.5;

+V2=(P1*V1)/P2;//Volume at 0.5atm in litres//

+printf('\nVolume at 0.5atm of H2 at 273.15K=V2=%fLitres',V2);

+MH2=2;//molecular weight of H2//

+DH2=MH2/V2;//density of H2 in gram per litre//

+printf('\nDensity of H2=DH2=%fgram per litre',DH2);

diff --git a/821/CH4/EX4.9/4_9.sce b/821/CH4/EX4.9/4_9.sce
new file mode 100755
index 000000000..e17a41d91
--- /dev/null
+++ b/821/CH4/EX4.9/4_9.sce
@@ -0,0 +1,10 @@
+printf('One mole of H2 Occupies 22.415Litres at 273.15K');

+printf('\nThe volume at 300K is obtained by considering the constancy of V/T');

+T1=273.15;

+V1=22.415;

+T2=300.15;

+V2=(V1/T1)*T2;//Volume at 300K in litres//

+printf('\nVolume of 1mole of H2 at 300.15K=V2=%fLitres',V2);

+MH2=2;//molecular weight of H2//

+DH2=MH2/V2;//density of H2 in gram per litre//

+printf('\nDensity of H2=DH2=%fgram per litre',DH2);