Added legrange polynomials to equations instead of simpsons rule.
BRF plotting function. Fixed the J term problem.