R=input('Please input R=');
whoM$ & p=0.1;
yUUg8xbpxF n=chaxun(R);
r2t|,%%N7 L=15;
oR``Jiob| c=3*10^8;
oSAO0h>0N h=6.626196*10^(-34);
NA=I7I@ k=1.3806505*10^(-23);
VXfp=JE B=198.95;
`]u!4pP" r=0.509*10^(-60);
O6
:GE'S d=532*10^(-9);
%$'YP for s=0:20
5*~]=(BE for T=0.1:50:300.1
k'(d$;Jgr if mod(s,2)==0
WUN|,P`b 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));
XA1gV>SJ y=diff('f','T');
aAT!$0H elseif mod(s,2)==1
`l%)0)T 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));
Q;h6F{i y=diff('f','T');
c] '-:= end
:gwM$2vv end
SJ
ay end
)qq5WShMJ plot3(y,s,T);
(4GDh% 怎么錯了,希望哪位大俠指點。
NKy