A Higher Order Iterative Method for Computing the Drazin Inverse

Algorithm 2

h  = (2. −  0)/(n + 1); Subscript[t, 0]  = 0; Subscript[w, 0]  =

         E∧3; Subscript[w, n + 1]  = Exp[−3];

Table[N[(Subscript[w, i −  1]  −  2 Subscript[w, i]  + Subscript[w, i + 1])/

      h∧2 −  3 Subscript[w, i]  +

     2 (Subscript[w, i + 1]  −  Subscript[w, i −  1])/(2 h) == 0], {i, 1, n}];