MATLAB: How to solve an equation with variable in integral interval

integralvariable interval values

Hi, I am trying to solve a system of equations with variable in integral interval.
fun = @(x) x*normpdf(x);
f1: (1-normcdf(q,3,1))*p/(1+r) == c*(1-normcdf((c*q-equity)*(1+r)/p,3,1))
f2: (1-normcdf((c*q-equity)*(1+r)/p))*(1+r)*(c*q-equity)+p*integral(fun,0,(c*q-equity)*(1+r)/p) == (c*q-equity)
q and r are the variables and p, c, and equity are parameters. When I tried to solve this with vpasolve, it could not solve this probably because a variable is also in the integral interval.
Can anybody please help me with this problem? Thanks in advance.

Best Answer