% Create tau objects.
[x, y] = tau('LegendreP', [0 10*pi], 100);

% Solve the problem.
a = tausolver(x, y, ............................................................% Tau variables.
              'y + diff(y, 2)=0', .............................% Ordinary differential equation.
              {'y(0)=1';'y''(10*pi)=0'}, ..........................................% Conditions.
              'exact_solution', 'cos(x)', ................................% Solution (optional).
              'pieces', 5, ..................................................% Number of pieces.
              'step', 0.1); .......................................% Step for graphical showing.