1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
import graph3;
import palette;
size(0,300);
currentprojection=perspective(3,-2,2);
real V(real r) {return r^4-r^2;}
real V(pair pos) {return V(abs(pos));}
real R=1/sqrt(2);
real z=-0.2;
bool active(pair pos) {return abs(pos) < R;}
bool above(pair pos) {return V(pos) >= z;}
pair a=(-1.5,-1);
pair b=(0.5,1);
real f=1.2;
draw(plane(f*(b.x-a.x,0,z),(0,f*(b.y-a.y),z),(a.x,a.y,z)),
lightgrey+opacity(0.5));
surface s=surface(V,a,b,40,Spline,active);
draw(s,mean(palette(s.map(new real(triple v) {
return above((v.x,v.y)) ? 1 : 0;}),
new pen[] {lightblue,lightgreen})),black);
xaxis3(Label("$\phi^\dagger\phi$",1),red,Arrow3);
zaxis3(Label("$V(\phi^\dagger\phi)$",1),0,0.3,red,Arrow3);
|
