function Q=f0(t); %forzante applicata nel nodo 0

%Q=10; %forzante costante
%return

load forzante;      %0=forzante costante
                    %1=flash
                    %2=campione unitario
                    %3=per condizioni di dirichlet
switch forzante
    case 0,
       Q=40*ones(length(t),1); 
    case 1,
        t=t.*5e-04;
        W=6.6296e+006;
        sig= 0.0106;
        Q=(W*t).*exp(-sqrt(t)/sig);
    case 2,
        load dati;
        load eulero;
        heat_val = 1e3;
        if eulero==0 %esplicito
            for j=1:length(t)
                if abs(t(j))<eps
                    Q(j)=1*heat_val;
                else
                    Q(j)=0;
                end
            end
        else %implicito
            for j=1:length(t)
                if abs(t(j))<eps
                    Q(j)=1*heat_val;
                else
                    Q(j)=0;
                end
            end
        end
    case 3,
       Q=zeros(length(t),1);
       Q(1)=40;
end
