A numerical method of hydrodynamic bearing analysis is presented which is simple in concept, yet capable of development to handle complex situations such as dynamic misalignment. It is similar to the finite difference solution of Reynolds equation, but incorporates a more realistic modelling of cavitation. The approach to a numerical solution is direct, and should facilitate a better ‘feel’ for the way in which the physical processes are modelled. Results produced with this analysis are compared with other published data for aligned crankshaft bearings and misaligned sterntube bearings.