How to construct the coefficient matrix(Pentagonal matrix) of five-point difference scheme?
6 views (last 30 days)
Show older comments
I want to solve an elliptic partial differential equation.
Step 1: I've the five difference scheme for this equation as follows
.
Then this scheme ends up with a linear system
, where A is a coefficient matrix.
Step 2: Construct A matrix with Matlab.
The matrix has the following form
So, how to construct the pentagonal matrix? I do not know how to write elements in pentagonal form, please give me some advice.
10 Comments
Accepted Answer
John D'Errico
on 8 Apr 2020
Ok. I'll give you a hint. Constructing a pentagonal matrix STARTS with the function spdiags. A real virtue of spdiags is you will end up with a sparse pentagonal matrix, which really is what you want when you are solving problems of that class anyway.
So do some serious reading in the help for spdiags. Try out the examples. Try out a few of your own. It is easier than it looks.
0 Comments
More Answers (0)
See Also
Categories
Find more on Startup and Shutdown in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!