利用matrix.h解一些高数考题

EquivInf(upoly("x")*Integral(Sin(Sin("x")*Sin("x"))*Deriv(Sin("x"))))

1/3x^4

#include "poly.h"
#define Num poly
#include "euclid.h"
using namespace std;
int main(){
    Matrix P=poly("t")*Plane(1,-3,-1,-4)+Plane(2,-1,1,-1);
    poly eq=dist2_dot_pl(Vector(0,0,0),P)-1;
    cout<<Factorization(eq)<<endl;
}

Limit(2*Sin("x")*Ln("1+x")-2*upoly("x")*Arctan("x")+upoly("x^3"),upoly("x^4"))

EquivInf(upoly("x")-Sin("x")+Integral(upoly("x^2")*Exp("x^2")))

1/2x^3

Limit(Exp("1/2x")-Pow("1-x",frac(-1,2)),upoly("x")*Tan("3x"))

-1/12

2x^2
2
x-1 2 
x^2+1 1

0
1ln|x-1|
-1(x-1)^-1
-1/2ln(x^2+1)+(0)/\4arctan((2x+0)/\4)
C

EquivInf(Integral(Sin("x^2")))

1/3x^3

diff ans=Deriv(Deriv(d))-2*Deriv(d)+d;
Matrix f=Matrix('C',to_v("2x"));

diff ans=Deriv(d);
Matrix f=Matrix('C',to_v("x^2-1"));

EquivInf(variableLimit("3",Exp("x"),Exp("2x")));

C0 * (1)e^-1x
(-5/4x)e^-1x

[(1/10)cos(1x)+(-1/5)sin(1x)]e^0x

#include "poly.h"
#define Num poly
#include "euclid.h"
using namespace std;
int main(){
    Line L1=make_pair(Vector(0,0,-1),Vector(1,-1,0));
    Line L2=make_pair(Vector(0,0,1),Vector(1,1,0));
    output(distLine(L1,L2));
}