clear;clc;close all
[p,c] =deal(8,2);
[mu,sigma] = deal(120,10);
[r,s] = deal(1,0.8);     
omega =@(q) (((1+s)*(r-s)+(2+r+s))*p*(1-arrayfun(@(t)quad(@(x)normpdf(x,mu,sigma),0,t),q))...
                 +(2+s)*(r-s)*p*q.*normpdf(q,mu,sigma)+((1+r)*(r-s)-(2+r+s))*c)./(2+r+s)*(r-s)-2.5;
x0=120;
qq=fzero(omega,x0)

题目在附近里面，请查看，求高手帮忙解答啊。非常感谢了