この作品の著作権者である私は、この作品を以下のライセンスで提供します。

このファイルはクリエイティブ・コモンズ 表示-継承 3.0 非移植、2.5 一般、2.0 一般、1.0 一般のライセンスのもとに利用を許諾されています。

あなたは以下の条件に従う場合に限り、自由に

• 共有 – 本作品を複製、頒布、展示、実演できます。
• 再構成 – 二次的著作物を作成できます。

あなたの従うべき条件は以下の通りです。

• 表示 – あなたは適切なクレジットを表示し、ライセンスへのリンクを提供し、変更があったらその旨を示さなければなりません。これらは合理的であればどのような方法で行っても構いませんが、許諾者があなたやあなたの利用行為を支持していると示唆するような方法は除きます。
• 継承 – もしあなたがこの作品をリミックスしたり、改変したり、加工した場合には、あなたはあなたの貢献部分を元の作品とこれと同一または互換性があるライセンスの下に頒布しなければなりません。

https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue

	このファイルは、 クリエイティブ・コモンズ 表示-非営利 3.0 非移植ライセンスのもとで利用許諾されています。
	あなたは以下の条件に従う場合に限り、自由に 共有 – 本作品を複製、頒布、展示、実演することができます。 再構成 – 二次的著作物を作成することができます。 あなたの従うべき条件は以下の通りです。 表示 – あなたは原著作者または許諾者が指定した方法でこの作品のクレジットを表示しなければなりません（ただし、彼らがあなたあるいはこの作品のあなたによる利用を保証していると示唆するような方法はいけません）。 非営利 – あなたはこの作品を営利目的で利用してはなりません。

この文書は、フリーソフトウェア財団発行のGNUフリー文書利用許諾書 (GNU Free Documentation License) 1.2またはそれ以降のバージョンの規約に基づき、複製や再配布、改変が許可されます。不可変更部分、表紙、背表紙はありません。このライセンスの複製は、GNUフリー文書利用許諾書という章に含まれています。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

あなたは上記のライセンスから、どれか一つ以上を選択できます。

# parametric with p&v
set parametric
set dummy p,v

# constant
N = 2.5
lambda = 0.12
gamma = 2.45
M= 1.40

p0 = 200

#function
NCL(p) = N - lambda * log(p)        # v
NCL_p(v) = exp((N - v)/lambda)      # p

CSL1(p) = gamma - lambda * log(p)   # v
CSL2(p) = M * p                     # q
CSL3(v) = M * exp((gamma-v)/lambda) # q

CSL1_p(v) = exp((gamma - v)/lambda) # p

CSL_CU_p(p0) = p0 * exp((gamma - N)/lambda)
CSL_CU_v(p0) = NCL(p0)
CSL_CU_q(p0) = CSL2(CSL_CU_p(p0))

CSL_CD_p(p0) = 1/(1-M/3.0) * p0
CSL_CD_v(p0) = CSL1(CSL_CD_p(p0))
CSL_CD_q(p0) = CSL2(CSL_CD_p(p0))

#other constant
vtics = 0.1
ptics = 100
qtics = 200

max_v = ceil(NCL(100)     /vtics)*vtics
max_p = ceil(CSL_CD_p(300)/ptics)*ptics
max_q = ceil(CSL2(max_p)  /qtics)*qtics

min_v = int(CSL1(max_p) * 0.98 /vtics)*vtics
min_p = 0
min_q = 0

max_x1 = max_p
max_x2 = max_v

min_x1 = -1
min_x2 = min_v

#range
set urange[min_x1:max_x1] # p
set vrange[min_x2:max_x2] # v

set xrange[min_v:max_v] # p
set yrange[min_p:max_p] # v
set zrange[min_q:max_q] # q

#ticslevel
set ticslevel 0

#label
set xlabel "v"
set ylabel "p'"
set zlabel "q"

#grid
set xtics vtics nomirror
set ytics ptics nomirror
set ztics qtics nomirror
set mxtics 2
set mytics 2
set mztics 2
set grid xtics ytics ztics mxtics mytics mztics front lt -1 lw 1, lt 0 lw 1

#border
unset border

# surface settings
set isosample 100
set view 60,120,1,1

# stress path
SP_p(p0, v) = p0 + (CSL_CD_p(p0) - p0) * (v  - NCL(p0))/(CSL_CD_v(p0) - NCL(p0))
SP_v(p0, v) = (v > NCL(p0)) ? NCL(p0) : ((v < CSL_CD_v(p0))? CSL_CD_v(p0) : v)
SP_q(p0, v) = CSL_CD_q(p0) * (v - NCL(p0)) / (CSL_CD_v(p0) - NCL(p0))

# plot
splot  SP_v(p, v), SP_p(p, SP_v(p, v)), SP_q(p, SP_v(p, v)) title "Roscoe" ls 2 lw 0.5, \
           NCL(p), p,       0 title "NCL" ls 1 lw 3, \
          CSL1(p), p,       0 notitle     ls 3 lw 2, \
            min_v, p, CSL2(p) notitle     ls 3 lw 2, \
                v, 0, CSL3(v) notitle     ls 3 lw 2, \
          CSL1(p), p, CSL2(p) title "CSL" ls 3 lw 3

# create gif

theta = 0
set terminal gif notransparent size 1280, 960 font "Times New Roman"

load "roscoe_loop.plt"

set output
set terminal wxt

set view 60,theta,1,1

set output sprintf("roscoe%3d.gif", theta)

replot

theta = theta + 5
if(theta<360) reread