matlabFunction をご使用ください
syms e w t
A=5;m=1;s=1;wA=1;c=10^-10;
%分母の関数を因数分解
eqn=w.^2+(-m-wA+(sqrt(2).*s+c)*1i).*w+(m-(sqrt(2).*s+c).*1i).*wA-sqrt(2).*A.*s==0;
sols=simplify(solve(eqn,w));
%部分分数分解
F=partfrac((w-m+(sqrt(2)+c)*1i)./((w-sols(1)).*(w-sols(2))),w,"FactorMode","full");
S1(t)=(1i./2*pi).*int(F*exp(-1i.*w.*t),w,-100,100);
T=0:0.10:1;
% tic
% S2=subs(S1,t,T);
% toc
F = matlabFunction(S1)
S3=zeros(length(T),1);
tic
for ii = 1:length(T)
S3(ii)= F(T(ii));
end
toc
S3
