clear all;clc %% the original conditions f=@(x)exp(-x)+x^2; a=0; b=1; Theta_error=0.15; %% execute FibonacciOpt Algorithm [x_opt,f_opt,stepNum]=FibonacciOpt(f,a,b,Theta_error); %% function figure x=a:0.01:b; y=exp(-x)+x.^2; plot(x,y,'k')