(* :Author: Deepak *)
(* :Date: 2021-01-04 *)

Needs["owenScalapinoFreeEnergy`"];

(* Populates the figures directory as ../figures *)
figuresDirectory = FileNameJoin[{
	ParentDirectory[
		DirectoryName[
			FileNameJoin[{
				Directory[],
				$ScriptCommandLine[[1]]
			}]
		]
	], "figures"
}];

figure[filename_] := FileNameJoin[{figuresDirectory, filename}];

label = AxesLabel -> {"\[CapitalDelta]", "Fs(\[CapitalDelta])"};

plot1 = ListPlot[{
	(* Usage for Fs *)
	(* Fs[delta, temp, mustar, debyeFreq, N0, V0] *)
	ParallelTable[{d, Fs[d, .05, 0, 100, 1, .2]}, {d, .05, 2, .1}],
	ParallelTable[{d, Fs[d, .05, .5, 100, 1, .2]}, {d, .05, 2, .1}],
	ParallelTable[{d, Fs[d, .05, 1, 100, 1, .2]}, {d, .05, 2, .1}]
}, ImageSize->Large, PlotLegends->{"\[Mu] = 0", "\[Mu] = .5", "\[Mu] = 1"}, label];
Export[figure["fig1.jpg"], plot1, ImageResolution -> 1200];