function [x,res,t]=Burger_decentre4pt(model)


dx=2*model.xmax/(model.nPt-1);
dt=model.courant*dx;
s=model.viscosite*dt/(dx*dx);
t=(1:model.ntemps)*dt;

i=2:(model.nPt-1);

[x,u]=conditions_initialles(model.xmax, model.nPt);
ud=0*u;
res=zeros(model.ntemps, model.nPt);

for j=1:model.ntemps
	ud(1)=u(1);
	ud(end)=u(end);
    
	ud(i)=u(i)+s*(u(i-1) - 2*u(i) + u(i+1))...
        -model.courant*(0.5*(F(u,i+1)-F(u,i-1))+model.q/3*(-3*F(u,i-1)+3*F(u,i)-F(u,i+1)+F(u,i-2)));
    
    
	res(j,:)=ud;
	u=ud;
end
