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 Annotations 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!