maple解常微分方程

factor(x^4+2*x^3-13*x^2-14*x+24);
                (x - 1) (x - 3) (x + 4) (x + 2)
expand((x^2-1)*(y-1)^3);
       2  3      2  2      2      3    2      2          
      x  y  - 3 x  y  + 3 x  y - y  - x  + 3 y  - 3 y + 1
gcd(133, 627);
                               19
lcm(21, 35, 45);
                              315
add(2*k-1, k = 1 .. 111);
                             12321
sum(k^3, k = 1 .. n);
              1        4   1        3   1        2
              - (n + 1)  - - (n + 1)  + - (n + 1) 
              4            2            4         
add(k^3, k = 1 .. n);
Error, unable to execute add
add(k^3, k = 1 .. 5);
                              225
solve({x-y = 5, 2*x+y = 1}, {x, y});
                        {x = 2, y = -3}
diff(x/sin(x), x$4);
           3                          4              2         
  24 cos(x)    20 cos(x)   24 x cos(x)    28 x cos(x)     5 x  
- ---------- - --------- + ------------ + ------------ + ------
         4            2            5              3      sin(x)
   sin(x)       sin(x)       sin(x)         sin(x)             
diff(x, x$1);
                               1
diff(x^2, x$1);
                              2 x
diff(x^2, x$2);
                               2

 maple求常微分方程（根据剪胀方程和相关联流动法则，求剑桥模型其势函数）：

f := (p-(1/2)*px)^2/((1/2)*px)^2+q^2/((1/2)*M*px)^2-1;
                               2             
                     /    1   \              
                   4 |p - - px|        2     
                     \    2   /     4 q      
                   ------------- + ------ - 1
                          2         2   2    
                        px         M  px     
simplify(f);
                                   2      2
                     4 p (p - px) M  + 4 q 
                     ----------------------
                              2   2        
                             M  px         

 maple求常微分方程（根据剪胀方程和相关联流动法则，求修正剑桥模型其势函数）：

eqqq := (D(y))(x)+((M*M)*x*x-y(x)*y(x))/(2*x*y(x)) = 0;
                             2  2       2    
                            M  x  - y(x)     
                  D(y)(x) + ------------- = 0
                              2 x y(x)       
dsolve(eqqq);
                         (1/2)                          (1/2)
         /  2  2        \               /  2  2        \     
  y(x) = \-M  x  + _C1 x/     , y(x) = -\-M  x  + _C1 x/