% Create tau objects.
[x, y] = tau('ChebyshevT', [-pi pi], 101);

% Specify problem, conditions and exact solution (if it is known).
sys = {'cos(x)*diff(y1)-(x^2-3*x^4)*y2+3*diff(y3)=cos(x)^2-cos(exp(x))*(-3*x^4+x^2)-6*sin(2*x)';
       'sin(x)*diff(y1)+diff(y2)-y3=sin(x)*cos(x)-cos(2*x)-sin(exp(x))*exp(x)';
       '-y1+exp(x)*y2+x^3*diff(y3)=cos(exp(x))*exp(x) - sin(x) - 2*x^3*sin(2*x)'};
cond = {'y1(-pi)=0';'y2(0)=cos(1)';'y3(pi)=1'};
es = {'sin(x)';'cos(exp(x))';'cos(2*x)'};

% Solve the problem.
a = tausolver(x, y, ............................................................% Tau variables.
              sys, ..........................................% System of differential equations.
              cond, ...............................................................% Conditions.
              'exact_solution', es, ......................................% Solution (optional).
              'spy', 1); ...............................................% Shows spy of T matrix.

Published with MATLAB® R2016b