store=[0 0];

for x=1:0.005:100
k=2*x^3/(x^2-1);
r=2*x^3/(1+x^2)^2;
store=[store
k  r];
end

store=store(2:end,:);

plot(store(:,1),store(:,2),'LineWidth',2)
 xlabel('k','FontSize', 24)
ylabel('r','FontSize', 24)

axis([0 40 0 0.8])