image thumbnail
from Adaptive residual subsampling for Radial Basis Functions by Alfa Heryudono
This simple codes solve initial-boundary value problems particularly Burgers' equation adaptively

Jburgers(t,u,epsilon,D1,D2)
function J = Jburgers(t,u,epsilon,D1,D2)
% Jacobian
J = epsilon*D2 - (diag(D1*u) + repmat(u,1,length(u)).*D1);

Contact us at files@mathworks.com