I am new to differential equations, any help would be great.(adsbygoogle = window.adsbygoogle || []).push({});

I have a ODE of the second order u''x = e^x over the domain [1, 1] where u'(0) = 0 is a Neumann boundary on the ODE. I am trying to approximate the solution using the finite differences method, I can do Dirichlet boundaries with finite differences with no problem however the Neumann boundaries are a problem.

The second-order finite difference is

(e^(x - h) - 2*e^(x) + e^(x + h)) / h^2

where h is the computed interval (change in x) across the domain.

How can you model the approximation so that the first derivative at u'(0) = 0 is taken into account. The values I am getting are nothing like the exact solution that I have computed. I am looking to learn this procedure so can anyone point me to the algorithm for this?

Thank you.

**Physics Forums | Science Articles, Homework Help, Discussion**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# ODE with Neumann Bound

**Physics Forums | Science Articles, Homework Help, Discussion**