Hi all,

I have the following function:

phiS := (C,q) -> w*(int(a*f*g(x), x=0..a*f) + int(x*g(x), x=a*f..((w+pi2)*q-pi2*C)/w) + int(q*g(x), x=((w+pi2)*q-pi2*C)/w..infinity)) - s*q - pi1*(f-C) - pi2*int((C-q)*g(x), x=((w+pi2)*q-pi2*C)/w..infinity);

when I type the command:

simplify(diff(phiS(C,q),C));

I get the correct output (I've differentiated phiS on paper): pi1- pi2*int(g(x), x=((w+pi2)*q-pi2*C)/w..infinity);

***************