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

% Solve the problem.
a = tausolver(x, y, ............................................................% Tau variables.
    {['diff(y1) + fred(y1, ''exp(x+t)'')+x*y2 ', ...
    '-(x+x^2)*volt(y2, ''sin(x-t)'')', ...
    '= cos(x)+exp(x)+x*(cos(x)-exp(x))-', ...
    '(x^2+x)*(cos(x)/2-exp(x)/2+sin(x)/2+(x*sin(x))/2)+', ...
    '(exp(1)*exp(x)*(exp(1)-cos(1)+sin(1)))/2']; ...
    ['x^2*y1 - x*volt(y1, ''2*x+sin(t)'')-3*diff(y2)' ...
    '-fred(y2, ''x-t'')', ...
    '= cos(1)+sin(1)+3*exp(x)+', ...
    '3*sin(x)-x*(x/2-2*x*(cos(x)-1)+2*x*(exp(x)-1)-', ...
    '(cos(x)*sin(x))/2-(exp(x)*(cos(x)-sin(x)))/2+1/2)-', ...
    'x*sin(1)+x^2*(exp(x)+sin(x))+x*(exp(1)-1)-2']}, .................................% Problem.
    {'y1(0)=1';'y2(0)=0'}); .......................................................% Conditions.