>

restart:

>	with(student):with(plots):

Warning, the name changecoords has been redefined

>	S1:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du1^3)));

>	S2:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du2^3)))

>

;

>	S3:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du3^3)))

>

;

>	S4:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du4^3)));

>	m_NaCl:=1e-15:

>	epsilon1:=0.001:

>	epsilon2:=0.01:

>	epsilon3:=0.1:

>	epsilon4:=1:

>	m1_SiO2:=m_NaCl/epsilon1:

>	m2_SiO2:=m_NaCl/epsilon2:

>	m3_SiO2:=m_NaCl/epsilon3:

>	m4_SiO2:=m_NaCl/epsilon4:

>	sigma:=72:

>	Mw:=0.018:

>	ns:=(2*m_NaCl)/0.58443:

>	rho:=1000000:

>	R:=8.3144;

>

T:=273;

>	du1:=(((6*m1_SiO2)/(Pi*2.65))^(1/3))/100:

>	du2:=(((6*m2_SiO2)/(Pi*2.65))^(1/3))/100:

>	du3:=(((6*m3_SiO2)/(Pi*2.65))^(1/3))/100:

>	du4:=(((6*m4_SiO2)/(Pi*2.65))^(1/3))/100:

>	p1:=semilogplot(S1-1,dp=0.916e-6..8e-6,color=blue):

>	p2:=semilogplot(S2-1,dp=4.49e-7..8e-6,color=yellow):

>	p3:=semilogplot(S3-1,dp=2.7e-7..8e-6,color=green):

>	p4:=semilogplot(S4-1,dp=2.3e-7..8e-6):

>	display(p1,p2,p3,p4);

>

restart:

>

>	with(student):with(plots):

Warning, the name changecoords has been redefined

>	S1:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du1^3))):

>	S2:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du2^3))):

>	S3:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du3^3))):

>	S4:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du4^3))):

>	m_NaCl:=1e-17:

>	epsilon1:=0.001:

>	epsilon2:=0.01:

>	epsilon3:=0.1:

>	epsilon4:=1:

>	m1_SiO2:=m_NaCl/epsilon1:

>	m2_SiO2:=m_NaCl/epsilon2:

>	m3_SiO2:=m_NaCl/epsilon3:

>	m4_SiO2:=m_NaCl/epsilon4:

>	sigma:=72:

>	Mw:=0.018:

>	ns:=(2*m_NaCl)/0.58443:

>	rho:=1000000:

>	R:=8.3144:

>

T:=273:

>	du1:=(((6*m1_SiO2)/(Pi*2.65))^(1/3))/100:

>	du2:=(((6*m2_SiO2)/(Pi*2.65))^(1/3))/100:

>	du3:=(((6*m3_SiO2)/(Pi*2.65))^(1/3))/100:

>	du4:=(((6*m4_SiO2)/(Pi*2.65))^(1/3))/100:

>	p1:=semilogplot(S1,dp=1.94e-7..2e-6,color=blue):

>	p2:=semilogplot(S2,dp=9.15e-8..2e-6,color=yellow):

>	p3:=semilogplot(S3,dp=4.55e-8..2e-6,color=green):

>	p4:=semilogplot(S4,dp=2.75e-8..2e-6):

>	display(p1,p2,p3,p4);

>

restart:

>	with(student):with(plots):

Warning, the name changecoords has been redefined

>	S1:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du1^3))):

>	S2:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du2^3))):

>	S3:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du3^3))):

>	S4:=exp(4*sigma*Mw/(rho*R*T*dp)-6*ns*Mw/(rho*Pi*(dp^3-du4^3))):

>	m_NaCl:=1e-19:

>	epsilon1:=0.001:

>	epsilon2:=0.01:

>	epsilon3:=0.1:

>	epsilon4:=1:

>	m1_SiO2:=m_NaCl/epsilon1:

>	m2_SiO2:=m_NaCl/epsilon2:

>	m3_SiO2:=m_NaCl/epsilon3:

>	m4_SiO2:=m_NaCl/epsilon4:

>	sigma:=72:

>	Mw:=0.018:

>	ns:=(2*m_NaCl)/0.58443:

>	rho:=1000000:

>	R:=8.3144:

>

T:=273:

>	du1:=(((6*m1_SiO2)/(Pi*2.65))^(1/3))/100:

>	du2:=(((6*m2_SiO2)/(Pi*2.65))^(1/3))/100:

>	du3:=(((6*m3_SiO2)/(Pi*2.65))^(1/3))/100:

>	du4:=(((6*m4_SiO2)/(Pi*2.65))^(1/3))/100:

>	p1:=semilogplot(S1,dp=4.166e-8..8e-7,color=blue):

>	p2:=semilogplot(S2,dp=1.94e-8..8e-7,color=yellow):

>	p3:=semilogplot(S3,dp=0.915e-8..8e-7,color=green):

>	p4:=semilogplot(S4,dp=4.55e-9..8e-7):

>	display(p1,p2,p3,p4);

>