2013-12-09 3 views
8

Я пытаюсь реализовать простой метод итерационного метода Ньютона с использованием Excel VB (я никогда не использовал VB.)Loop без Do ошибки

я получаю ошибку loop without a Do, и я не могу понять, что я я делаю неправильно здесь.

Я пытаюсь найти корни функции z^3 - z^2 - (B^2 + B - A)z - A*B, называемой коэффициентом сжимаемости.

Мой источник MSN

Function zCalculation(ByVal temp As Double, ByVal press As Double) As Double 

Dim tempCr As Double 
Dim pressCr As Double 
Dim A As Double 
Dim B As Double 

tempCr = temp/238.5 

pressCr = press/547.424092 

A = pressCr/tempCr 
A = A/(9 * (2^(1/3) - 1)) 
B = pressCr/tempCr 
B = B * (2^(1/3) - 1)/3 



Dim zNot As Double 
Dim counter As Integer 
counter = 0 
zNot = 1# 

Do 
    counter = counter + 1 

    zNot = zNot - (zNot^3 - zNot^2 - (B^2 + B - A) * zNot - A * B)/(3 * zNot^2 - 2 * zNot - (B^2 + B - A)) 
    If counter > 1000 Then 
     Exit Do 

Loop Until eval(zNot, A, B) < 0.000001 


zCalculation = zNot 


End Function 

перерыв

Function eval(ByVal z As Double, ByVal A As Double, ByVal B As Double) As Double 

    eval = z^3 - z^2 - (B^2 + B - A) * z - A * B 

End Function 

ответ

12

Вам нужно:

End If 

в вашем коде.

+0

Dang ..... это все. Благодаря! Не могу поверить, что я пропустил это. –

+3

Не ваша ошибка ............ Сообщение об ошибке Excel в этом случае неверно. –

+1

не могу поверить, что я действительно искал в stackoverflow перед проверкой этого – Tascalator

1

Вы можете попробовать:

Function zCalculation(ByVal temp As Double, ByVal press As Double) As Double 

    Dim tempCr As Double 
    Dim pressCr As Double 
    Dim A As Double 
    Dim B As Double 

    tempCr = temp/238.5 

    pressCr = press/0.546789 

    A = pressCr/tempCr 
    A = A/(9 * (2^(1/3) - 1)) 
    B = pressCr/tempCr 
    B = B * (2^(1/3) - 1)/3 



    Dim zNot As Double 
    Dim counter As Integer 
    counter = 0 
    zNot = 1# 

    Do 
     counter = counter + 1 

     zNot = zNot - (zNot^3 + zNot^2 - (B^2 + B - A) * zNot - A * B)/(3 * zNot^2 + 2 * zNot - (B^2 + B - A)) 
     If counter > 1000 Then 
     Exit Do 
     End if ' <--- Here 

    Loop Until eval(zNot, A, B) < 0.000001 

    zCalculation = zNot 
End Function 
+0

, хотя у вас было это правильно, а Гэри выложил первый. Благодаря! –

+0

Без проблем человек. Я рад помочь вам. – Makah

-2
Sub datacalculationsandformat() 
Dim row As Integer 
row = 2 
Do While Cells(row, 2) <> "" 
Cells(row, 3).Value = Cells(row, 2).Value * 0.3 
Cells(row, 4).Value = Cells(row, 2) * 0.1 
Cells(row, 5).Value = Cells(row, 2).Value + Cells(row, 3).Value + Cells(row, 4).Value 
If Cells(row, 5).Value >= 8000 Then 
Worksheets("Sheet1").Cells(row, 5).Font.Bold = True 
row = row + 1 
Loop