Path: news.mathworks.com!newsfeed-00.mathworks.com!oleane.net!oleane!news.ecp.fr!aioe.org!not-for-mail
From: dpb <none@non.net>
Newsgroups: comp.soft-sys.math.maple,sci.math.num-analysis,comp.soft-sys.matlab
Subject: Re: How to take 200th order derivative, and then evaluate it numerically?
Date: Wed, 18 Mar 2009 15:01:22 -0500
Organization: Aioe.org NNTP Server
Lines: 19
Message-ID: <gprk0t$g8$1@aioe.org>
References: <f6ebcf2b-4782-4984-aebc-7cdb5b28ccee@p2g2000prf.googlegroups.com> <rbisrael.20090318173245$27f5@news.acm.uiuc.edu>
NNTP-Posting-Host: ClqEiWWk0dzWRGf3aneEDw.user.aioe.org
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
X-Complaints-To: abuse@aioe.org
NNTP-Posting-Date: Wed, 18 Mar 2009 20:02:38 +0000 (UTC)
X-Notice: Filtered by postfilter v. 0.7.7
Cancel-Lock: sha1:DIOUnZAv4nMU0HujA9f85SuPYqM=
User-Agent: Thunderbird 2.0.0.17 (Windows/20080914)
Xref: news.mathworks.com comp.soft-sys.math.maple:24556 sci.math.num-analysis:107040 comp.soft-sys.matlab:525971


Robert Israel wrote:
...
> Perhaps what you might try is computing a Taylor series up to order x^200.  
> There still might be some numerical difficulties, ... [snip for effect]

Think???  :)

>  For example, in Maple:
> 
>   f:= x^3*sin(x)/(x - cos(x));
>   Digits:= 20;
>   200! * coeff( series(f, x=0.3, 201), x-0.3, 200);
> 
>     -.90246797086526191165e446

Hmmmm....any idea how many of those digits are actually meaningful in 
front of that 10**446 exponent???

--