RL=[1:1:20]';
i=-8./(12+RL);
p=i.^2.*RL;
plot(RL,p);