l = 1;
g = 35;

x = linspace(0.05,3,1000);

for l = 1:4
    v = l*(l+1)*x.^-2 - g*exp(-x)./x;
    plot(x,v);
    hold on;
end

grid on;
xlabel('x');
ylabel('V_{eff}');
ylim([-150 400]);