clear;close;
data=[459 362 624 542 509 584 433 748 815 505 ...
   612 452 434 982 640 742 565 706 593 680 ...
   926 653 164 487 734 608 428 1153 593 844 ...
   527 552 513 781 474 388 824 538 862 659 ...
   775 859 755 649 697 515 628 954 771 609 ...
   402 960 885 610 292 837 473 677 358 638 ...
   699 634 555 570 84 416 606 1062 484 120 ...
   447 654 564 339 280 246 687 539 790 581 ...
   621 724 531 512 577 496 468 499 544 645 ...
   764 558 378 765 666 763 217 715 310 851 ]';
%hist(data,20);
xbar=mean(data);sigma=std(data);
x=0.1:0.1:1200;
minfee=inf;
for n=100:50:200
   for m=220:20:300
      y=lathe2fun(x,n,m).*normpdf(x,xbar,sigma);
      f=trapz(x,y);
      %f=quad8('lathefun2',0,1400,[],[],n,m);
      if f<minfee
         minfee=f;
         n0=n;m0=m;
      end
   end
end
n0,m0,minfee