1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
"""
How to fill an area
===================
"""
# sphinx_gallery_thumbnail_number = 2
# %%
# %%
# In this example, we show how to fill specified area with a given color.
# %%
import openturns as ot
from openturns import viewer
from matplotlib import pylab as plt
ot.Log.Show(ot.Log.NONE)
# %%
# We generate a sample from the standard Normal distribution.
# %%
dist = ot.Normal()
graph = dist.drawPDF()
view = viewer.View(graph)
# %%
def linearSample(xmin, xmax, npoints):
"""Returns a sample created from a regular grid
from xmin to xmax with npoints points."""
step = (xmax - xmin) / (npoints - 1)
rg = ot.RegularGrid(xmin, step, npoints)
vertices = rg.getVertices()
return vertices
# %%
a = 1.0
b = 2.0
# %%
nplot = 100 # Number of points in the plot
x = linearSample(a, b, nplot)
y = dist.computePDF(x)
# %%
# Compute the bounds to fill: the lower vertical bound is 0 and the upper vertical bound is the PDF.
# %%
vLow = [0.0] * nplot
vUp = [y[i, 0] for i in range(nplot)]
# %%
area = dist.computeCDF(b) - dist.computeCDF(a)
# %%
boundsPoly = ot.Polygon.FillBetween(x.asPoint(), vLow, vUp)
graph = dist.drawPDF()
graph.add(boundsPoly)
graph.setTitle("Area = %.3f" % (area))
graph.setLegends([""])
view = viewer.View(graph)
plt.show()
# %%
# The CDF difference is equal to the area under the curve.
|
