微分方程数值解作业 2

Problem 1

To solve (2.4)

We can always obtain the equation at non-boundary positions,

The above FDE approximation has an error of

Question (a)

The first condition

and in Eqn. becomes

For the second condition

substitute with an one-sided approximation

we have

And the final equation becomes

The whole FDE system in matrix form

Question (b)

Therefore

FDE system in matrix form

Problem 2

Question (a)

Use central difference to approximate with

Substitute back

Question (b)

Given the differerntial equation

Simply susbtitute with in the above equation, and discard the higher order terms, we have

Question (c)

When

we will be able to seperate from , and the difference equation becomes

Collect terms on the left hand side

The above equation is linear, and could be written in matrix form. Denote

Assume the boundary conditions are and , we can obtain Equation at and

Finally, the FDE system is expressed in explicit matrix form

Question (d)

Therefore, when satisfies , Equation (2.4) will not contain term, and can be solve by the method in Question (b). That requires to be a constant.

Problem 3

At , use a third-ordered forward difference to approximiate and

And use a forth-ordered forward difference to approximiate

Substitute back to the equation we have one equation