matlab code for differential quadrature method
-1) / (N-1)); % Chebyshev points x = (a + b)/2 + (b - a)/2 x; % Scale to domain % Compute weights w1 = computeWeights(x); w2 = w1 w1; % Define source term Q = @(x) -pi^2 sin(pi x); % Setup RHS rhs = Q(x)'; % Boundary conditions T_a = 0; T_b = 0; % Modify system for boundary cond