I have one surface for which I have an analytical solution, and another represented by a list of 3D points. My question is: how can I plot the part of one that lies outside the other?
The first surface:
R = 1;
Ω = Disk[{0, 0}, {R, R}];
θ = 1;
G = 1;
Φ = 1/2 G θ R^2 (1 - x^2/R^2 - y^2/R^2);
Plot3D[Φ, {x, y} ∈ Ω,
PlotStyle -> None,
PlotTheme -> "Detailed", Mesh -> {25},
AxesLabel -> {"x", "y", "ϕ(x,y)"},
LabelStyle -> Directive[FontFamily -> "Courier New"]]
The second surface:
circle =
{{(2 - Sqrt[3])/2, 0.5, 0}, {0, 1, 0},
{(2 - Sqrt[3])/2, 1.5, 0}, {0.5, (2 - Sqrt[3])/2, 0},
{0.5, 0.5, 0.3774047358083551`}, {0.5, 1, 0.4599364905389034`},
{0.5, 1.5, 0.3774047358083551`}, {0.5, (2 + Sqrt[3])/2, 0},
{1, 0.`, 0}, {1, 0.5, 0.4599364905389034`},
{1, 1, 0.5849364905389033`}, {1, 1.5, 0.4599364905389034`},
{1, 2, 0}, {1.5, (2 - Sqrt[3])/2, 0},
{1.5, 0.5, 0.3774047358083551`}, {1.5, 1, 0.4599364905389034`},
{1.5, 1.5, 0.3774047358083551`}, {1.5, (2 + Sqrt[3])/2, 0},
{(2 + Sqrt[3])/2, 0.5, 0}, {2, 1, 0}, {(2 + Sqrt[3])/2, 1.5, 0}}
ListPlot3D[circle,
PlotStyle -> None,
PlotTheme -> "Detailed",
Mesh -> {3, 3},
AxesLabel -> {"x", "y", "ϕ(x,y)"},
InterpolationOrder -> 2,
LabelStyle -> Directive[FontFamily -> "Courier New"]]
kor
? $\endgroup$Φ
over the intersection of the domains? Something else, perhaps? $\endgroup$