Now it works! Thank you so much!!

Now it works! Thank you so much!!

I was just trying to get some code I saw online to work, because the scale I get when I run the code is different from the output in the webpage from which I took the code, that's why I thought the scale might be the problem. I tried semilogplot, instead of plot but the code didn't work. The code in the website was written in maple V 4 and I use Maple 10, can that be the problem? Here is the picture: http://www.rose-hulman.edu/Users/groups/packets/HTML/modern/lghtbulb1.gif And here is the code: -------------------------- #Clear all variables. restart; #Load the plots package with(plots): #Planck's radiation formula in eV u:=6.2415E18*8*Pi*h*nu^3/(c^3*exp(h*nu/(k*T))-1); #define constants h:=6.626*10^(-34); k:=1.381*10^(-23); c:=2.998*10^8; #Display the function is decimal form. evalf(u); #reset variable T:='T': #Generate functions of u for different values of T for q from 2 to 7 do u.q:=subs(T=q*1000,u); od; #Plot these function on one graph. plot({u.(2..8)},nu=0..20*10^14,title=`Energy Distribution Plots for T=2000K-7000K`); ---------------------- I would appreciate it if you can tell me whether the code gives an output as what in the picture or if it's my maple prog that's missing up. Thank you in advance,

I was just trying to get some code I saw online to work, because the scale I get when I run the code is different from the output in the webpage from which I took the code, that's why I thought the scale might be the problem. I tried semilogplot, instead of plot but the code didn't work. The code in the website was written in maple V 4 and I use Maple 10, can that be the problem? Here is the picture: http://www.rose-hulman.edu/Users/groups/packets/HTML/modern/lghtbulb1.gif And here is the code: -------------------------- #Clear all variables. restart; #Load the plots package with(plots): #Planck's radiation formula in eV u:=6.2415E18*8*Pi*h*nu^3/(c^3*exp(h*nu/(k*T))-1); #define constants h:=6.626*10^(-34); k:=1.381*10^(-23); c:=2.998*10^8; #Display the function is decimal form. evalf(u); #reset variable T:='T': #Generate functions of u for different values of T for q from 2 to 7 do u.q:=subs(T=q*1000,u); od; #Plot these function on one graph. plot({u.(2..8)},nu=0..20*10^14,title=`Energy Distribution Plots for T=2000K-7000K`); ---------------------- I would appreciate it if you can tell me whether the code gives an output as what in the picture or if it's my maple prog that's missing up. Thank you in advance,