3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %%% Dispersion bei erweiterter Regularisierung der BL %%%
5 %%% u_t + f(u)_{x} = epsilon (D(u)u_{x})_{x} + epsilon^{2}*tau*u_{xxt} %%%
6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7 % eigentlich muesste es
8 % u_t + f(u)_{x} = epsilon [D(u)(u_{x} + epsilon*tau*u_{xt})]_{x}
10 % WIRD SO IN D(U)allg VERWENDET
19 D(k) = (u(k)^2*(1-u(k))^2)/(u(k)^2+M*(1-u(k))^2); % M=2 bis jetzt fest, entsprechend aendern!!!