Complete writing DirectForcingSolver.

 * Write functions generateA, generateBNQ and calculateRN.
 * Bug-fix: Change > to >= in function outside.
 * Bug-fix: Calculate coeffsX and coeffsY correctly.
 * Add constructor to DirectForcingSolver.
 * Provide the option to choose linear or no interpolation.