YES
O(n^2)
TRS:
 {
             sqr(0()) -> 0(),
             sqr(s()) -> s(),
             terms(N) -> cons(recip(sqr(N))),
             dbl(0()) -> 0(),
             dbl(s()) -> s(),
          add(0(), X) -> X,
          add(s(), Y) -> s(),
        first(0(), X) -> nil(),
  first(s(), cons(Y)) -> cons(Y)
 }
 DUP: We consider a non-duplicating system.
  Trs:
   {
               sqr(0()) -> 0(),
               sqr(s()) -> s(),
               terms(N) -> cons(recip(sqr(N))),
               dbl(0()) -> 0(),
               dbl(s()) -> s(),
            add(0(), X) -> X,
            add(s(), Y) -> s(),
          first(0(), X) -> nil(),
    first(s(), cons(Y)) -> cons(Y)
   }
  Matrix Interpretation:
   Interpretation class: triangular
           [X3]  [X1]    [1 0][X3]   [1 0][X1]   [1]
   [first]([X2], [X0]) = [0 0][X2] + [0 0][X0] + [0]
   
           [0]
   [nil] = [0]
   
         [X3]  [X1]    [1 0][X3]   [1 0][X1]   [1]
   [add]([X2], [X0]) = [0 0][X2] + [0 1][X0] + [1]
   
         [X1]    [1 1][X1]   [1]
   [dbl]([X0]) = [0 0][X0] + [1]
   
         [0]
   [s] = [1]
   
         [0]
   [0] = [1]
   
           [X1]    [1 1][X1]   [1]
   [terms]([X0]) = [0 0][X0] + [0]
   
         [X1]    [1 1][X1]   [0]
   [sqr]([X0]) = [0 1][X0] + [0]
   
           [X1]    [1 0][X1]   [0]
   [recip]([X0]) = [0 0][X0] + [0]
   
          [X1]    [1 0][X1]   [0]
   [cons]([X0]) = [0 0][X0] + [0]
   
   
   Qed