function y=my_func(x)
%This function returns the value of the function: 1/(1+3*cos(x)*cos(x))
%at a given value of x
