Private Sub CommandButton1_Click()

CommandButton1.Caption = "Zn"

CommandButton1.BackColor = vbCyan

CommandButton1.Font = vbbold

Dim L As Double

Dim he As Double

Dim Zn1 As Double

Dim n As Double

Dim corr As Double

Dim Delta As Double

Dim Zn2 As Double

Dim cotZn1 As Double

Dim cotZn2 As Double

Dim error1 As Double

Dim error2 As Double

L = Sheet1.Range("B2")

he = Sheet1.Range("A2")

m = Val(Sheet1.Range("C2"))

Sheet1.Range("C4") = "Zn error"

Sheet1.Range("C4").Interior.ColorIndex = 8

For n = 0 To m - 1

Sheets("sheet1").Cells(n + 5, 1) = "Z" + Str(n + 1)

Sheet1.Cells(n + 5, 1).Interior.ColorIndex = 8

Zn1 = n * 3.1415926

Sheets("sheet1").Cells(n + 5, 2) = 0.0001

For corr = 1 To 1000

Delta = 0.001 * 3.1415926

Zn1 = Zn1 + Delta

Zn2 = Sheets("sheet1").Cells(n + 5, 2)
cotZn1 = 1 / Tan(Zn1)

cotZn2 = 1 / Tan(Zn2)

error1 = (cotZn1 - (Zn1 / (he * L))) ^ 2

error2 = (cotZn2 - (Zn2 / (he * L))) ^ 2

If error1 < error2 Then

Sheets("sheet1").Cells(n + 5, 2) = Val(Zn1)

Sheets("sheet1").Cells(n + 5, 3) = Val(error1)

End If

Next corr

Next n

End Sub

Private Sub CommandButton2_Click()

CommandButton2.Caption = "find Alpha"

CommandButton2.BackColor = vbGreen

CommandButton2.Font = vbbold

Dim lastrow, Alpha, i, t, L, x, Ua, Uo, m, j, Zn, part1, part2, part3, sum, u, Error, sumerror As Double

lastrow = Sheets("sheet1").Cells(Rows.Count, 9).End(xlUp).Row

Sheets("sheet1").Cells(lastrow + 1, 17) = 10000

For Alpha = 0.00000001 To 0.01 Step 0.00000001

If Alpha > 0.000001 Then GoTo limit

sumerror = 0

For i = 2 To lastrow

t = Sheets("sheet1").Cells(i, 10)

L = Sheet1.Range("B2")

x = Sheet1.Range("F4")

Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

m = Val(Sheet1.Range("C2"))
sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i, 16) = u

Erro = ((Sheets("sheet1").Cells(i, 16) - Sheets("sheet1").Cells(i, 9)) ^ (2))

Sheets("sheet1").Cells(i, 17) = Erro

sumerror = sumerror + Erro

Sheet1.Range("S2") = sumerror

Next i

L = Sheets("sheet1").Cells(lastrow + 1, 17)

If Sheet1.Range("S2") < L Then

Sheets("sheet1").Cells(lastrow + 1, 17) = sumerror

Sheet1.Range("N2") = Alpha

End If

Next Alpha

limit:

For Alpha = 0.000001 To 0.01 Step 0.000001

If Alpha > 0.0001 Then GoTo limit1

sumerror = 0

For i = 2 To lastrow

t = Sheets("sheet1").Cells(i, 10)

L = Sheet1.Range("B2")

x = Sheet1.Range("F4")
Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i, 16) = u

Erro = ((Sheets("sheet1").Cells(i, 16) - Sheets("sheet1").Cells(i, 9)) ^ (2))

Sheets("sheet1").Cells(i, 17) = Erro

sumerror = sumerror + Erro

Sheet1.Range("S2") = sumerror

Next i

L = Sheets("sheet1").Cells(lastrow + 1, 17)

If Sheet1.Range("S2") < L Then

Sheets("sheet1").Cells(lastrow + 1, 17) = sumerror

Sheet1.Range("N2") = Alpha

End If

Next Alpha

limit1:

For Alpha = 0.0001 To 0.01 Step 0.0001

sumerror = 0

For i = 2 To lastrow

t = Sheets("sheet1").Cells(i, 10)
L = Sheet1.Range("B2")

x = Sheet1.Range("F4")

Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i, 16) = u

Erro = ((Sheets("sheet1").Cells(i, 16) - Sheets("sheet1").Cells(i, 9)) ^ (2))

Sheets("sheet1").Cells(i, 17) = Erro

sumerror = sumerror + Erro

Sheet1.Range("S2") = sumerror

Next i

L = Sheets("sheet1").Cells(lastrow + 1, 17)

If Sheet1.Range("S2") < L Then

Sheets("sheet1").Cells(lastrow + 1, 17) = sumerror

Sheet1.Range("N2") = Alpha

End If

Next Alpha

End Sub

Private Sub CommandButton3_Click()

CommandButton3.BackColor = vbGreen

CommandButton3.Font = vbbold

Dim ac, lower, upper, lastrow, Alpha, i, t, L, x, Ua, Uo, m, j, Zn, part1, part2, part3, sum, u, Erro, sumerror
As Double

lastrow = Sheets("sheet1").Cells(Rows.Count, 9).End(xlUp).Row

Sheets("sheet1").Cells(lastrow + 1, 17) = 10000

ac = Sheet1.Range("N2")

lower = ac / 100

upper = ac * 3

For Alpha = lower To upper Step lower

sumerror = 0

For i = 2 To lastrow

t = Sheets("sheet1").Cells(i, 10)

L = Sheet1.Range("B2")

x = Sheet1.Range("F4")

Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i, 16) = u
Erro = ((Sheets("sheet1").Cells(i, 16) - Sheets("sheet1").Cells(i, 9)) ^ (2))

Sheets("sheet1").Cells(i, 17) = Erro

sumerror = sumerror + Erro

Sheet1.Range("S2") = sumerror

Next i

L = Sheets("sheet1").Cells(lastrow + 1, 17)

If Sheet1.Range("S2") < L Then

Sheets("sheet1").Cells(lastrow + 1, 17) = sumerror

Sheet1.Range("N2") = Alpha

lastrow = Sheets("sheet1").Cells(Rows.Count, 9).End(xlUp).Row

sumerror = 0

For i = 2 To lastrow

t = Sheets("sheet1").Cells(i, 10)

L = Sheet1.Range("B2")

x = Sheet1.Range("F4")

Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i, 11) = u

Erro = ((Sheets("sheet1").Cells(i, 11) - Sheets("sheet1").Cells(i, 9)) ^ (2))

sumerror = sumerror + Erro

Sheets("sheet1").Cells(lastrow + 1, 12) = sumerror

Next i

End If

Next Alpha

Sheets("sheet1").Cells(lastrow + 1, 17) = Sheet1.Range("S2")

End Sub

Private Sub CommandButton4_Click()

CommandButton4.Caption = "Draw Suction vs. Time"

CommandButton4.BackColor = vbYellow

CommandButton4.Font = vbbold

Dim Alpha, i, t, L, x, Ua, Uo, m, j, Zn, part1, part2, part3, sum, u As Double

Alpha = Sheet1.Range("N2")

L = Sheet1.Range("B2")

x = Sheet1.Range("F4")

Ua = Sheet1.Range("F2")

Uo = Sheet1.Range("G2")

t=0

For i = 1 To 10

t = t + 100

Sheets("sheet1").Cells(i + 25, 6) = t

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i + 25, 7) = u

Next i

t = 1000

For i = 1 To 9

t = t + 1000

Sheets("sheet1").Cells(i + 35, 6) = t

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m

Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i + 35, 7) = u

Next i

t = 10000

For i = 1 To 9

t = t + 10000

Sheets("sheet1").Cells(i + 44, 6) = t

m = Val(Sheet1.Range("C2"))

sum = 0

For j = 1 To m
Zn = Sheets("sheet1").Cells(j + 4, 2)

part1 = (2 * (Uo - Ua) * Sin(Zn)) / (Zn + Sin(Zn) * Cos(Zn))

part2 = Exp(-1 * ((Zn) ^ (2)) * Alpha * t / ((L) ^ (2)))

part3 = Cos(Zn * x / L)

sum = sum + part1 * part2 * part3

Next j

u = Ua + sum

Sheets("sheet1").Cells(i + 44, 7) = u

Next i

End Sub