用matlab牛顿法求解x^2-2x-1=0

dsolve('Dy-2*y/(x+1)=(x+1)^5/2')ans=(C2*exp((2*t)/(x+1)))/4-(3*x)/2-(15*x^2)/4-5*x^3-(15*x^4)/4-(3*x

原因在这一句：elsefprintf('x0=%fxe=%fk=%d\n',x0,xe,k)其中xe=%f只能输出的精度是0.000000,默认6位仅只小数点后面包含6位.而实际xe=0.000000

esp表达错了吧,e的指数幂应该用exp表示.另外你这个n怎么定义的?我不大懂啊.给你分享个牛顿迭代法的程序吧,你参考下看能不能解决你的问题：

#include#includevoidmain(){floatx,x0,f,f1;x0=0.5;do{f=x0*x0*x0-x0*x0-1;f1=3*x0*x0-2*x0;x=x0-f/f1;x0=

源程序如下：clearclcN=100;x=2;forii=0:Nxl=x;num=ii;x=x-(x^3-3*x+1)/(3*x^2-3);ifabs(x-xl)

functionroot=NewtonRoot(f,a,b,eps)%用牛顿法求方程的一个根%方程的表达式：f%区间的左端点：a%区间的右端点;b%根的精度：eps%求得的根：rootif(nargi

我是用C得到结果:2.1155229/*牛顿迭代法解方程组的解x0为迭代的初值,n为迭代次数,jingdu为精度function为求根代数式,d2functoin为其导数返回最终符合一定精度的根*/d

functionx=NewRaph(x0,f,fprime,tol)x(1)=x0;i=1;gx=x(i)-f(x(i))/fprime(x(i));while(abs(x(i)-gx)>tol)gx

兄弟,西理工的?

采用第一个.首先你的两个代码的计算过程和方法以及步骤是一致的.只不过第二个将k==N放在循环内部判断是没有必要的.放在while外面,可以节省点计算量.如果你要求结果精度高一些的话,你调用：x=nan

贴上来,或者发到

function x = newton_method(fun,x0,e,N)if nargin<5    N=50

Newton-Raphson求解非线性方程组matlab源程序matlab程序如下：functionhom[P,iter,err]=newton('f','JF',[7.8e-001;4.9e-001

[r,n]=newton调用即可

这样吧,你不要直接通过solve求解啊,可以画一个图像啊,以x的值为横坐标,行列式值为纵坐标,或者纵坐标取一下对数（如果变化范围太大）,然后可以得出大致得到解的范围；最后可以通过其他方法计算更加精确的

x1=x-func1_1(x)/func1_1_1(x);是点除再问：Error:File:func1_1.mLine:1Column:22TheinputcharacterisnotvalidinM

是这个样子的,你应该是直接运行的m文件当然提示你x没有定义了可以在命令行输入x0=...;newtoneqs(x0);就可以了提示一下,如果你者几个函数放在一个m文件中,应该是newtoneqs(x0

定义函数functiony=nd(x)y=0.036-((x/2090.7).^(1/0.1585))-x/182000functiony=nd0(x)y=-(1/0.1585)*(x/2090.7)

程序：clcclearx=zeros(1,25);x(1)=2;formatlongfori=1:24;x(i+1)=x(i)-(x(i)^3-3*x(i)+2)/(3*x(i)^2-3);enddi

程序：A=[1 2 3;1 3 1; 0 1 2];B=[1 0; 0 1;-1 0];X