def METHOD_NAME(m):
    return sum(
        h[i]
        * sum(
            (
                alpha1 * (m.c[i, j] - c_des) ** 2
                + alpha2 * (m.t[i, j] - t_des) ** 2
                + alpha3 * (m.u[i] - u_des) ** 2
            )
            * m.a[j, ncp - 1]
            for j in CP
        )
        for i in range(2, nfe)
    ) + h[0] * sum(
        (
            alpha1
            * (
                (m.c_init_var + time * h[0] * sum(m.a[k, j] * m.cdot[0, k] for k in CP))
                - c_des
            )
            ** 2
            + alpha2
            * (
                (m.t_init_var + time * h[0] * sum(m.a[k, j] * m.tdot[0, k] for k in CP))
                - t_des
            )
            ** 2
            + alpha3 * (m.u[0] - u_des) ** 2
        )
        * m.a[j, ncp - 1]
        for j in CP
    )