MATLAB: Geometry function – domain sketch

Question

Hi, I am trying to sketch the domain for PDE solver. It is similar to the example function CIRCLEFUNCTION.m.

My domain has four edges. Left and right edges are parallel to the y-axis. Top and bottom edges are curves:

I have following so far:

________________________

function [x,y]= dom(bs,s)

if nargin = = 0

x=4;

end

if nargin = = 1

A=[ 0, 1, 1, 0 ;

1, 1, 0 , 0 ;

1, 1, 1, 1 ;

0, 0, 0, 0]

I don't know where to go from there. I read examples online, but they are complicated.

You could choose curves connecting top vertices and bottom vertices to be simple quadratic curves, vertices are (0,0), (1,0), (1,1), (0,1).

Thank you!

Sandy

Answer 1

Try pdepoly

x = 0:0.1:1;
y = x.^2-x+1;
x = [0 x 1];
y = [0 y 0];
pdepoly(x,y)