function [P] = DiffprotDataset(PCTRL,PTRT,a,b,k,m0,t0)
N=size(PCTRL,1);
P=zeros(N,3);
ll1=0;
ll2=0;
ll3=0;
for i=1:N
    [DP1,DP2,DP3]=DiffProtMatrix(PCTRL(1,:),PTRT(i,:),a,b,k,m0,t0);
    P(i,1)=DP1.PO1;
    P(i,2)=DP2.PO1;
    P(i,3)=DP3.PO1;
    ll1=ll1+DP1.ll;
    ll2=ll2+DP2.ll;
    ll3=ll3+DP3.ll;
end
end