Numerically difficult since intermediate values end up close to a local maximum of cos