Methods 100 Implementation (used by Prof. Dr. Bernd Fischer to solve the problem) 1. Initialization:   Start with:   0 = u Equation 13: Initialization 2. Iteration Force field: ( ) ( ) ( ) (  ) (  ) (  )  grad (  ) T R T = - - - - f  u x    u  x x x    u  x Equation 14: Force   Finite difference approximation of the partial differential equation: ( ) grad div (  ) new new m l m a + + = u u f  u Equation 15: Finite differ- ence approximation Update u: new = u    u Equation 16: Update dis- placement field u 3. End