%Turbular viscosity
clear
clc
clf
[dist1,visc1]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',8,'headerlines',2);
[dist2,visc2]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',8,'headerlines',11);
[dist3,visc3]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',7,'headerlines',20);
[dist4,visc4]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',17,'headerlines',28);
[dist5,visc5]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',15,'headerlines',46);
[dist6,visc6]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%f%f',16,'headerlines',62);
hold on
plot(dist1,visc1,'.r')
plot(dist2,visc2,'.g')
plot(dist3,visc3,'.b')
plot(dist4,visc4,'.m')
plot(dist5,visc5,'.k')
plot(dist6,visc6,'.c')

xy=[dist1,visc1;dist2,visc2;dist3,visc3;dist4,visc4;dist5,visc5;dist6,visc6];

[x, indeks] = sort(xy(:, 1));
xy1 = xy(indeks, :);

X=xy1(:,1);
Y=xy1(:,2);
P=polyfit(X,Y,5);
Y1=polyval(P,X);

plot(X,Y1);

hold off
xlabel('distance')
ylabel('velocity')
[lan11,lan12,lan13]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',1);
[lan21,lan22,lan23]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',10);
[lan31,lan32,lan33]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',19);
[lan41,lan42,lan43]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',27);
[lan51,lan52,lan53]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',45);
[lan61,lan62,lan63]=textread('H:\chwila\TS\MatLab\turb_visc.dat','%s%s%s',1,'headerlines',61);
lan1=strcat(lan11,lan12,lan13);
lan2=strcat(lan21,lan22,lan23);
lan3=strcat(lan31,lan32,lan33);
lan4=strcat(lan41,lan42,lan43);
lan5=strcat(lan51,lan52,lan53);
lan6=strcat(lan61,lan62,lan63);
legend([lan1,lan2,lan3,lan4,lan5,lan6],'location','SouthEast')
