#pragma TextEncoding = "UTF-8"
#pragma rtGlobals=3		// Use modern global access method and strict wave access.
#pragma rtglobals=2

Function VelocityProfile()
	Wave ParameterWave
	Wave positionsy
	Wave positionsx
	Wave velocityP 

Make/O/N=200, positionsy; SetScale /P x 0, .0000001, positionsy; positionsy=x
Make/O/N=800, positionsx; SetScale /P x -.00002, .00000005, positionsx; positionsx=x
Make/O/N=1600, velocityP
 
velocityP = 0

	Variable Pressure = ParameterWave[0]
	Variable viscosity = ParameterWave[1]
	Variable H = ParameterWave[2]
	Variable L = ParameterWave[3]
	Variable w = ParameterWave[4]
	Variable n0 = ParameterWave[5]
	Variable n_max = ParameterWave[6]
	Variable step = ParameterWave[7]
	
variable n 


for (n=1; n<n_max; n+=2)
velocityP += (Pressure/(viscosity*L))*4*(H^2)/(pi)^3*(1/n^3)*(1-(cosh(n*pi*positionsx/H))/(cosh(n*pi*w/(2*H)))*sin(n*pi*(positionsy)/H))

endfor 


end