def METHOD_NAME(self):
    """
    Case from R. Simpson's PhD thesis.
    His wing tip displacements: 15.627927627927626, 3.3021978021978025
    I always get higher deflection when using my gravity implementation
    instead of his. His results with gravity are not validated.***
    :return:
    """
    import sharpy.sharpy_main
    solver_path = os.path.abspath(os.path.dirname(os.path.realpath(__file__)) +
                                  '/smith_g_4deg/smith_g_4deg.sharpy')
    sharpy.sharpy_main.main(['', solver_path])
    # read output and compare
    output_path = os.path.dirname(solver_path) + '/output/smith_g_4deg/WriteVariablesTime/'
    pos_data = np.genfromtxt(output_path + 'struct_pos_node20.dat')
    self.assertAlmostEqual((pos_data[2] - 15.55)/15.55, 0.00, 2)
    self.assertAlmostEqual((pos_data[3] - 3.671)/3.671, 0.00, 2)
    # results:
    # N = 10 elements
    # M = 15 elements
    # full wake:
    # Nrollup = 100
    # Mstar = 80
    # pos last node of the wing
    # 0.05668 15.5422 3.53971
    # total forces:
    # tstep | fx_st | fy_st | fz_st
    #   0   | 3.229 | -1.059e-3| 3.766e2
    # 7500 seconds