$\frac{d^2 \theta}{dx^2} (1 + \beta \theta) + \beta \left(\frac{d \theta}{d x}\right)^2 - m^2 \theta = 0$
Boundary Conditions $\theta=100$ at $x = 0$, $\frac{d\theta}{dx} = 0$ at $x = 2$
$\beta$ and $m$ are constants. Please help me solve this numerically (using finite difference). The squared term is really complicating things! Thank You!