Question

The same governing equation for the temperature distribution with time on a 2D

square plate measuring  1  unit by 1 unit is given as

 

∂T/∂t = ∂²T/∂x² + ∂²T/∂y² ,

 

In this case, the boundary conditions are given as the Dirichlet type for 3 sides of the plate and reflected as follows (Fig. 2),

 

0  ≤  x  ≤  1.0, y = 0,            T = 0.0

0  ≤  x  ≤  1.0, y = 1.0,         T = 1.0

0  ≤  y  ≤  1.0, x = 0,            T = 0.0

 

and the Neumann boundary condition for

 

0  ≤  y  ≤  1.0, x = 1.0,

 

is given as

 

∂T/∂x = 0.0.

 

Obtain the temperature contour plot on the square plate with time, say at=0.01, 0.1 and at steady state. (You can provide contours at other times too to depict the convergence of the results at steady state.)

How to write a Matlab code. I'm a newbie