syms z;
int(sin(z),z)
int(z*cos(z^2),z)

2.sym z
f=z/(z+1);
r = taylor(f, z, 1, 'Order', 8)

[R,P,K]=residue([1,-3,2],[1,-7,12])

[R,P,K]=residue([1,0],[1,0,0,0,-1])
S=2pii*sum(R)