<<TwoDiff.m

TestTwoDiff[level_Integer] :=
	Block[{a, b, a1, a2, a3, b1, b2, b3, 
		cputime, zeta = N[0.5-0.5*I,16], diff},
	a = {{1001, 999}, {999, 1001}};
	b = {{I*1001, 999}, {-999, I*1001}};
	a1 = 1; a2 = -1; a3 = -I;
	b1 = (1001*I+999)/(-999+1001*I);
	b2 = (-1001*I+999)/(999+1001*I);
	b3 = (-1001*I+999)/(-999*I+1001);
	cputime = TimeUsed[];
	diff = TwoDiff[{a1, a2, a3}, {b1, b2, b3}, zeta, level, {a,b}];
	cputime = TimeUsed[] - cputime;
	{cputime, diff}]

