马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
急求大神帮忙啊!!!!!如何把二分法和牛顿法编入同一个m文件中进行调用?????以下是代码
二分法:function [x,err,yc]=resecm(f,a,b,eps)
ya=feval('li6_1fun',a);
yb=feval('li6_1fun',b);
if yb==0;
x=b;
end
if ya*yb<0,
disp('[a,b]是有根区间');
end
max1=1+round((log(b-a)-log(eps))/log(2));
for k=1:max1
x=(a+b)/2;
yc=feval('li6_1fun',x);
if yc==0;
a=x;
b=x;
break;
elseif yb*yc>0
b=x;
yb=yc;
else a=x;
ya=yc;
end
if (b-a)<eps,
break;
end
end
k
x=(a+b)/2
err=abs(b-a)
yc=feval('li6_1fun',x)
function y=li6_1fun(x);
y=x^3-x-1;
牛顿法:function [ x,err,k,y ] = Newtoniter(f,df,x0,tol,n)
y(x1)=x0;
for k=1:n
x=x0-feval('li6_4fun',x0)/feval('dli6_4fun',x0);
err=abs(x-x0);
x0=x;
if(err<tol)|(y==0),
break;
end
end
function y=li6_4fun(x);
y=x^3-3*x+2;
function y=dli6_4fun(x);
y=3*x^2-3;
这两个程序单独为一个M文件时没问题,但是放到同一个M文件中时,就不知道怎么编主函数了,跪求解答啊!!!!十分感谢
|