feval(@assignin,'caller','score',5000); % msec

a=[100 99 100 10 9 2 2 10 30 70]; d=[100 100 99 9 10 1 2 2 30 80]; y_c=[0.8079031789315619 0.7888693135658454 0.8230449788340404 0.5580697529719042 0.3798720048109818 0.4166666666666667 0.10609567901234569 0.9901146432872121 0.633266311153744 0.5011352886279803]; tsum=0; for i=1:length(a) ta=clock; y=risk_prob(a(i), d(i)); t1=etime(clock,ta)*1000; % time in msec tsum=tsum+t1; assert(abs(y - y_c(i)) <= 1e-6,sprintf('A=%i D=%i Expect=%.9f pwin=%.9f',a(i),d(i),y_c(i),y)) fprintf('A %3i D %3i Time(msec) %7.3f\n',a(i),d(i),t1); end feval( @assignin,'caller','score',floor(min( 5000,tsum )) );