In this talk, we present the numerical analysis of a bond-based peridynamic model. Model is nonlinear and is characterized by two wells. One well, near zero strain, correspond to the linear elastic interaction whereas second well corresponds to the softening of material. We present finite element analysis when the solution is in C1 and C2. We also analyze the finite difference approximation when the solution is in Hölder space. We consider Hölder space as it allows fields which do not have well-defined derivatives. We show O(h2/ε2) rate of convergence when the solution is in C2 and O(hγ/ ε2) when the solution is in Hölder space. h is the size of mesh, ε is the size of the horizon, and γ є (0,1] is the Hölder exponent. We show that the instability in the peridynamic model is due to softening of bonds and occurs when sufficiently large number of bonds have strain above critical strain. We also consider the model in one dimension. We show that the solution of peridynamic model converges to the solution of elastodynamic. The same is shown for the approximation of peridynamic model. We present some numerical results to verify the claims. This is a joint work with Dr. Robert Lipton.