這本是個作業，但因為各種原因及出錯，這個簡單的**寫了2個多小時，內心不爽。因此寫下部落格記錄。
**實現如下：
clear;
clc;
syms u v
f1 = 6 * u^3 + u * v - 3 * v^3 - 4;
f2 = u^2- 18 * u * v^2 + 16 * v^3 + 1 ;
f=[f1 f2 ];
df0=[diff(f,u);diff(f,v)];
df = df0.';
x0=[-1 -1]; %%%as shown in the book, the solution we find by newton's way is based on the first estimate
n=200;
% 求解
for i=1:n; % 運算次數
p=subs(f,,);% 為subs函式的賦值運算的矩陣
q=subs(df,,);% 為subs函式的賦值運算 2*2的矩陣
x =x0 - grass(q,p);% 得到的一組新數值矩陣 理解為 新的x0
if norm(x-x0)
Python例項 多變數非線性方程求解
coding utf 8 多變數非線性方程求解 import sympy import scipy from scipy import optimize import numpy as np from matplotlib import pyplot as plt from pylab import...
matlab利用牛頓法求解非線性方程01
在如下鏈結中，本人看到了牛頓法求解非線性方程的乙個程式，給人耳目一新的感覺。但不知何故，本人未能執行該博主寫的程式，於是我對原程式做了一定修改，程式得以執行，如下 牛頓法程式 function gen newton f,x,tol f為函式，x0為初值，tol為指定允差，如果預設，預設為1e 6 i...
MATLAB牛頓法求解非線性方程組2
牛頓法主程式 clear clc format x0 0 0 迭代初始值 eps 0.00001 定位精度要求 for i 1 10 f double subs fun x0 df double subs dfun x0 得到雅克比矩陣 x x0 f df if abs x x0 eps break...