close all
clear all
% ***** read u
load u_v_time_4nodes_hybrid.dat
u1=u_v_time_4nodes_hybrid(:,1);
v1=u_v_time_4nodes_hybrid(:,2);
u2=u_v_time_4nodes_hybrid(:,3);
v2=u_v_time_4nodes_hybrid(:,4);
u3=u_v_time_4nodes_hybrid(:,5);
v3=u_v_time_4nodes_hybrid(:,6);
u4=u_v_time_4nodes_hybrid(:,7);
v4=u_v_time_4nodes_hybrid(:,8);

% time step=0.0049
dt=0.0049;
n=length(u1);
% compute time array
t=dt:dt:n*dt;

%%%%%%%%%%%%%%%% plotting section %%%%%%%%%%%%%%%%%%%%%%%%%%
% plot u
plot(t,u2)
hold
plot(t,u4,'r--')
xlabel('t')
ylabel('u')
handle=gca;
set(handle,'fontsize',[20])
%print u_time_hybrid.ps -deps
%
