equation = @(x,y) 1e-7*diff(y, 2)-2*x*(cos(x)-0.8)*diff(y)+(cos(x)-0.8)*y;
domain = [-1 1];
conditions = @(y) {y(-1) - 1; y(1) - 1};
options = tau.settings('degree', 100, 'pieces', 50, 'basis', 'LegendreP');
problem = tau.problem(equation, domain, conditions, options);
% solution via tau method
[yn, info] = tau.solve(problem);
plot(yn);
