[Phi] = [Psi] E^(-(x^2/(2 [Delta]^2)))

 netcharDe = D[D[[Phi], x], x];

passDens = Integrate[
distfuncMax, {v, -Infinity, -Sqrt[[Phi][[Psi], [Delta],
    x]]}] + Integrate[
distfuncMax, {v, Sqrt[[Phi][[Psi], [Delta], x]], Infinity}];

trapcharDens = netcharDe - passDens + 1;