%DURCHFLUSS Q

load StrangK_fertig.txt
who
ergebnis = cumsum(StrangK_fertig)
Q = ergebnis(1:4)
save('durchfluss.txt', 'Q','-ascii')

%FLIESSGESCHWINDIGKEIT v

load StrangL_fertig.txt
durchmesser = StrangL_fertig(:,2)
A = pi/4 * (durchmesser.^2)
load durchfluss.txt
Q = durchfluss


v = (Q/1000) ./ (A/1000000)

save ('fliessgeschwindigkeit.txt', 'v', '-ascii')

%LAMBDA

sigma = 1.206 * 10^-6
lambda1 = 0
load fliessgeschwindigkeit.txt
load StrangL_fertig.txt
v = fliessgeschwindigkeit(1:4,1)
k = StrangL_fertig(:,4)
d = StrangL_fertig(:,2) / 1000
e = d.*v
f = k./d
j = f'
Re = e / sigma


    
   lam = - 2 * log10((2.51/(Re.*lam.^0.5 )) + f .* (1/3.71))

    
    

lam2 = 1./(lam.^0.5)
lambda = lam2'



%REIBUNGSGEFAELLE

g = 9.81
I = lambda .* 1./d .* ((v.^2)./(2*g))

%REIBUNGSVERLUSTHÖHE
laenge = StrangL_fertig(:,3)
hr = I .* laenge


%EXPORT
titel1 = 'N' 
titel2 = 'Re'
titel3 = 'Laenge'
titel4 = 'hr'
titel5 = 'v'
titel6 = 'Q'

N = StrangL_fertig(:,1)
Z = [N Re lambda hr v Q ]


dat = fopen('ergebnis.txt', 'w')
fprintf(dat, '%12s%12s\t%12s\t%12s\t%12s\t%12s\n', titel1, titel2, titel3, titel4, titel5, titel6)
fclose(dat);

dlmwrite ('ergebnis.txt', Z , 'delimiter','\t', 'precision','%12.4f', '-append')