How to solve a 2D PDE with backward Euler?

12 次查看(过去 30 天)
Hi, i have to solve the 2D heat equation:
∂T/∂t = α∇^2 T = α(∂^2T/∂x^2 + ∂^2T/∂y^2)
It is given that at the 4 boundaries the T is 0. I've already written some code, but it doesn't works. Can someone review my code and help me? The error occurs at the end, when the new T are put in a matrix. The code has been added as an attachment.

采纳的回答

M
M 2019-1-14
编辑:M 2019-1-14
At line 58, you are trying to access the 11th element of the Number matrix, but its dimensions are 10 x 10, so you get an error.

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Geometry and Mesh 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by