使用 cvxpy 优化 M/M/k/k 排队系统

我想制定一个排队网络（M/M/k/k）优化问题，以服务器数量（k）为cvxpy中的变量。由于目标函数（即阻塞概率）是凸的，所以我很乐意在 cvxpy 中解决它，但是，没有可用的原子函数来帮助我。

我意识到有一个在this post上设计排队系统的例子，当变量变成服务器数量时，就很难解决了。这个问题有点类似于cvxpy website的优化问题（5），不同之处在于到达/服务率是恒定的，而服务器的数量是优化变量。

如何在cvxpy中实现这个问题？任何想法或提示表示赞赏。