f(x) = sin(x)
@interact
def _(width=pi/12,start=0,end=pi,RAM=['None','Left','Right']):
show(f)
graph = plot(f,(x,start,end), color='red')
if RAM == 'Left':
graph += LRAMarea(f,width,start,end)
print 'Approximate area:',LRAM(f,width,start,end).n()
elif RAM == 'Right':
graph += RRAMarea(f,width,start,end)
print 'Approximate area:',RRAM(f,width,start,end).n()
if RAM != 'None': print 'Actual area:',integral(f,x,start,end).n()
show(graph,aspect_ratio=1)