% Create tau objects.
[x, y] = tau('ChebyshevT', [-1 1], 6);

% Solve the problem.
a = tausolver(x, y, ............................................................% Tau variables.
            {'y2-1*diff(y1)=0'; ...................................................% Equation 1.
             'y3-1*diff(y2)=0'; ...................................................% Equation 2.
             '(x^2+1)*diff(y3)-(x^2+3*x)*diff(y2)+5*x*y2-5*y1=60*x^2-10'}, ........% Equation 3.
             {'y1(-1)=4';'y2(1)=2';'y3(0)=0'}, ....................................% Conditions.
             'exact_solution', {'x^5-3*x+2';'5*x^4-3';'20*x^3'}, .........% Solution (optional).
             'spy', 1); ................................................% Shows spy of T matrix.