function rhs = ic_variation_rhs(t,y)
rhs = sin(y);