function Vel(r)

global F H_bilr ID IEN Node r Nml_Vct

Velocity=H_bilr\F;

dt=10^(-5);
d_r=Velocity*dt;

Sz1=size(ID);
IDR=zeros(Sz1(1),Sz1(2));
for i=1:Sz1(2)
    for j=1:Sz1(1)
        P=ID(j,i);
        if P>0
            IDR(j,i)=d_r(P);
        end
    end
end

NND=Node(:,2:4)+IDR';
Node=[Node(:,1) NND];


scrsz = get(0,'ScreenSize');
figure('Position',[1 scrsz(4)/3 scrsz(3)/3 scrsz(4)/3]);
PIC=trisurf(IEN(:,2:4),Node(:,2),Node(:,3),Node(:,4));
colormap(gray)
axis equal
axis([-60*10^(-4) 60*10^(-4) -60*10^(-4) 60*10^(-4) -20*10^(-4) 20*10^(-4)])

ext=num2str(r);
Name=['T' ext];
saveas(PIC,Name,'jpg')

disp('Vel OK!')

