Как я могу построить потенциал Леннард Джонса, как этот, с gnuplot?

Question

Я начинаю использовать gnuplot, и я пытаюсь построить потенциал LJ как этот, который у меня есть, но только с 12-6. Также я хотел бы показать на графике значения сигмы и эпсильона.

сигма = 2,6

epsillon = 5

Это, как я хотел бы выглядеть следующим образом:

Это код, который я сделал.

set xlabel "Distance r_0 ({\305})" 
set ylabel "Energy E (J/mol) 
f(r,sigma,epsilon)=4*epsilon*((sigma/r)**12-(sigma/r)**6) 
plot[2:5.2] [-6:12] f(x,2.6,5) title '12-6 LJ' 

Буду признателен за любую помощь и совет. Я студент-самоучка и извиняюсь, если этот вопрос тривиален.

Большое спасибо заранее

Answer 1

Вот первая попытка:

f(r,sigma,epsilon)=4*epsilon*((sigma/r)**12-(sigma/r)**6) 

sigma=2.6 
epsilon=5 

set arrow 1 from first 3,0 to first 3, f(3,sigma,epsilon) 
set arrow 2 from first 3,0 to first 2.8, f(2.8,sigma,epsilon) 
set label 1 "(r_0, {/Symbol e}_0)" at first 3,0 center offset 0,1 

set xlabel "Distance r_0 ({\305})" 
set ylabel "Energy E (J/mol) 

set xr [2:5.2] 
set yr [-6:12] 

set key reverse Left at 4,5 

plot for [epsi= epsilon-1:epsilon+1:1 ] f(x,sigma,epsi) title sprintf("12-6 LJ {/Symbol s}=%.1f {/Symbol e}=%.1f",sigma,epsi) 

set term pngcairo enhanced 
set out "LJ.png" 
replot 

и вот результат: