ⅆ2ⅆ⁡z2⁡cs⁡(z,k)=(1+k′2)⁢cs⁡(z,k)+2⁢cs3⁡(z,k).