function out=peaplants(in)

pYY=in(1);
pYg=in(2);
pgg=in(3);


qYY=pYY.^2+pYY.*pYg+pYg.^2/4;
qYg=pYY.*pYg+pYY.*pgg+pYg.^2/2+pYg.*pgg+pYY.*pgg;
qgg=pYg.^2/4+pYg.*pgg+pgg.^2;

out(1)=qYY;
out(2)=qYg;
out(3)=qgg;