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