Problem $3$

A function $u$ is defined in the domain $($shaded$)$

$u=\sqrt[\phantom{2}]{x^2+(y-3{)}^2}$

Write a MATLAB program to draw a contour plot of $u$ in the domain, as shown