工程结构优化设计编程作业1. 用黄金分割法求方程在区间 -1,1 上的解。在 MATLAB 中没有专门的函数实现黄金分割法求解线性方程,可通过编写glodf.m 函数实现黄金分割法求解,其代码如下:function x=glodf(f,a,b,eps) if(nargin=3)eps=1.0e-4; endf1=subs(sym(f),findsym(sym(f),a);f2=subs(sym(f),findsym(sym(f),b); if(f1=0)x=a; end if(f2=0)x=b; end if(f1*f20)disp(两端点函数值乘积大于0! ); return;elset1=a+(b-a)*0.382; t2=a+(b-a)*0.618;f_1=subs(sym(f),findsym(sym(f),t1);f_2=subs(sym(f),findsym(sym(f),t2); tol=abs(t1-t2);while(toleps) if(f_1*f_20)a=t1