The laminar-turbulent transition of spherical Couette flow for case of clearance ratio beta = 0.14 was investigated by laboratory experiment. Calculating the correlation dimension and drawing the Poincare section, it is revealed that the flow field traces a scenario as follows: steady state, periodic state, quasi-periodic state, chaos, periodic state, steady state, periodic state, and chaos. A finite-difference method for solving three-dimensional, time-dependent, incompressible Navier-Stokes equations in spherical polar coordinates is also presented. Based on a new algorithm, a higher accurate numerical code has been developed, and it is demonstrated that the initial-boundary numerical code is valid for studying the spherical Couette flow problems.