return (
        amp *
            (((1 - aa) * Math.sin(phase)) /
                (1 + aa - a2 * Math.cos(beta * phase))) +
        dc