{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart():with(stats );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "normalize:=x->(x*10); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "normalize(0.763);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "N:=10:unidata:=map(normalize ,[random[uniform](N)]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 " mean:=l->sum(l[ii],ii=1..nops(l))/nops(l);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 21 "evalf(mean(unidata));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "sqr:=x->x^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "sdist:=l->sqrt(mean(map(sqr,l))-mean(l)^2);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 15 "sdist(unidata);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 31 "A:=5;sigma:=evalf(10/sqrt(12));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "factor(int(x,x=a..b)/(b-a)); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "factor(int((x-(b+a)/2)^ 2,x=a..b)/(b-a));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 150 "M:=10 00:N:=10: p1:=proc(newdata) local k,unidata: for k from 1 to M do unid ata:=map(normalize,[random[uniform](N)]):newdata[k]:=mean(unidata):od: end;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "newdata:='newdata' :p1(newdata);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "newdata:=[seq(newdata[i],i=1 ..M)];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 126 "mycount:=proc(l, a,b) local i,k; k:=0: for i from 1 to nops(l) do if (a " 0 "" {MPLTEXT 1 0 23 "mycount(newdata,5,5.2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "h:=0.025;a:=0;b:=10;Range:=(b-a)/h;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "L := [seq(mycount(newdata,a+i*h,a+(i+1)*h),i=0..Range)];" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "plot([seq([a+i*h,mycount(new data,a+i*h,a+(i+1)*h)/M/h],i=0..Range)]);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "plot(exp(-(x- A)^2/(2*(sigma/sqrt(N))^2))/sqrt(2*Pi),x=0..10);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 " " }}}}{MARK "26" 0 }{VIEWOPTS 1 1 0 1 1 1803 }