I have this code. What is the best way to solve these king of trigonometric equations?
A= 1
B = 0.6
from scipy.optimize import fsolve
from scipy import optimize
def func(x):
return((((A**2) + (B**2)*np.cos(A*x)*np.cos(B*x)) + ( (A**3) *np.sinh(A*x)*np.cos(B*x)) - ((B**3) *np.sin(B*x)*np.cosh(A*x))))
root = fsolve(func, 21.5)