plot 2 abs(sec(x/2)) sin^4(x/4) , plotrange((0,pi/4), (0,0.004)), aspectratio 0.5