clear all
load pippdata.dat
x1=cos(pi*pippdata(1:19,1)/180);
s1=pippdata(1:19,2);
x2=cos(pi*pippdata(20:46,1)/180);
s2=pippdata(20:46,2);
xf=-1:0.1:1;
sf=min(s2)*(1+3*xf.^2);
plot(xf,sf,'-.',x1,s1,'-',x2,s2,'ok','MarkerFaceColor','black','MarkerSize',6.5)
axis([-1 1 0 31])
[legh,objh]=legend('1+3cos^2\theta','Full Fit','Data',0);
set(gca,'Fontsize',15')
xlabel('cos\theta','FontSize',20)
ylabel('d\sigma/d\Omega (mb/sr)','FontSize',20)
print -deps pippelas.eps