This matlab code simulates Rayleigh Bernard convection using the Chorin scheme from computational fluid dynamics.