声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1106|回复: 2

[编程技巧] [求助] 跪求大神解答!!关于二分法和牛顿法如何编入同一个M文件

[复制链接]
发表于 2013-10-5 11:28 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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文件中时,就不知道怎么编主函数了,跪求解答啊!!!!十分感谢
回复
分享到:

使用道具 举报

 楼主| 发表于 2013-10-5 22:53 | 显示全部楼层
没有大神来帮帮忙吗?
 楼主| 发表于 2013-10-6 19:54 | 显示全部楼层
求人帮忙啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 23:49 , Processed in 0.053847 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表