Mathematica, 8364 83 90 bytes
i=1;v=Tr[#*x^Range@Length@#]&;While[1<2,Echo@Limit[D[v@#Do[Echo@Limit[D[#/v@#2#2/i!&@@Fold[x#+#2&]/@#,{x,i}],x->0];i++]&>0],{i,∞}]&
Thanks to @ngenisis and @Jenny_mathy !
Take input as two lists.
Need Alt+.
to terminate the execution to see the result. Frontend may crash due to rapid output.
Limit[D[]]
way is 2 83 bytes shorter thanversion SeriesCofficient[]
(@Jenny_mathy):
I remember there is a built-in to replace v
but I cannot find it.
i=1;v=Tr[#*x^Range@Length@#]&;While[1<2,Echo@Limit[D[v@#/v@#2/i!,{x,i}],x->0];i++]&