double yf (double a, double b){
                               y =(exp(x)+exp(-x))/2;
                               return y;
                               }