How do I plot a point at the x intercept on this graph?
1 view (last 30 days)
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
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')
0 Comments
See Also
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!