def myTwistFunction(Epsilon):
    """User-defined function describing the variation of twist as a function
    of the leading edge coordinate."""
    RootTwist = 0
    TipTwist = -2
    return RootTwist + Epsilon * TipTwist