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
|
import graph;
real f(real x) {return x^3-x+2;}
pair F(real x) {return (x,f(x));}
void rectangle(real a, real b, real c, real h(real,real))
{
real height=(a < c && c < b) ? f(c) : h(f(a),f(b));
pair p=(a,0), q=(b,height);
path g=box(p,q);
fill(g,lightgray);
draw(g);
}
void partition(real a, real b, real c, real h(real,real))
{
rectangle(a,a+.4,c,h);
rectangle(a+.4,a+.6,c,h);
rectangle(a+.6,a+1.2,c,h);
rectangle(a+1.2,a+1.6,c,h);
rectangle(a+1.6,a+1.8,c,h);
rectangle(a+1.8,b,c,h);
draw((a,0)--(F(a)));
draw((b,0)--(F(b)));
draw(graph(f,a,b,operator ..),red);
draw((a,0)--(b,0));
labelx("$a$",a);
labelx("$b$",b);
}
|
