System of ODE Boundary Value Problem with 2nd Order Backward Difference

Click For Summary
SUMMARY

The discussion focuses on solving a system of ordinary differential equations (ODEs) using the Second Order Backward Difference formula in MATLAB. The equations presented are: du/dx = 998u + 1998v and dv/dx = -999u - 1999v, with boundary conditions u(0) = 1 and v(0) = 0 for the interval 0 < x < 10. The user encounters difficulties with the k-2 index in the backward difference equation, suggesting the potential use of a fictitious point. Ultimately, they consider employing the Runge-Kutta 4th order (RK4) method to initiate the numerical solution.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with MATLAB programming
  • Knowledge of numerical methods, specifically the Second Order Backward Difference formula
  • Experience with the Runge-Kutta 4th order (RK4) method
NEXT STEPS
  • Implement the Second Order Backward Difference formula in MATLAB
  • Explore the use of fictitious or ghost points in numerical methods
  • Learn about the Runge-Kutta 4th order (RK4) method for solving ODEs
  • Investigate Neumann boundary conditions and their applications in numerical solutions
USEFUL FOR

Mathematicians, engineers, and students working on numerical methods for solving boundary value problems in ordinary differential equations.

teknodude
Messages
150
Reaction score
0
[tex]{\frac {{\it du}}{{\it dx}}}=998\,u+1998\,v[/tex]
[tex]{\frac {{\it dv}}{{\it dx}}}=-999\,u-1999\,v[/tex]
[tex]u \left( 0 \right) =1[/tex]
[tex]v \left( 0 \right) =0[/tex]
0<x<10
Second Order Backward Difference formula
[tex]{\frac {f_{{k-2}}-4\,f_{{k-1}}+3f_{{k}}}{h}}[/tex]

I'm trying solve this numerically in matlab, but can't seem to figure out what to do with the k-2 indice in the 2nd order backward difference equation, because it is outside the boundary. I was thinking of using a ficticious or ghost point, but I thought that only applies if a neuman boundary condition is given. The way i think of it, I have 4 unknowns and only 2 equations.

EDIT: ok after thinking about it. I think I have to use another numerical method to start the process, like RK4.
 
Last edited:
Physics news on Phys.org
I think that's a good idea xD
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 0 ·
Replies
0
Views
3K
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
0
Views
2K