MATLABA割线法求方程根的论文...明天就要交了.有谁能帮帮小弟啊!!!!

怎么帮你？
割线法函数(求解一元方程,MATLAB)


%x(k+1)=x(k)-(f(x(k))*(x(k)-x(k-1)))/(f(x(k))-f(x(k-1)))
function y=ger(x0,x1)
x2=x1-fc(x1)*(x1-x0)/(fc(x1)-fc(x0));
n=1;
while(abs(x1-x0)>=1e-4)&(n<100000000)
    x0=x1;
    x1=x2;
    x2=x1-fc(x1)*(x1-x0)/(fc(x1)-fc(x0));
    n=n+1;
end
x2
n
%用割线法求方程f(x)=x^3-3*x-1=0在x0=2附近的根。误差限为1e-4,取x0=2,x1=1.9。
%解：编制函数文件：
%fc.m
%function y=fc(x)
%y=x.^3-3*x-1;
%>>ger(2,1.9)
%x2=
%   1.8794
%n=
%   4

fc.m
function y=fc(x)
y=x.^3-3*x-1;
这样可以吗

function y=gexian(x0,x1) %割线法

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)*(x1-x0)/(fun(x1)-fun(x0)); %迭代运算

n=n+1;

end

y=x2

n

function y=fun(x)

y=x^4+2*x^2-x-3;

%y=1-x-sin(x);
这些网上都有

关键是我们那神经病老师要我们写论文...我头都晕了..哪位大哥有经验的帮帮小弟吧..跪求了

http://60.12.106.38/viewthread.php?tid=37390&highlight= 论坛上的这个链接有24篇论文可是我权限不够.哪位好心人帮我下一下啊.我QQ7067643 !!!