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

% Solve the problem.
a = tausolver( ...
    x, ..............................................................% Independent tau variable.
    y, ................................................................% Dependent tau variable.
    {['y1-volt(y1,''sin(x-t)-1'')-volt(y2,''1-t*cos(x)'')', ...
    '=-0.5*(x-2)*sin(x)-x*cos(x)^2+(sin(x)+2)*cos(x)-1']; ...
    'y2-volt(y1)-volt(y2,''x-t'')=-x+sin(x)'}, .......................................% Problem.
    {'no'}, .......................................................................% Conditions.
    'exact_solution', {'cos(x)';'sin(x)'}); ...................................% Exact solution.

% Data to compare: from L. H. Yang (2013) and M. Rabbani (2007).
Data = [ 0.0 0 6.93348e-11 0 3.60316e-11
0.1 1.37735e-4 4.53518e-09 1.52721e-4 2.75123e-08
0.2 9.27188e-4 8.84879e-09 1.14715e-3 3.10611e-08
0.3 2.67117e-3 1.28253e-08 3.71248e-3 3.53307e-08
0.4 5.45507e-3 1.65442e-08 8.57201e-3 4.03402e-08
0.5 9.22670e-3 2.00881e-08 1.64412e-2 4.61209e-08
0.6 1.38644e-2 2.35657e-09 2.78243e-2 5.27214e-08
0.7 1.92960e-2 2.71160e-08 4.25337e-2 6.02041e-08
0.8 2.56349e-2 3.09302e-08 5.91212e-2 6.86601e-08
0.9 3.31574e-2 3.52645e-08 7.48883e-2 7.82029e-08
1.0 4.19808e-2 3.67322e-08 8.70896e-2 1.02387e-07];

figure(2), hold on;
plot(Data(:, 1), Data(:, 2), '-^', 'Linewidth', 1.6)
plot(Data(:, 1), Data(:, 4), '-^', 'Linewidth', 1.6)
plot(Data(:, 1), Data(:, 3), '-s', 'Linewidth', 1.6)
plot(Data(:, 1), Data(:, 5), '-s', 'Linewidth', 1.6)
legend('y1 Tau Toolbox','y2 Tau Toolbox', ...
       'y1 Rabbani (2007)','y2 Rabbani (2007)', ...
       'y1 Yang (2013)','y2 Yang (2013)')