%input
syms a x;
f = exp(-5*x) * sin(3*x + pi/3);
T = taylor(f,x,'Order', 5); pretty(T) % 此处展开到前 5 项
T2 = taylor(f,x,a,'Order', 5); pretty(T2) % 此处展开到前 5 项
%input
%第一种方式
syms n k;
s = symsum(1/(1+k*pi/n^2),k,1,n) % 先求和
limit(s/n,n,Inf) % 再求极限
%第二种方式
syms k n
limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)
%output
s =
(n^2*psi(1 + n + n^2/pi))/pi - (n^2*psi(1 + n^2/pi))/pi
ans =
1
x=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2];
y=[0,2.208,3.206,3.444,3.241,2.816,2.311,1.81, 1.36,0.982,0.679,0.447,0.277];
[dy1,dx1]=diff_ctr(y,x(2)-x(1),1);
[dy2,dx2]=diff_ctr(y,x(2)-x(1),2);
[dy3,dx3]=diff_ctr(y,x(2)-x(1),3);
[dy4,dx4]=diff_ctr(y,x(2)-x(1),4);
figure,subplot(221),plot(dx1+x(1),dy1,'r-');subplot(222),plot(dx2+x(1),dy2,'g-');
subplot(223),plot(dx3+x(1),dy3,'b-');subplot(224),plot(dx4+x(1),dy4,'k-');
x = [0:0.1:1.2]';
y = [0,2.208,3.206,3.444,3.241,2.816,2.311,1.81,1.36,0.982,0.679,0.447,0.277]';
S = trapz(x,y)
输出
S =
2.2642
format long
f = @(x)exp(x).*log(x);
q = quadgk(f,0,1)
输出
q =
-1.317902162414081
f = @(x,y,z) y.*sin(x) + z.*cos(x);
q1 = integral3(f,0,pi,0,1,-1,1)
% q1 = integral3(f,0,pi,0,1,-1,1,'AbsTol', 0,'RelTol',1e-9)
输出:
q1 =
2.0000