
k=1.11/355;
ks=k^2;
s=[0:.01:6]
Rp=4*ks^3*s.^3 + 8*ks^2*s.^3 + 12*ks^2*s.^2 + 4*ks*s.^3 - 20*ks*s.^2 + 12*ks*s - s + 4;

figure(1)

plot(s,Rp,s,zeros(1,length(s)),'--');
