residue1.m

薛定宇高等应用数学问题的MATLAB求解例子源码

function f=residue1(F,s)
f=sym(0); if nargin==1, syms s; end
[num,den]=numden(F); 
x=solve(den); [x0,ii]=sort(double(x)); x=[x(ii);rand(1)];
kvec=find(diff(double(x))~=0); ee=x(kvec);
kvec=[kvec(1); diff(kvec(:,1))];
for i=1:length(kvec), 
    for j=1:kvec(i), m=kvec(i); s0=ee(i);
        k=limit(diff(F*(s-s0)^m,s,j-1),s,s0);
        f=f+k/(s-s0)^(m-j+1)/factorial(j-1);
    end
end