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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
set terminal cairolatex pdf color standalone
set output "p-laplacian-fixed-point-p=1,5-err-lp.tex"
set size square 0.9
set colors classic
set logscale xy
set key bottom
set xrange [1e-2:1e0]
set yrange [1e-9:1e-1]
set xlabel '[c]{\large $h$}'
set label '[l]{\large $\|u-u_h\|_{0,p,\Omega}$}' at graph 0.08,0.92
set xtics (\
"[c]{$10^{-2}$}" 1e-2, \
"[c]{$10^{-1}$}" 1e-1, \
"[c]{$1$}" 1)
set ytics (\
"[r]{$10^{-8}$}" 1e-8, \
"[r]{$10^{-6}$}" 1e-6, \
"[r]{$10^{-4}$}" 1e-4, \
"[r]{$10^{-2}$}" 1e-2 )
graph_ratio = 2.0/8.0
# triangle a droite, pente +2
slope_A = graph_ratio*2.0
xA = 0.15
yA = 0.64
dxA = 0.10
dyA = dxA*slope_A
set label "[l]{\\scriptsize $2=k+1$}" at graph xA+dxA+0.02, yA+0.5*dyA right
set arrow from graph xA, yA to graph xA+dxA, yA nohead
set arrow from graph xA+dxA, yA to graph xA+dxA, yA+dyA nohead
set arrow from graph xA+dxA, yA+dyA to graph xA, yA nohead
# triangle a gauche, pente +1
slope_B = graph_ratio*3.0
xB = 0.15
yB = 0.35
dxB = 0.10
dyB = dxB*slope_B
set label "[r]{\\scriptsize $3$}" at graph xB-0.02, yB+0.5*dyB right
set arrow from graph xB, yB to graph xB, yB+dyB nohead
set arrow from graph xB, yB+dyB to graph xB+dxB, yB+dyB nohead
set arrow from graph xB+dxB, yB+dyB to graph xB, yB nohead
# triangle a droite, pente +4
slope_C = graph_ratio*4.0
xC = 0.25
yC = 0.14
dxC = 0.10
dyC = dxC*slope_C
set label "[l]{\\scriptsize $4$}" at graph xC+dxC+0.02, yC+0.5*dyC right
set arrow from graph xC, yC to graph xC+dxC, yC nohead
set arrow from graph xC+dxC, yC to graph xC+dxC, yC+dyC nohead
set arrow from graph xC+dxC, yC+dyC to graph xC, yC nohead
plot \
'p-laplacian-fixed-point-p=1.5-P1-err.gdat' \
u (1/$1):3 t '[r]{$k=1$}' w lp lw 2, \
'p-laplacian-fixed-point-p=1.5-P2-err.gdat' \
u (1/$1):3 t '[r]{$k=2$}' w lp lw 2 lc '#008800', \
'p-laplacian-fixed-point-p=1.5-P3-err.gdat' \
u (1/$1):3 t '[r]{$k=3$}' w lp lw 2
# 0.4*x, 0.2*x**2, 0.1*x**3, 0.2*x**4
#pause -1 "<retour>"
|
