请在上述条件的基础上，写出一个使用最小二乘法对三维曲线进行曲线拟合，最后输出拟合的三维坐标的程序...

用最小二乘法对三维曲线进行曲线拟合的程序如下：

import numpy as np
from scipy.optimize import leastsq

def func(p, x):
    a, b, c = p
    return a*x[0]**2 + b*x[1]**2 + c*x[2]**2

def residuals(p, y, x):
    return y - func(p, x)

x = np.array([[1, 2,