如何在scala中编写圆拟合程序

import breeze.linalg._
import breeze.optimize._

object CircleFitting {
  def main(args: Array[String]): Unit = {
    // 输入数据点坐标
    val points = DenseMatrix(
      (1.0, 2.0),
      (3.0, 4.0),
      (5.0, 6.0),
      (7.0, 8.0)
    )

    // 定义圆心和半径作为优化变量
    val initialGuess = DenseVector(0.0, 0.0, 1.0)

    // 定义优化目标函数
    val objective = new DiffFunction[DenseVector[Double]] {
      def calculate(x: DenseVector[Double]): (Double, DenseVector[Double]) = {
        val center = x(0 to 1)
        val radius = x(2)
        val residuals = points(*, ::) - center
        val distances = sqrt(sum(residuals :* residuals, Axis._1))
        val error = sum(pow(distances - radius, 2))
        val gradient = DenseVector(
          -2.0 * sum(residuals(::, 0) * (distances - radius) / distances),
          -2.0 * sum(residuals(::, 1) * (distances - radius) / distances),
          -2.0 * sum(distances - radius)
        )
        (error, gradient)
      }
    }

    // 运行优化算法
    val result = LBFGS[DenseVector[Double]]().minimize(objective, initialGuess)

    // 提取优化结果
    val center = result.x(0 to 1)
    val radius = result.x(2)

    // 打印结果
    println(s"圆心坐标：$center")
    println(s"半径：$radius")
  }
}