       Re: Re: DSolve and N do not commute

• To: mathgroup at smc.vnet.net
• Subject: [mg40242] Re: [mg40213] Re: [mg40184] DSolve and N do not commute
• From: Selwyn Hollis <selwynh at earthlink.net>
• Date: Thu, 27 Mar 2003 06:51:02 -0500 (EST)
• Sender: owner-wri-mathgroup at wolfram.com

```Very strange... No such "bug" in Version 4.1.5 for Mac OS X:

eq := (f''[r] + 8 2^(2/3)/3 f'[r] == -8 2^(1/3)/3 f[r])

DSolve[eq, f[r], r]

{{f[r] -> C/E^((2/3)*2^(2/3)*r) + C/E^(2*2^(2/3)*r)}}

DSolve[N[eq], f[r], r]

{{f[r] -> C/E^(3.1748*r) + C/E^(1.05827*r)}}

-----
Selwyn Hollis
http://www.math.armstrong.edu/faculty/hollis

On Wednesday, March 26, 2003, at 02:42  AM, Dr Bob wrote:

> It's a bug.  Here's the general solution:
>
> <snip>

> On Tue, 25 Mar 2003 03:04:52 -0500 (EST), arkadas ozakin
>
>> Take the following differential equation,
>>
>> eq = (f''[r] + 8 2^(2/3)/3 f'[r] == -8 2^(1/3)/3 f[r])
>>
>> DSolve[eq, f[r], r] gives something like
>>
>> f[r] -> C Exp[k1 r] + C Exp[k2 r]
>>
>> where k1 and k2 are two different numbers.
>>
>> However,
>>
>> DSolve[N[eq],f[r],r] gives something like
>>
>> f[r] -> (C + C r) Exp[k2num r]
>>
>> where k2num is the numerical version of the constant k2 above.
>>
>> I don't know how DSolve handles numerical constants, but the
>> discrepancy between the two results wasn't something I expected
>> (that's why it took me quite a bit of time to figure out what was
>> going wrong in my longish Mathematica notebook...)
>>
>> Does anyone know why this happened? Is this a bug, or am I doing
>> something wrong? Any suggestions for avoiding similar things in the
>> future?
>>
>>
>>
>
>
>
> --
> majort at cox-internet.com
> Bobby R. Treat
>
>
>

```

• Prev by Date: Re: matrix differential equations with NDSolve
• Next by Date: Re: Q: (how?) notebooks with I/O loop
• Previous by thread: Re: DSolve and N do not commute
• Next by thread: Re: Re: Re: DSolve and N do not commute