cauchys := proc (a, b, d) local j, poc; global c; if nargs = 4 then poc := 
args[4] else poc := 1 end if; for j from poc to d+poc-1 do c[j] := simplify(
sum(subs(n = i,a)*subs(n = j-i+poc,b),i = poc .. j)); lprint(evaln(c[j]) = c[j
]) end do; c[n] := simplify(sum(subs(n = i,a)*subs(n = n-i+poc,b),i = poc .. n
)); lprint(evaln(c[n]) = c[n]); sum(c[n],n = poc .. infinity) end proc;
dirichs := proc (a, b, d) local j, poc; global c; if nargs = 4 then poc := 
args[4] else poc := 1 end if; c[poc] := subs(n = poc,a)*subs(n = poc,b); 
lprint(evaln(c[poc]) = c[poc]); for j from poc+1 to d+poc-1 do c[j] := 
simplify(sum(subs(n = i,a)*subs(n = j,b)+subs(n = j,a)*subs(n = i,b),i = poc 
.. j-1)+subs(n = j,a)*subs(n = j,b)); lprint(evaln(c[j]) = c[j]) end do; c[n]
:= simplify(sum(subs(n = i,a)*b+a*subs(n = i,b),i = poc .. n-1)+a*b); lprint(
evaln(c[n]) = c[n]); sum(c[n],n = poc .. infinity) end proc;
