The secant method is a root-finding algorithm (like the Newton-Raphson method) for solving an equation of the form f(x) = 0. Given two initial values x0 and x1, the successive terms are defined recursively as follows.

x_n =x_n?1 ?f(x_n?1)*(x_n?1 ?xn?2)/(f(x_n?1) ? f(x_n?2)) , n=2,3,....

Implement the secant method as a R function in the following form

secant <- function(f,x0,x1,n,err) {...} where f is a function to evaluate f(x), x0 and x1 are the values of x₀ and x₁ respectively, and err is the error tolerance in absolute value of f(x). Try to solve

3e^?y0.5 + 3e^?y1.0 + 3e^?y1.5 + 3e^?y2.0 = 98.39 using your function.