[Date Index] [Thread Index] [Author Index]

Re: Bump Function

• To: mathgroup at smc.vnet.net
• Subject: [mg31994] Re: [mg31967] Bump Function
• From: Daniel Lichtblau <danl at wolfram.com>
• Date: Sun, 16 Dec 2001 03:44:29 -0500 (EST)
• References: <200112140921.EAA03644@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

Ajitkumar wrote:
> 
> Hi,
> 
> Could  somebody tell me how to plot the graph of bump functions. For
> example, if I want to draw the graph of the function say,
> 
> f(x,y)  = x^2+y^2 for  x^2+y^2<1
>         = 2  for x^2+y^2 >2
> 
> Thanks
> Ajit Kumar
> Dept of Maths
> University of Mumbai
> Vidyanagari, Kalina
> Mumbai 400 098
> India
> Phone No. +91 652 6683
> Email- ajitkumar at math.mu.ac.in
>        ajit_kumara at hotmail.com
>  http://math.mu.ac.in/faculty/scholars/ajit.html


Below is, I think, a standard sort of "bump" function wherein we take
the value 1 for 0<=r<=1, 0 for r>=2, and smoothly move from the plateau
to 0 when 1<r<2. Here I use 'r' to indicate radial distance. We start
with a helper function that rises smoothly from 0 and falls back
smoothly to zero.

b[r_] := Exp[-1/(r*(1-r))]
normalization = NIntegrate[b[s], {s,0,1}];

By integration we use the helper to fall smoothly from 1 to 0 as r goes
from 1 to 2.

bigB[r_] /; r>=2 := 0
bigB[r_] /; r<=1 := 1
bigB[r_] /; 1<r<2 := NIntegrate[b[s], {s,r-1,1}] / normalization

Now define the bump function in {x,y} using radial distance
Sqrt[x^2+y^2].

bump[x_,y_] := bigB[Sqrt[x^2+y^2]]

You can plot it as below.

Plot3D[bump[x,y], {x,-3,3}, {y,-3,3}, PlotPoints->50]


Daniel Lichtblau
Wolfram Research

• References:
  • Bump Function
    • From: Ajitkumar <ajitkumar@math.mu.ac.in>