var('x, y')
g=sqrt(-log(x));
R=plot(sqrt(-log(x)), [0, 1]);
S=plot(x, [0, 1.8], linestyle='--', rgbcolor=(1,0,0));
T1=polygon([(0,0),(0,2)]+[(0,2),(e^(-4),2)]+[(y,g.subs_expr(x==y)) for y in [e^(-4),e^(-4)+0.05,..,1]]+[(1,0),(0,0)], alpha=0.1);
show(R+S+T1, aspect_ratio=1, figsize=[5,5], axes_labels=('y','x'));
reset()