plot y=x+arcsin(1.5sin ( pi/3-arcsin(sinx/1.5)))-pi/3 for 0.4875<x<pi/2