Numerically finding roots of multiple variables in an integral equation

Question

I am trying to solve for Hb and Et in the following integral equation:

All the other constants are known and J is an array. So essentially I am trying to find Hb and