3146.sagews

from sage.plot.polygon import Polygon
g = Graphics();
x = var('x');
radius = 1;
g = plot(sqrt(radius**2 - x**2), xmin=-1, xmax=1, rgbcolor=(0,.8,.2));
g += polygon2d([(0,0), (1,0), (radius*cos(pi/3), radius*sin(pi/3))], rgbcolor=(.8,0,.2));
g += polygon2d([(0,0), (radius*cos(pi/3), radius*sin(pi/3)), (radius*cos(2*pi/3), radius*sin(2*pi/3))], rgbcolor=(.85,0,.22));
g += polygon2d([(0,0), (radius*cos(2*pi/3), radius*sin(2*pi/3)), (radius*cos(pi), radius*sin(pi))], rgbcolor=(.87,0,.25));
show(g, aspect_ratio=1);