How to solve a 2D PDE with backward Euler?

6 ビュー (過去 30 日間)
Danny Helwegen
Danny Helwegen 2019 年 1 月 13 日
編集済み: M 2019 年 1 月 14 日
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 Center および File Exchange で Boundary Conditions についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by