Plot PDE triangular mesh
pdemesh(p,e,t) pdemesh(p,e,t,u) h = pdemesh(p,e,t) h = pdemesh(p,e,t,u)
pdemesh(p,e,t,u) plots PDE node or triangle data u using a mesh plot. If u is a column vector, node data is assumed. If u is a row vector, triangle data is assumed. This command plots substantially faster than the pdesurf command.
The geometry of the PDE problem is given by the mesh data p, e, and t. For details on the mesh data representation, see initmesh.
This command is just shorthand for the calls
If you want to have more control over your mesh plot, use pdeplot instead of pdemesh.
Plot the mesh for the geometry of the L-shaped membrane.
[p,e,t] = initmesh('lshapeg'); [p,e,t] = refinemesh('lshapeg',p,e,t); pdemesh(p,e,t)
Now solve Poisson's equation –Δu = 1 over the geometry defined by the L-shaped membrane. Use Dirichlet boundary conditions u = 0 on ∂Ω, and plot the result.
u = assempde('lshapeb',p,e,t,1,0,1); pdemesh(p,e,t,u)