> dsolve({diff(y(t), t) = y(t)^3/(t + 1), y(0) = 1}, y(t));
                                  1           

                 y(t) = ----------------------

                                         (1/2)

                        (1 - 2 ln(t + 1))     


;
> dsolve({D(D(y))(t) - 3*D(y)(t) + 2*y(t) = 0, y(0) = 2, D(y)(0) = 3}, y(t));
                    y(t) = exp(2 t) + exp(t)


;
> with(LinearAlgebra);
> A := Matrix([[3, -4], [4, 3]]);
> evalf(Eigenvalues(A));
> evalf(Eigenvectors(A))[2];
                           [3  &uminus0;4]

                      A := [             ]

                           [4      3     ]


                     [3. + 4. ⅈ]

                     [                    ]

                     [3. - 4. ⅈ]


             [ⅈ  &uminus0;ⅈ]

             [                                   ]

             [     1.                1.          ]


;
> 
;