A computational approach to quasi-static propagation of brittle fracture