Pour résoudre symboliquement une équation, ou un système d'équation on utilise la fonction solve(eq,var) qui résoud l'équation eq par rapport à la variable var.
>x^3-1/2*x^2*a+13/3*x^2 = 13/6*x*a+10/3*x-5/3*a: eq:=%:
>solve(eq,x);
Si on ne peut pas calculer de solution annalytique, on peut rechercher les racines numériquement avec fsolve.
Pour des systèmes linéaires, on utilise la fonction lsolve.
Pour des équations différentielles, on utilise la fonction dsolve.
>>diff(y(t),t$2)+5*diff(y(t),t)+6*y(t)=0; eqd:=%:
>>dsolve(eqd,y(t));
les constantes d'intégrations et sont déterminées par les conditions initiales.