Mathematica: Plotting piecewise functions

  • Mathematica
  • Thread starter Niles
  • Start date
  • #1
1,868
0
Hi guys

Take a look at

Code:
DensityPlot[If[x == 0 && y == 0, 1, 0], {x, -1, 1}, {y, -1, 1}]

This just gives me a uniform plot. Why does the peak at the origin not show up?

Sincerely,
Niles.
 

Answers and Replies

  • #2
422
1
There are two main reasons why you'll have trouble with attempting to plot functions like that using DensityPlot.

  • The function is zero everywhere but at a single point. DensityPlot[] is unsuitable for functions like this, as a few minutes thought should convince you.
  • Even for other function with jump discontinuities, the number of plot points and the maximum recursion depth you use for DensityPlot[] can have a significant effect on the approximation. For instance, if you try to use

    Code:
    DensityPlot[HeavisideTheta[x, y], {x, -1, 1}, {y, -1, 1}]

    you'll see that the resulting density plot exhibits approximation artifacts around the jump discontinuity in the 2D Heaviside function. This can be controlled using PlotPoints and MaxRecursion, but you should be aware that DensityPlot[] is more suited to functions that are at least [itex]C^1[/itex] everywhere in the region of interest.
 
  • #3
1,868
0
There are two main reasons why you'll have trouble with attempting to plot functions like that using DensityPlot.

The function is zero everywhere but at a single point. DensityPlot[] is unsuitable for functions like this, as a few minutes thought should convince you.

You are correct; I need to think a little more. Thanks; I'll try out your suggestions.
 

Related Threads on Mathematica: Plotting piecewise functions

Replies
5
Views
2K
Replies
5
Views
4K
Replies
2
Views
2K
Replies
3
Views
1K
Replies
2
Views
2K
Replies
0
Views
4K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
1
Views
2K
Top