  The exponential B-spline basis function set is used to develop a collocation
method for some initial boundary value problems (IBVPs) to the Gardner
equation. The Gardner equation has two nonlinear terms, namely quadratic and
cubic ones. The order reduction of the equation is resulted in a coupled system
of PDEs that enables the exponential B-splines to be implemented. The system is
integrated in time by Crank-Nicolson implicit method. The validity of the
method is investigated by calculating the discrete maximum error norm and
observing the absolute relative changes of the conservation laws at the end of
the simulations.
