leg.py now has correct J formula, includes area, and works fine with good plot of R_load vs P