3.3积分问题的解析解

>> syms x;
y=sin(x)/(x^2+4*x+3);
y1=diff(y),y0=int(y1)
 
y1 =
 
cos(x)/(x^2 + 4*x + 3) - (sin(x)*(2*x + 4))/(x^2 + 4*x + 3)^2
 
 
y0 =
 
sin(x)/(x^2 + 4*x + 3)
 
>> y4=diff(y,4),y0=int(int(int(int(y4)))),simplify(y0)
 
y4 =
 
sin(x)/(x^2 + 4*x + 3) + (12*sin(x))/(x^2 + 4*x + 3)^2 + (24*sin(x))/(x^2 + 4*x + 3)^3 - (24*cos(x)*(2*x + 4)^3)/(x^2 + 4*x + 3)^4 - (12*sin(x)*(2*x + 4)^2)/(x^2 + 4*x + 3)^3 - (48*sin(x)*(2*x + 4)^2)/(x^2 + 4*x + 3)^4 + (24*sin(x)*(2*x + 4)^4)/(x^2 + 4*x + 3)^5 + (4*cos(x)*(2*x + 4))/(x^2 + 4*x + 3)^2 + (16*cos(x)*(2*x + 4))/(x^2 + 4*x + 3)^3 + (8*cos(x)*(8*x + 16))/(x^2 + 4*x + 3)^3 - (6*sin(x)*(2*x + 4)*(8*x + 16))/(x^2 + 4*x + 3)^4
 
 
y0 =
 
sin(x)/(x^2 + 4*x + 3)
 
 
ans =
 
sin(x)/(x^2 + 4*x + 3)

>> syms x;
>> int(exp(-x^2/2))
 
ans =
 
(2^(1/2)*pi^(1/2)*erf((2^(1/2)*x)/2))/2

3.3.2

>> syms x;
I1=int(exp(-x^2/2),x,0,1.5),vpa(I1),I2=int(exp(-x^2/2),x,0,inf)
 
I1 =
 
(2^(1/2)*pi^(1/2)*erf((3*2^(1/2))/4))/2
 
 
ans =
 
1.0858533176660165697024190765423
 
 
I2 =
 
(2^(1/2)*pi^(1/2))/2

>> syms x t ;
f(x)=(-2*x^2+1)/(2*x^2-3*x+1)^2;
I1=int(f)
 
I1(x) =
 
1/(x - 1) - 1/(2*x - 1)
 
>> I=I1(exp(-2*t))-I1(cos(t))
 
I =
 
1/(2*cos(t) - 1) - 1/(cos(t) - 1) - 1/(2*exp(-2*t) - 1) + 1/(exp(-2*t) - 1)

>> syms x;
>> f=1/x/sqrt(1-log(x)^2);
>> I=int(f,x,1,2*exp(sym(1)))
 
I =
 
asin(log(2) + 1)

>> syms x y z;

>> f0=-4*z*exp(- y*x^2 - z^2)*(cos(x^2*y) - sin(x^2*y) + 4*x^4*y^2*cos(x^2*y) + 4*x^4*y^2*sin(x^2*y) - 10*x^2*y*cos(x^2*y));

>> f1=int(f0,z),f1=int(f1,y),f1=int(f0,x),f1=int(f1,x)
 
f1 =
 
exp(- y*x^2 - z^2)*(2*cos(x^2*y) - 2*sin(x^2*y) + 8*x^4*y^2*cos(x^2*y) + 8*x^4*y^2*sin(x^2*y) - 20*x^2*y*cos(x^2*y))
 
 
f1 =
 
-2*y*exp(- y*x^2 - z^2)*(sin(x^2*y) - cos(x^2*y) + 4*x^2*y*cos(x^2*y))
 
 
f1 =
 
4*x*z*exp(- y*x^2 - z^2)*(sin(x^2*y) - cos(x^2*y) + 2*x^2*y*cos(x^2*y))
 
 
f1 =
 
-2*x^2*z*exp(- y*x^2 - z^2)*(cos(x^2*y) - sin(x^2*y))
 
 
 
>> f2=int(f0,z),f2=int(f2,x),f2=int(f2,x),f2=int(f2,y)
 
f2 =
 
exp(- y*x^2 - z^2)*(2*cos(x^2*y) - 2*sin(x^2*y) + 8*x^4*y^2*cos(x^2*y) + 8*x^4*y^2*sin(x^2*y) - 20*x^2*y*cos(x^2*y))
 
 
f2 =
 
-2*x*exp(- y*x^2 - z^2)*(sin(x^2*y) - cos(x^2*y) + 2*x^2*y*cos(x^2*y))
 
 
f2 =
 
x^2*exp(- y*x^2 - z^2)*(cos(x^2*y) - sin(x^2*y))
 
 
f2 =
 
sin(x^2*y)*exp(- y*x^2 - z^2)

 
>> syms x y z; 
int (int(int(4*x*z*exp(-x^2*y-z^2),x,0,2),y,0,pi),z,0,pi)
 
ans =
 
-(exp(-pi^2) - 1)*(eulergamma + log(4*pi) - ei(-4*pi))
 
>> vpa(ans)
 
ans =
 
3.1080794020854127228346146476714