comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Integration qpint

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Integration qpint [message #81061]

Sun, 12 August 2012 10:42

Gompie
Messages: 76
Registered: August 2012

Member

Hi
Can anyone help. I am using qpint1d to integrate but i get
"Program caused arithmetic error: Floating overflow:"

The code is
C1 = 1.1910427E-12
c2=1.43883
T=300
c=c1*(c2^(-4))*T^4
llimit=c2*2380/T
ulimit=c2*2940/T
print,"X=",llimit,ulimit
lin= c*qpint1d('x^3/(EXP(X)-1)', llimit, +inf, /expr)
rin=c*qpint1d('x^3/(EXP(X)-1)', ulimit, +inf, /expr)
print,abs(lin-rin)
end

This outputs a number ( i am not sure if it is correct) but it also says floating point error.
Thanks in advance
Gompie

Report message to a moderator

Previous Topic:	RE: How to add trend line over Coyote graphics plot
Next Topic:	Re: Integration qpint

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 14:00:59 PDT 2025

Total time taken to generate the page: 0.00497 seconds