R=input('Please input R='); -_>c P
p=0.1; asPD>j c
n=chaxun(R); d 'x;]#S
L=15; <0v'IHlZ8
c=3*10^8; RM%lhDFY
h=6.626196*10^(-34); PM%./
k=1.3806505*10^(-23); >tib21*
B=198.95; eA{,=,v)
r=0.509*10^(-60); m_\CK5T_
d=532*10^(-9); YJrK oK}
for s=0:20 //H+S
q66
for T=0.1:50:300.1 =wS:)%u
if mod(s,2)==0 Og30&a!~F
f=p*n*L*6*((2*h*c)/(9*k*T))*((3/2)*(s+1)*(s+2)/(2*s+3))*(64*(pi^2)/15)*((c/d-(4*s+6)*B)^4)*r*exp(-(h*c*B/(k*T))*s(s+1)); #z~D1Zl
y=diff('f','T'); YwB5Zqr
elseif mod(s,2)==1 .}Bb
:*@
f=p*n*L*3*((2*h*c)/(9*k*T))*((3/2)*(s+1)*(s+2)/(2*s+3))*(64*(pi^2)/15)*((c/d-(4*s+6)*B)^4)*r*exp(-(h*c*B/(k*T))*s(s+1)); K8284A8v
y=diff('f','T'); O}IS{/^7
end ~aA+L-s|
end 6kIq6rWF9
end f4!^0%l
plot3(y,s,T); *zz/U
(9D
怎么错了,希望哪位大侠指点。 2z )h,<D