From: "Stefan Koch" <s.koch@seh-zuelpich.de>
To: <gsl-discuss@sourceware.cygnus.com>
Subject: hypergeometric function
Date: Wed, 19 Dec 2001 13:20:00 -0000 [thread overview]
Message-ID: <NEBBIBGCPKOPEIODNFDIAEIECAAA.s.koch@seh-zuelpich.de> (raw)
Hello,
I'm a new user of the gsl-library. Everything works well. But when I tried
to look at the hypergeometric function, I recognized a discontinuity (jump)
in the graphic-drawing of the function. I think this should not be so. I
used the following code segment:
void Test(void)
{
double AktVal, x;
int n, p, LoopX, MidX, MidY, CalcSize, gsl_error;
CSize sizeTotal;
gsl_sf_result result;
sizeTotal = GetTotalSize();
MidX = sizeTotal.cx/2.0;
MidY = sizeTotal.cy/2.0;
n = 3;
p= 8;
CalcSize = 1000;
x = 0;
gsl_error = gsl_sf_hyperg_1F1_int_e(-n, p+1, x, &result) ;
if(gsl_error) AfxMessageBox(gsl_strerror (gsl_error));
AktVal = result.val;
pDC->MoveTo(0, MidY - AktVal);
for(LoopX = 0; LoopX < CalcSize; LoopX++ )
{
x += 0.1;
gsl_error = gsl_sf_hyperg_1F1_int_e(-n, p+1, x, &result);
if(gsl_error) AfxMessageBox(gsl_strerror (gsl_error));
AktVal = 100.0*result.val;
pDC->LineTo( LoopX, MidY - AktVal);
}
}
(Visual C++ 6.0, NT 4.0)
Can anybody tell me if there's something wrong with my code or with the
function? No error message occured.
Thanks for your help.
Best regards
Stefan Koch
mailto:s.koch@seh-zuelpich.de
WARNING: multiple messages have this Message-ID
From: "Stefan Koch" <s.koch@seh-zuelpich.de>
To: <gsl-discuss@sourceware.cygnus.com>
Subject: hypergeometric function
Date: Sat, 01 Dec 2001 23:20:00 -0000 [thread overview]
Message-ID: <NEBBIBGCPKOPEIODNFDIAEIECAAA.s.koch@seh-zuelpich.de> (raw)
Message-ID: <20011201232000.cqFZrUSiqUADVD5TSfVOMkX-wtdkMQi_JS-IPG2F7jY@z> (raw)
Hello,
I'm a new user of the gsl-library. Everything works well. But when I tried
to look at the hypergeometric function, I recognized a discontinuity (jump)
in the graphic-drawing of the function. I think this should not be so. I
used the following code segment:
void Test(void)
{
double AktVal, x;
int n, p, LoopX, MidX, MidY, CalcSize, gsl_error;
CSize sizeTotal;
gsl_sf_result result;
sizeTotal = GetTotalSize();
MidX = sizeTotal.cx/2.0;
MidY = sizeTotal.cy/2.0;
n = 3;
p= 8;
CalcSize = 1000;
x = 0;
gsl_error = gsl_sf_hyperg_1F1_int_e(-n, p+1, x, &result) ;
if(gsl_error) AfxMessageBox(gsl_strerror (gsl_error));
AktVal = result.val;
pDC->MoveTo(0, MidY - AktVal);
for(LoopX = 0; LoopX < CalcSize; LoopX++ )
{
x += 0.1;
gsl_error = gsl_sf_hyperg_1F1_int_e(-n, p+1, x, &result);
if(gsl_error) AfxMessageBox(gsl_strerror (gsl_error));
AktVal = 100.0*result.val;
pDC->LineTo( LoopX, MidY - AktVal);
}
}
(Visual C++ 6.0, NT 4.0)
Can anybody tell me if there's something wrong with my code or with the
function? No error message occured.
Thanks for your help.
Best regards
Stefan Koch
mailto:s.koch@seh-zuelpich.de
next reply other threads:[~2001-12-19 13:20 UTC|newest]
Thread overview: 21+ messages / expand[flat|nested] mbox.gz Atom feed top
2001-12-19 13:20 Stefan Koch [this message]
2001-12-01 23:20 ` Stefan Koch
2001-12-19 13:20 ` Dan, Ho-Jin
2001-12-02 18:35 ` Dan, Ho-Jin
2001-12-19 13:20 ` Dan, Ho-Jin
2001-12-19 13:20 ` Brian Gough
2001-12-03 8:21 ` Brian Gough
[not found] ` <NEBBIBGCPKOPEIODNFDIEEIHCAAA.s.koch@seh-zuelpich.de>
2001-12-19 13:20 ` AW: " Brian Gough
2001-12-03 15:27 ` Brian Gough
2001-12-19 13:20 ` J.J. Gomez Cadenas
2001-12-05 2:07 ` J.J. Gomez Cadenas
2001-12-19 13:20 ` Brian Gough
2001-12-05 12:20 ` Brian Gough
2001-12-19 13:20 ` J.J. Gomez Cadenas
2001-12-05 12:20 ` J.J. Gomez Cadenas
2001-12-19 13:20 ` Brian Gough
2001-12-06 8:34 ` Brian Gough
2001-12-19 13:20 ` J.J. Gomez Cadenas
2001-12-06 9:15 ` J.J. Gomez Cadenas
2001-12-19 13:20 ` J.J. Gomez Cadenas
2001-12-19 13:20 ` J.J. Gomez Cadenas
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=NEBBIBGCPKOPEIODNFDIAEIECAAA.s.koch@seh-zuelpich.de \
--to=s.koch@seh-zuelpich.de \
--cc=gsl-discuss@sourceware.cygnus.com \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).