Help with a loglogplot in maple please

  • Context: Maple 
  • Thread starter Thread starter elarson89
  • Start date Start date
  • Tags Tags
    Maple
Click For Summary
SUMMARY

The discussion focuses on creating a log-log plot of Planck's Law using Maple software. The user encountered an empty plot error due to incorrect constants, specifically the value of Boltzmann's constant. The correct value is k:=1.381*10^(-23), which resolves the issue by ensuring the exponential function behaves correctly within the specified domain. The user successfully corrected the code and received assistance from the community.

PREREQUISITES
  • Familiarity with Maple software for mathematical plotting
  • Understanding of Planck's Law in physics
  • Knowledge of logarithmic functions and their applications in plotting
  • Basic grasp of physical constants such as Boltzmann's constant
NEXT STEPS
  • Explore advanced plotting techniques in Maple, including customization options
  • Study the implications of Planck's Law in thermal radiation
  • Learn about the behavior of exponential functions in mathematical modeling
  • Investigate common errors in mathematical plotting and how to troubleshoot them
USEFUL FOR

Students in physics or engineering, educators teaching thermodynamics, and anyone interested in using Maple for scientific plotting and analysis.

elarson89
Messages
19
Reaction score
0
Help with a loglogplot in maple please!

Hi all,

I have some physics homework that I'm trying to finish tonight. I have it all done and everything, and now all i have to do is make a plot. I'm trying to do a loglogplot of Planck's Law. Here's my code:


with(plots);

c:=2.998*10^8;
k:=1.381*10^(23);
h:=6.626*10^(-34);
t:=4500;

Plancks_Law:=x->(8*Pi*h/(c^3))*(x^3/(exp(h*x/(k*t))-1));



and all I'm doing to try to plot is, loglogplot(Plancks_Law(x),x=1..10000);

I don't care what the domain is, but anything i try gives me an empty plot error. Any suggestions??

Thanks
 
Physics news on Phys.org


I think you have some constants wrong.
When I evaluate h / (k t) with your numbers, I get something like 10-60 while it should be like 10-14.
As a consequence, unless you take x = 1060, your exponential will be very close to 1, so your denominator becomes 1 - 1 = 0.

Follow-up: try k:=1.381*10^(-23) instead of k:=1.381*10^(23).
 


genius! I had boltzmann's constant to be of order 10^23, not 10^-23. Thanks a ton!
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 0 ·
Replies
0
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 9 ·
Replies
9
Views
5K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
9K