% m28chebfun.m - Chebfun solution of the same equation L = chebop(0,3); L.op = @(t,u) diff(u) - u*tanh(exp(t)*sin(5*t)); L.lbc = 1; u = L\0 plot(u) err = u(3) - 1.403737560405910; title(['err = ' num2str(err)],FS,16) % For further explorations try length(u) and plot(u,'.'). % Then repeat with L.domain = [0 6]. Such use of a global % polynomial of very high degree is pretty extreme. % For a more local representation, rerun after "splitting on".