f:=(x,mu,sigma) -> 1/sqrt(2*Pi)/sigma*exp(-1/2*((x-mu)/sigma)^2);

NiM+SSJmRzYiZio2JUkieEdGJUkjbXVHRiVJJnNpZ21hR0YlRiU2JEkpb3BlcmF0b3JHRiVJJmFycm93R0YlRiUqKC1JJXNxcnRHRiU2IywkSSNQaUdJKnByb3RlY3RlZEdGNCIiIyEiIjkmRjYtSSRleHBHRiU2IywkKiYsJjkkIiIiOSVGNkY1RjchIiMjRjZGNUY/RiVGJUYl

assume(sigma>0);ut:=int(-exp(-lambda*(q*(x-p)+cash))*f(x,mu,sigma),x=-infinity..infinity);

NiM+SSN1dEc2IiwkLUkkZXhwRzYkSSpwcm90ZWN0ZWRHRipJKF9zeXNsaWJHRiU2IywkKiZJJ2xhbWJkYUdGJSIiIiwqKiZJInFHRiVGMEkicEdGJUYwIiIjSSVjYXNoR0YlISIjKihGL0YwRjNGNUknc2lnbWF8aXJHRiVGNUYwKiZGM0YwSSNtdUdGJUYwRjdGMCNGMEY1ISIi

diff(ut,q);

NiMsJCooSSdsYW1iZGFHNiIiIiIsKEkicEdGJiIiIyooRiVGJ0kicUdGJkYnSSdzaWdtYXxpckdGJkYqRipJI211R0YmISIjRictSSRleHBHNiRJKnByb3RlY3RlZEdGM0koX3N5c2xpYkdGJjYjLCQqJkYlRicsKiomRixGJ0YpRidGKkklY2FzaEdGJkYvKihGJUYnRixGKkYtRipGJyomRixGJ0YuRidGL0YnI0YnRipGJyMhIiJGKg==

solve(%,q);

NiMqKCwmSSJwRzYiISIiSSNtdUdGJiIiIkYpSSdsYW1iZGFHRiZGJ0knc2lnbWF8aXJHRiYhIiM=