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

% Specify ode and conditions.
ode = '1e-7*diff(y, 2)-2*x*(cos(x)-0.8)*diff(y)+(cos(x)-0.8)*y=0';
conditions = {'y(-1)=1';'y(1)=1'};

% Solve the problem.
a = tausolver(x, y, ............................................................% Tau variables.
              ode, ............................................% Ordinary differential equation.
              conditions, .........................................................% Conditions.
              'pieces', 50); ................................................% Number of pieces.