2Â¥: Originally posted by xiuyouxu at 2012-04-13 17:23:46:
½¨Òé²ÉÓÃÅ£¶Ùµü´ú·¨£¬²»¹ý»ý·Ö²¿·ÖµÄ¼ÆË㻹ÊDZȽÏÂé·³£¬Ã»°ì·¨£¬Ö»ÄÜÓÃÊýÖµ»ý·Ö·½·¨Çó½â¡£

LoopConstant Q=[1000,5000,10000,50000];
Constant c1=3.741832E-16,c2=1.438786E-2, sigma=5.6688e-008,
         lambda1=2.5*1e-6;
Function Q=sigma*T^4-int(c1*x^(-5)/(exp(c2/(x*T))-1),x=lambda1,1E8);

c1:=3.741832E-16:
>  c2:=1.438786E-2:
>  sigma:=5.6688e-008:            ##˹ÌØÅ˳£Êý¦Ò
>  lambda1:=2.5*1e-6:              #»ý·ÖÏÂÏÞ ¦Ë1
> eqn:='Q=sigma*T^4-int(c1*lambda^(-5)/(exp(c2/(lambda*T)-1)),lambda=lambda1..infinity)';

Q:=1000:
> sol:=T=fsolve(Q=sigma*T^4-int(c1*lambda^(-5)/(exp(c2/(lambda*T)-1)),lambda=lambda1..infinity),T);
> eval(eqn,sol);
> eval(eqn,T= -364.440988765583);