message: Add example for tf.linalg.triangular_solve
diff: @@ -52,5 +52,34 @@ The output is a tensor of shape `[..., M, K]`. If `adjoint` is
  If `adjoint` is `False` then the strictly then the innermost matrices in
  `output` satisfy matrix equations
  `adjoint(matrix[..., i, k]) * output[..., k, j] = rhs[..., i, j]`.
  +
  +Example:
  +```python
  +
  +a = tf.constant([[3, 0, 0, 0],
  + [2, 1, 0, 0],
  + [1, 0, 1, 0],
  + [1, 1, 1, 1]], dtype=tf.float32)
  +
  +b = tf.constant([[4],
  + [2],
  + [4],
  + [2]], dtype=tf.float32)
  +
  +x = tf.linalg.triangular_solve(a, b, lower=True)
  +x
  +# <tf.Tensor: id=257, shape=(4, 1), dtype=float32, numpy=
  +# array([[ 1.3333334 ],
  +# [-0.66666675],
  +# [ 2.6666665 ],
  +# [-1.3333331 ]], dtype=float32)>
  +
  +a@x
  +# <tf.Tensor: id=263, shape=(4, 1), dtype=float32, numpy=
  +# array([[4. ],
  +# [2. ],
  +# [4. ],
  +# [1.9999999]], dtype=float32)>
  +```
  END
  }
  