30 Apr 2009
05 Apr 2012)
Numerical computation with functions instead of numbers.
function pass = chebop_linop
% Test the chebop linop method.
% Basically test that the direct construction with diff and diag is equivalent
% to constructing a chebop through an anonymous function.
% Asgeir Birkisson, December 2010
% Create a linop, u -> u''+x*u, u'(0) = u'(2) = 0
d = domain(0,2);
x = chebfun(@(x) x, d);
L1 = diff(d,2)+diag(x);
L1.bc = 'neumann';
% Create a chebop representing the same operator
N = chebop(d);
N.op = @(u) diff(u,2)+x.*u;
N.bc = 'neumann';
% Convert to linop
L2 = linop(N);
% Compare the results
pass = norm(L1(10)-L2(10)) == 0;