结合上述部分,以下是完整的牛顿迭代法代码实现: importnumpyasnp# 定义方程的函数defF(X):x,y=Xreturnnp.array([x**2+y**2-4,x**2-y-1])# 定义雅可比矩阵defJ(X):x,y=Xreturnnp.array([[2*x,2*y],[2*x,-1]])# 初始化猜测值X0=np.array([1.0,1.0])# 初始猜测值tolerance=1e-6max_...
牛顿法求解非线性方程组 importnumpyasnpimportmathfromsympyimport*defsolve_NR(U,X):# define symbol variablex=symbols('x')y=symbols('y')u1=U[0,0]u2=U[1,0]# define objective functionf1=x*cos(y)+y**2-u1*cos(x)-1f2=x-y+3*u2+3# take derivative of x_1 and x_2grad_11=dif...