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

% Solve the problem.
a = tausolver(x, y, .............................................................% Tau variables.
              'x*diff(y, 2)+2*diff(y)-(6*x+4*x^3)*y=0' , .......% Ordinary differential equation.
              {'y(0)=1';'y''(0)=0'}, ...............................................% Conditions.
              'exact_solution', 'exp(x^2)'); ..............................% Solution (optional).