我试图通过在 Matlab 中绘制给定的势函数来找到三体系统的 5 个拉格朗日点。唯一的问题是我不太擅长编程。任何帮助将不胜感激。我想知道的是如何让这段代码给我一个像样的等高线图:
function Lagrange(a)
x = ( -10000: 1 : 10000);
y = ( -10000: 1 : 10000);
Potential = zeros(length(x));
for i = 1: length(x)
for j = 1 : length(y)
Potential(i,j) = ( 1 - a ) / sqrt( ( x(i) - a )^2 + y(j)^2) + a / sqrt( ( x(i) + 1 - a )^2 + y(j)^2 ) + ( x(i)^2 + y(j)^2 ) / 2 ;
end
j = 1;
end
contour(Potential);
xlabel('X axis');
ylabel('Y axis');
zlabel('Z axis');