thePSolution=DSolve[{p11'[t]==2 p12[t],p12'[t]==p22[t]-b p12[t],p22'[t]==-2 b p22[t]+(s^2),p11[0]==0,p12[0]==0,p22[0]==(s^2)},{p11,p12,p22},{t}]