Question

关于Matlab函数调用，我是初学者，求大神赐教呀

function [x_star,k]=bisect1(fun,a,b,ep)
%主函数，每次运行时要输入fun=inline('x^3+x-1');[x_star,k]=(fun,1,1.5,0.0005)
%但是老师要求编一个M文件直接就得出结果
if nargin<4 ep=1e-5;end
fa=feval(fun,a);fb=feval(fun,b);
if fa*fb>0 x_star=[fa,fb]; k=0;return;end
k=1;
while abs(b-a)/2>ep
x=(a+b)/2;fx=feval(fun,x);
if fa*fx<0
b=x;fb=fx;
else
a=x;fa=fx;
end
k=k+1;
end
x_star=(a+b)/2;
%我编写的文件如下：
function y=fun(x);y=('sin(x)-x^2/4');[x_star,k]=bisect1(fun,1.5,2)
%运行时我直接输入fun ,j结果为：
fun=inline function;
fun(x)=x^3-x-1

Answer 1

不明白你为啥要写下面的m文件。是不是这样：
function [x_star,k]=bisect1(fun,a,b,ep)
%主函数，每次运行时要输入fun=inline('x^3+x-1');[x_star,k]=(fun,1,1.5,0.0005)
%但是老师要求编一个M文件直接就得出结果
if nargin<4 ep=1e-5;end
fa=feval(fun,a);fb=feval(fun,b);
if fa*fb>0 x_star=[fa,fb]; k=0;return;end
k=1;
while abs(b-a)/2>ep
x=(a+b)/2;fx=feval(fun,x);
if fa*fx<0
b=x;fb=fx;
else
a=x;fa=fx;
end
k=k+1;
end
x_star=(a+b)/2;

然后再命令窗口里写
fun=inline('x^3+x-1');
[x_star,k]=(fun,1,1.5,0.0005)
就好了吧。

结果：
>> [x_star,k]=baidu7(fun1,1,1.5,2)
x_star =
1.0000 3.8750
k =
0