ExpReturn =[ -0.0072 0.0026 -0.0054];
ExpCovariance =[0.0251 0.0134 0.0134;
0.0134 0.0239 0.0085;
0.0134 0.0085 0.0178];
portopt(ExpReturn,ExpCovariance)
hold on
Weights=rand(4000,3);
Total=sum(Weights,2);
Total=Total(:,ones(3,1));
Weights=Weights./Total;
[PortRisk,PortReturn]=portstats(ExpReturn,ExpCovariance, Weights);
plot(PortRisk,PortReturn,'.r')