Using Inigo Quilez's palette function. Reference
vec3 palette(in float t, in vec3 a, in vec3 b, in vec3 c, in vec3 d) { return a + b*cos(6.28318530718*(c*t+d)); }