How to input crank nicolson into matlab?

<<

>>

I am at a loss on how to code these to solve in the Crank Nicolson equation. I have solved the equations, but cannot code it into matlab.

Answer 1

[Conrad Suen]

Answer 2

[Conrad Suen]

This is the code and errors I am getting.

Answer 3

[simopera6]

Hi Conrad, If you are trying to solve by crank Nicolson method, this is not the way to do it. You have to solve it by tri-diagonal method as there are minimum 3 unknowns for the next time step. There are many videos on YouTube which can explain this. Hope this helps. If you can post a code after doing this, we can have a look at it. Best,

Answer 4

[Bjorn Gustavsson]

Crank-Nicolson is for partial differential equations, this seems to be ordinary differential equations. Perhaps try matlabs ode23 et al.

HTH