SL4

                             
                              plot2d(16*sin(x)^3,13*cos(y)-5*cos(2*y)-2*cos(3*y)-cos(4*y))