Geological faults have frictional heterogeneity. Here, we investigate the slip behavior of such heterogeneous faults in the simplest situation, simulating anti-plane shear slip along an infinite linear fault in a 2D elastic model space, subjected to constant loading by external stresses characterized by stiffness and plate-motion velocity. The fault is discretized into many sub-faults, each of which has either a velocity-weakening (VW) or velocity-strengthening (VS) friction law.