This guy rocks. F_j=A[-n ̂_j p+(n ̂_j sin⁡(alpha)-V̂_rel)(τ/cos(alpha)) ]