[wenroutianshi]

a=1
b=7
while b-a>0.01
a1=a+0.382*(b-a)
a2=a+0.618*(b-a)
if fun1(a1)<fun1(a2)
b=a2
a2=a1
a1=a+0.382*(b-a)
else
a=a1
a1=a2
a2=a+0.618*(b-a)
end
x=(b+a)/2
end
y=fun1(x)