【matlab基础知识代码】（八）不定积分的求解

>> syms x; y=sin(x)/(x^2+4*x+3); y1=diff(y); y0=int(y1)
 
y0 =
 
sin(x)/(x^2 + 4*x + 3)

>>  y4=diff(y,4); y0=int(int(int(int(y4)))); simplify(y0)
 
ans =
 
sin(x)/(x^2 + 4*x + 3)

 

 

>> syms a x;f=int(x^3*cos(a*x)^2);simplify(f)
 
ans =
 
((3*sin(a*x)^2)/8 + (a^3*x^3*sin(2*a*x))/4 - (3*a^2*x^2*(2*sin(a*x)^2 - 1))/8 - (3*a*x*sin(2*a*x))/8)/a^4 + x^4/8

>>  f1=x^4/8+(x^3/(4*a)-3*x/(8*a^3))*sin(2*a*x)+(3*x^2/(8*a^2)-3/(16*a^4))*cos(2*a*x); simplify(f-f1) 
 
ans =
 
3/(16*a^4)

 

 

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

>> syms a x; int(x*sin(a*x^4)*exp(x^2/2))
 
ans =
 
int(x*sin(a*x^4)*exp(x^2/2), x)

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

 

>> 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(-x^2*y-z^2)*(cos(x^2*y)-10*cos(x^2*y)*y*x^2+4*sin(x^2*y)*x^4*y^2+4*cos(x^2*y)*x^4*y^2-sin(x^2*y)); f1=int(f0,z); f1=int(f1,y); f1=int(f1,x); f1=simplify(int(f1,x))
 
f1 =
 
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),z,0,pi),y,0,pi),x,0,2)
 
ans =
 
-(2*exp(-pi^2) - 2)*(eulergamma/2 + log(2) + log(pi)/2 - ei(-4*pi)/2)

 >> vpa(ans)
 
ans =
 
3.1080794020854127228346146476714