|
- function ff=gold(f,a,b,t)
- tic
- x=a+0.382*(b-a);
- y=a+0.618*(b-a);
- f1=f(x);
- f2=f(y);
- count=0;
- while abs(a-b)>t
- if f1<f2
- b=y;
- y=x;
- f2=f1;
- x=a+0.382*(b-a);
- f1=f(x);
- else
- a=x;
- x=y;
- f1=f2;
- y=a+0.618*(b-a);
- f2=f(y);
- end
- count=count+1;
- end
- x=0.5*(a+b),f=f(x),count
- toc
复制代码
==========================================
命令窗口中执行:
f=@(x) (x.^2-1).^2;
gold(f,0,2,.1)
==========================================
关于上面你的错误,最好先看看最基本的MATLAB有关M文件的内容 |
|