代码供参考
f = @(x) [x(1)-0.7*sin(x(1))-0.2*cos(x(2)); x(1)-0.7*cos(x(1))+0.2*sin(x(2)) ]; x0 = [0,0];%初值不同结果也会有差异 fsolve(f,x0)
结果:
ans = 0.588460884110645 -0.031023352250812