@Kitonum I am looking for about this for a long time

^^ It is really nice bro

@acer 

zero_one := proc(A, x, n, widthgoal)
local Astar, var, L, Aprime, Lprime, Aprime2, Lprime2, y;
1 Astar := expand(subs(y = x+1,expand(y^n*subs(x = 1/y,A))));
2 var := polyvariations(Astar,x,n);
3 if var = 0 then
4 return []
elif var = 1 then
5 if 1 <= widthgoal then
6 return [[0, 1]]
else
7 return midpoint(A,x,n,widthgoal)
end if
end if;
8 Aprime := subs(x = 1/2*x,A);
9 if subs(x = 1,Aprime) = 0 then
10 L := [[1/2, 1/2]];
11 divide(A,2*x-1,'Aprime');
12 Aprime := subs(x = 1/2*x,Aprime)
else
13 L := []
end if;
14 Aprime := expand(2^n*Aprime);
15 Lprime := zero_one(Aprime,x,n,2*widthgoal);
16 L := [op(L), op(map((x, y) -> map(y,x),Lprime,z -> 1/2*z))];
17 Aprime2 := expand(subs(x = x+1,Aprime));
18 Lprime2 := zero_one(Aprime2,x,n,2*widthgoal);
19 L := [op(L), op(map((x, y) -> map(y,x),Lprime2,z -> 1/2*z+1/2))]
end proc

@acer It is another Proc and doesnt affected here. Hope you can help me fix my code.

Thank a lot.

It is still running good in Maple 16 :(