R=input('Please input R=');
z'qVEHc) p=0.1;
=O;SXzgE n=chaxun(R);
\OU+Kl< L=15;
N6J$z\
P c=3*10^8;
4]B3C\
v h=6.626196*10^(-34);
5pok%g
k=1.3806505*10^(-23);
*b)b#p B=198.95;
R82Y&s; r=0.509*10^(-60);
pt4xUu{ d=532*10^(-9);
*cf"l for s=0:20
vfv5ex( for T=0.1:50:300.1
r6$=|Yto if mod(s,2)==0
%7d"()L 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));
20moX7L y=diff('f','T');
ob00(?;H elseif mod(s,2)==1
*n*y!z 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));
mN