function f=func3(delu1,p)

mtau2=p(1);
Mr2=p(2);
del2=p(3);
Dku1=p(4);
Mrm1=p(5);
Dkm1=p(6);
Dr1x=p(7);
grate=p(8);

Mru1     = (1-mtau2)*(Mr2-del2)+delu1;
Mxu1     = Mru1*Dku1+Mrm1*Dkm1-Dr1x;
f = delu1  - Mxu1/Dku1 +grate;


