Several methods have been proposed in the literature for solving the Black-Scholes equation for European Options. The method proposed in the current study achieves spectral accuracy in both space and time. The method is based on minimization of a functional given in terms of the sum of squares of the residuals in the partial differential equation and initial condition in different Sobolev norms, and a term which measures the jump in the function and its derivatives across inter-element boundaries in appropriate fractional Sobolev norms. To obtain values of the solution and its derivatives the initial condition is mollified and the computed solution is post processed. Error estimates are obtained for this method. Specific numerical examples are given to show the efficiency of this method. © 2015 Elsevier Ltd.