In order to compute linear wave transformation over a single step bottom, both variational approximation and eigenfunction expansion method are used. Both numerical results are in good agreement for reflection and transmission coefficients, surface displacement respectively. However x velocity profiles at the boundary of step are seen to be different to each other even though x velocity matching condition is used.