How do I plot a point at the x intercept on this graph?
Show older comments
e=1;
sigma=0.154;
r=linspace(0.154,1,100);
U = @(e,sigma,r) 4*e*((sigma./r).^12)-((sigma./r).^6);
norm_r = r./sigma;
figure
plot(norm_r,U(e,sigma,r))
title('Lennard-Jones Potential')
xlabel('Distance (nm)')
ylabel('Potential Energy')
2 Comments
dpb
on 20 Feb 2018
Not sure precisely what you want to add where???
Rhianne Tallarico
on 20 Feb 2018
Answers (1)
Are Mjaavatten
on 20 Feb 2018
Let U be a function of r only and find the solution to U(r) = 0:
e=1;
sigma=0.154;
r=linspace(0.154,1,100);
U = @(r) 4*e*((sigma./r).^12)-((sigma./r).^6);
norm_r = r./sigma;
figure
plot(norm_r,U(r))
title('Lennard-Jones Potential')
xlabel('Distance (nm)')
ylabel('Potential Energy')
norm_r_0 = fzero(U,1)/sigma;
hold on; plot(norm_r_0,0,'*r')
Categories
Find more on Quantum Mechanics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!