Path: news.mathworks.com!newsfeed-00.mathworks.com!newsfeed2.dallas1.level3.net!news.level3.com!postnews.google.com!g30g2000yqc.googlegroups.com!not-for-mail
From: "moiseev.igor" <moiseev.igor@gmail.com>
Newsgroups: comp.soft-sys.matlab
Subject: Re: MATLAB computes Incomplete Elliptic Integrals incorrectly?
Date: Tue, 27 Apr 2010 15:36:18 -0700 (PDT)
Organization: http://groups.google.com
Lines: 19
Message-ID: <3238d66a-fc90-4ed9-96c3-c2314a02c9c8@g30g2000yqc.googlegroups.com>
References: <hqujk9$7o2$1@fred.mathworks.com>
NNTP-Posting-Host: 151.65.42.177
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
X-Trace: posting.google.com 1272407778 10589 127.0.0.1 (27 Apr 2010 22:36:18 GMT)
X-Complaints-To: groups-abuse@google.com
NNTP-Posting-Date: Tue, 27 Apr 2010 22:36:18 +0000 (UTC)
Complaints-To: groups-abuse@google.com
Injection-Info: g30g2000yqc.googlegroups.com; posting-host=151.65.42.177; 
	posting-account=v1AK6woAAABv3opmjcojTn7lnklV0HSz
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/5.0 (X11; U; Linux i686; en-US) AppleWebKit/533.4 
	(KHTML, like Gecko) Chrome/5.0.375.19 Safari/533.4,gzip(gfe)
Xref: news.mathworks.com comp.soft-sys.matlab:630432

You may try to use the "elliptic12" function from the elliptic package
http://elliptic.googlecode.com/

Just get the single file elliptic12.m, put it in your PATH and it
gives the results for incomplete ellip. integrals F and E.

The only difference with Maple is that it calculates the integrals
with different to maple representation,
check the sample code for "asin" and square of "sqrt((1-.6^2)/(1-.
4^2))" to obtain the same results!

> > mfun('EllipticE',sqrt(1-.4^2),sqrt((1-.6^2)/(1-.4^2)))
> > ans =  0.9863

[F,E] = elliptic12(asin(sqrt(1-.4^2)),(1-.5^2)/(1-.4^2))
F =
   1.474960073006188
E =
   0.949580324404286