y=13cos(arcsin((x/16)^1/3))-5cos(2arcsin((x/16)^1/3))-2cos(3arcsin((x/16)^1/3))-cos(4arcsin((x/16)^1/3))