Re: [Newbie] Interpreting output
- To: mathgroup at smc.vnet.net
- Subject: [mg56651] Re: [mg56642] [Newbie] Interpreting output
- From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
- Date: Sun, 1 May 2005 03:13:33 -0400 (EDT)
- References: <200505010446.AAA13098@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
On 1 May 2005, at 13:46, MC wrote:
> Hi everybody....
>
> I'm having troubles in interpreting the output of Mathematica.
>
> My problem is to define a function F(x,y,z) such that:
> F = (whatever) if x>c
> F=(whatever) if x<=c
>
> I achievied this using /;
>
> When I derive F with respect to one of the variables, the result is a
> mess.
> Just an example with 1 variable:
>
> f=x/;x>0;
> f=-x^2/;x<0;
>
> D[f,x]=-2Condition(1,0)[x^2,x<0]+Condition(1,0)[x^2,x<0]Less(1,0)[x,0]
>
> Even in this very simple case, the result is a real mess:
> I really do not know how should I read that...and I'm wondering if
> I have defined the function in the right way.
>
> Anybody so kind to explain me what's goung on?
>
> Thank you very much for your patience
>
>
What is going on is simply that you that Mathematica's notion of
derivative does not work with piecewise functions defined by means of
pattern matching. (Actually, there is nothing in the documentation that
suggests that it might.)
In order for differentiation to be possible you have to construct a
piecewise function using a construction that "understands" the notion
of a derivative. In Mathematica 5.1 this is easy:
f[x_] := Piecewise[{{x^2, x < 2}, {x^3, x >= 2}}]
f'[1]
2
f'[3]
27
f'[2]
Indeterminate
In earlier version you could try doing this by means of the UnitStep
function:
g[x_] := x^2*UnitStep[2 - x] + x^3*UnitStep[x - 2]
this will work almost the same except for the singular value
g'[1]
2
g'[3]
27
g'[1]
2
g'[3]
27
g'[2]
4*DiracDelta[0] + 16
Note also that both approaches work well when you construct a
differentiable everywhere piecewise function:
f[x_] := Piecewise[{{(x-1)^2, x < 1}, {(x-1)^3, x >= 1}}]
f'[1]
0
g[x_]:=(x-1)^2*UnitStep[1-x]+(x-1)^3*UnitStep[x-1]
g'[1]
0
However, this is achieved in a quite different way, which we can see
from
f'[x]
Piecewise[{{2*(x - 1), x < 1}, {0, x == 1}}, 3*(x - 1)^2]
and
g'[x]
DiracDelta[x - 2]*(x - 2)^3 - DiracDelta[x - 2]*
(x - 2)^2 + 3*UnitStep[x - 2]*(x - 2)^2 +
2*UnitStep[2 - x]*(x - 2)
The second answer is given as a "generalized function" and although it
seems "correct" I would not rely on this type of use of generalised
functions (distributions) in anything but the most simple types of
situations. The Piecewise approach looks more promising.
Andrzej Kozlowski
Chiba, Japan
http://www.akikoz.net/andrzej/index.html
http://www.mimuw.edu.pl/~akoz/
- References:
- [Newbie] Interpreting output
- From: "MC" <daemmerung@enttaeuschung.com>
- [Newbie] Interpreting output