Wo=10, T=500, p=100, f(x) follows the normal distribution N (350,100);
This is a function image, obviously an increase function, but it will be reduced to 0.
This is the code
p=100;
wo=10;
T=500;
f1=@(x,mu)(x.*p.*normpdf(x,350,10));
Ewr=@(mu) quadgk(@(x) f1(x,mu),0,(mu-1/wo)*T)
B=0:1:1000;
Ewr_mu=arrayfun(@(mu) Ewr(mu),B);
plot(B,Ewr_mu)
grid on
title('Image of expected profit')
xlabel('mu')
ylabel('E(wr)')
Anyone can help me??Thanks!!!
Best Answer