clear all
disp('Metodo FB per l''equazione del trasporto con funzione iniziale g(x) in [0,1]  ')
disp('e valore 0 sull''intervallo [1,4]')
k=input('k=');   %passo temporale
N=input('N=');  %numero di passi temporali
J=input('J=');  %numero di passi spaziali in [0,1]
h=1/J;      %passo spaziale
t=0:k:1;
x=0:h:1;
g=(4*x.*(1-x)).^2;
g=[g,zeros(1,3*J)];
%b=(t.*(1-t)).^2;
C=k/h
s=length(g);
A=eye(s)-C*eye(s)+C*diag(ones(s-1,1),-1);
u(:,1)=g';
for n=2:N
    u(:,n)=A*u(:,n-1);
    %u(1,n)=u(1,n)+C*b(n-1);
    n=n+1;
end
mesh(u)