| Re: qsimp, qromb, qromo [message #39621] |
Thu, 03 June 2004 14:17  |
mmeron
Messages: 44
Registered: October 2003
|
Member |
|
|
In article <b3eff7b0.0406031146.1be4bb29@posting.google.com>, gnarloo@libero.it (michele) writes:
> Dear all,
> I'm new to this newsgroup so first of all, hello everybody!
>
> The problem i encountered is the following:
>
> I tried to use
> qsimp and qromb to integrate a couple of functions (scattering
> problems, nothing special) but the system crashes..
>
> I tried to integrate cos(x) between 0 and pi (the result is obviously
> zero)
> and these functions give out some kind of error message but not the
> result..
>
> so i turned to write the code myself as i used to do when i was an
> undrgraduate student using old riemann sums and old fortran, but as
> this approach is a bit time-consuming i have one question:
>
> is there any trick I can use or these functions only sometimes work?
> (for example if you integrate sin(x) on the same interval as above
> they give out
> the right result.
>
> Personally I can't believe the problem is in the algorithm (Simpson
> rule works, Romberg method I don't know it)...
>
I'm not much familiar with IDL's implementation of Romberg's method
since I'm using my own. So, from my own experience, I can attest that
Romberg, when properly implemented, works beautifully.
Mati Meron | "When you argue with a fool,
meron@cars.uchicago.edu | chances are he is doing just the same"
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
