数值分析作业:matlab牛顿法、二分法和割线法 txt.txt
2006-10-15 评论:8 下载:3
下载
大小:1.68KB
数值分析作业:matlab牛顿法、二分法和割线法 txt
节选段落一:
数值分析作业:matlab牛顿法、二分法和割线法
By tian0114xiao 发表于 2006-10-12 21:10:00
http://bbs.sciei.com/ 科研中国
function y=newton(x0) %牛顿法
x1=x0-fun(x0)/dfun(x0);
n=1;
while (abs(x1-x0)>=1.0e-4) & (n<=100000000)
x0=x1;
x1=x0-fun(x0)/dfun(x0);节选段落二:
n=n+1;
end
y=x1
n
function y=fun(x)
y=x^4+2*x^2-x-3;
function y=dfun(x)
y=4*x^3+4*x-1;
function y=erfen(f,a,b,esp) %二分法
if nargin<4 esp=1.0e-4;
end
if feval(f,a)*feval(f,b)<0
n=1;
c=(a+b)/2;节选段落三:
x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0)); %根据初始XO和X1求X2
n=1;
while (abs(x1-x0)>=1.0e-4) & (n<=100000000) %判断两个条件截止
x0=x1; %将x1赋给x0
x1=x2; %将x2赋给x1
x2=x1-fun(x1)
数值分析作业:matlab牛顿法、二分法和割线法
By tian0114xiao 发表于 2006-10-12 21:10:00
http://bbs.sciei.com/ 科研中国
function y=newton(x0) %牛顿法
x1=x0-fun(x0)/dfun(x0);
n=1;
while (abs(x1-x0)>=1.0e-4) & (n<=100000000)
x0=x1;
x1=x0-fun(x0)/dfun(x0);节选段落二:
n=n+1;
end
y=x1
n
function y=fun(x)
y=x^4+2*x^2-x-3;
function y=dfun(x)
y=4*x^3+4*x-1;
function y=erfen(f,a,b,esp) %二分法
if nargin<4 esp=1.0e-4;
end
if feval(f,a)*feval(f,b)<0
n=1;
c=(a+b)/2;节选段落三:
x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0)); %根据初始XO和X1求X2
n=1;
while (abs(x1-x0)>=1.0e-4) & (n<=100000000) %判断两个条件截止
x0=x1; %将x1赋给x0
x1=x2; %将x2赋给x1
x2=x1-fun(x1)
查看更多评论 >