Introduction - If you have any usage issues, please Google them yourself
input: phi is the objective function, a, b is the search interval of two endpoints
delta, epsilon, respectively is the independent variable and function value of allowable error
Output: s, phis is approximate minimum point and minimum value respectively, G is n times 4 matrix