На 40 строке выдает ошибку "Неверный вызов функции"

CLS
PRINT " Ras4et xolodnou prokatku na stane 2030 no MeToduke KopoLeBa"
FOR i = 1 TO 80
PRINT "=";
NEXT i
 
PRINT " napametpbl staHa:"
DIM n AS INTEGER
INPUT "N=", n
DIM r, sigma0, sigma1, sigma02, a1, n1, b
INPUT "R=", r
INPUT "sigma0(nepednee Hat9}|{EHue)=", sigma0
INPUT "sigma1(3adHee Hat9}|{EHue)=", sigma1
 
PRINT ""
PRINT " napametpbl ctaJIu:"
INPUT "sigma02(ucxodHoe)=", sigma02
INPUT "a=", a1
INPUT "n1=", n1
INPUT "b=", b
 
DIM h(n), dh(n), hcp(n), esum(n), sigmat(n), Tcp(n), Dz(n), l(n), mvh(n), mvuh(n), Ksi(n), Pcp(n), P(n)
 
INPUT "h0(Ha4aJIbHa9 ToJIshina)=", h(0)
PRINT "O6}|{aTu9 no KJIeT9M:"
FOR i = 1 TO n
INPUT "E=", E(i)
NEXT i
 
PRINT "koef treHu9 no kJIeT9M"
FOR i = 1 TO n
INPUT "mu=", mu(i)
NEXT i
 
FOR i = 1 TO n
h(i) = h(i - 1) * E(i)
dh(i) = h(i - 1) - h(i)
hcp(i) = (h(i) + h(i - 1)) / 2
esum(i) = (h(0) - h(i)) * 100 / h(0)
sigmat(i) = sigma02 + (a1 * esum(i) ^ n1)
NEXT i
 
FOR i = 1 TO n
sigmat(0) = sigma02
Tcp(i) = 1.15 * (sigmat(i) + sigmat(i - 1)) / 2
sigmat(0) = sigma02
Dz(i) = sigmat(i - 1) / sigmat(i)
 
p0 = 0
x2 = 0
10 : l(i) = SQR(r * dh(i) + x2(i) ^ 2) + x2(i)
mvh(i - 1) = 2 * mu(i) * l(i) / (h(i - 1))
mvuh(i) = 2 * mu(i) * l(i) / (h(i))
Ksi(i) = .5 * (1 - E(i)) / (1 - .5 * E(i)) * (1 - (1 / mvh(i - 1)) * ln(Dz(i) - ((sigma1) / (Tcp(i))) / (1 - (sigma0) / (Tcp(i)))))
Pcp(i) = Tcp(i) * ((Dz(i) - ((sigma1) / (Tcp(i))) * (1 / mvuh(i)) * (EXP(mvuh(i) * (Ksi(i))) - 1) + (1 - (sigma0) / (Tcp(i))) * (EXP(mvh(i - 1) * (1 - Ksi(i)) - 1))))
 
IF ABS(Pcp(i) - p0) / Pcp(i) > .01 THEN
x2 = Pcp(i) * r / 95000
p0 = Pcp(i): GOTO 10
END IF
P(i) = Pcp(i) * b * l(i) / (10 ^ 6)
NEXT i
 
FOR i = 1 TO 80
PRINT "=";
NEXT i
PRINT ""
 
PRINT " h(mm) dh(mm)"
FOR i = 1 TO n
PRINT ; h(i); dh(i)
NEXT i
PRINT ""
PRINT "      mvh     mvuh      l(mm)     sigmat(H/mm2)    Ksi      Pcp(H/mm2)      P(MH)"
FOR i = 1 TO n
PRINT ; mvh(i - 1); mvuh(i); "     "; l(i); sigmat(i); "      "; Ksi; Pcp(i); P(i)
NEXT i
END