#C10.txt: motivating Gibbs-Sampling Help:=proc(): print(`IterS(f,x,x0,eps,N); IterS((1-x^2)/2,x,1,10^(-10),100)`): print(`MIterS(F,var,var0,eps)`): print(`MIterS([1+3*x+y,2+3*x-y],[x,y],[1/2,1/2],10^(-10),100)`); print(` MIterS([(1+x+y)/(2+x+y), (5+x+y)/(2+3*x+y)],[x,y],[1,1],10^(-7),100);`): print(`GenBiPoisson(L1,L2,n,N);`): print(`GenBiPoisson(1.2,3.4,10,20);`): end: with(Statistics): Digits:=30: #Solving the eq. x=f(x), starting at x=x0 #(y+1)^2=2 , y^2+2y+1=2, y=(1-y^2)/2 IterS:=proc(f,x,x0,eps,N) local x1,x2,co: x1:=x0: x2:=evalf(subs(x=x1,f)): co:=0: while evalf(abs(x2-x1))>eps and coeps and co