My talk is composed of three parts. The first part is on high order residual distribution (RD) schemes for steady state hyperbolic conservation laws. High order RD schemes are conservative schemes that overcome the restriction of mesh sizes in high order finite difference schemes, and yet have comparable computational costs. It has a broad range of applications from Navier-Stokes equations to semiconductor simulations. I will present the design of the scheme, a Lax-Wendroff type theorem and the numerical results. In the second part, I will discuss the applications in systems biology. The modeling of the two biological systems--cell polarization and multi-stage cell lineages, and the computational aspect will be discussed. New efficient numerical schemes for both time evolution and steady state reaction-diffusion equations that arise in many biological models will be presented in the third part.