не могу понять на что ругается

cls:rem Lame Variant 1 Student MPO-18-01 Abdrakhmanov Rinat
defdbl a,b,e,n,r,s:open "result.txt" for output as#1 
asm=5:bsm=15:k=1.05:E=2.1e11:nu=0.3:sigt=0.25e8:J=1 
Input "Vvedite J" J
a=asm/100:b=bsm/100:sigdop=sigt/k:pmax=sigdop*(b^2-a^2)/2/b^2 
R$="## ####.# ###.# ####.# ####.# ###.### -#.##^^^^ -#.##^^^^ -#.##^^^^ #.##^^^^" 
? " r sigr sigfi sigz sigMises u epsr epsfi epsz epsMises " 
? " sm MPa MPa MPa MPa mm" 
? #1, " r sigr sigfi sigz sigMises u epsr epsfi epsz epsMises " 
? #1, " sm MPa MPa MPa MPa mm" 
shag=(b-a)/10:for r=a to b+shag/2 step shag:gosub NDS 
? using R$;r*1e2,sigr/1e6,sigfi/1e6,sigz/1e6,sigmis/1e6,u*1e3,epsr,epsfi,epsz,epsmis 
? #1, using R$;r*1e2,sigr/1e6,sigfi/1e6,sigz/1e6,sigmis/1e6,u*1e3,epsr,epsfi,epsz,epsmis 
next:gosub Vyvod:close 1:END 
NDS: 
pa=pmax:pb=0:AA=(pa*a^2-pb*b^2)/(b^2-a^2):BB=(pa-pb)*a^2*b^2/(b^2-a^2) 
sigr=AA-BB/r^2:sigfi=AA+BB/r^2: 
if J=1 then sigz=0 
if J=2 then sigz=2*nu*AA 
if J=3 then sigz=(pa*a^2-pb*b^2)/(b^2-a^2) 
u=1/E*(AA*(1-nu)*r+BB*(1+nu)/r-nu*r*sigz) 
epsr=1/E*(AA*(1-nu)-BB*(1+nu)/r^2-nu*sigz) 
epsr1=1/E*(sigr-nu*(sigfi+sigz)) 
epsfi=1/E*(AA*(1-nu)+BB*(1+nu)/r^2-nu*sigz): 
epsfi1=1/E*(sigfi-nu*(sigr+sigz)) 
epsz=1/E*(sigz-2*nu*AA):epsz1=1/E*(sigz-nu*(sigr+sigfi)) 
epsmis=1/sqr(2)*sqr((epsr-epsfi)^2+(epsfi-epsz)^2+(epsz-epsr)^2) 
er=E*(epsr+epsfi+epsz)-(1-2*nu)*(sigr+sigfi+sigz) 
sigm=(sigr+sigfi+sigz)/3 
sigmis=sqr(1.5)*sqr((sigr-sigm)^2+(sigfi-sigm)^2+(sigz-sigm)^2) 
sigmis1=1/sqr(2)*sqr((sigr-sigfi)^2+(sigfi-sigz)^2+(sigz-sigr)^2) 
RETURN 
Vyvod: 
? using " sigdop = ####";sigdop/1e6;:?" MPa" 
? using " pmax = ####.#";pmax/1e5;:? "atm" 
SELECT CASE J 
CASE 1 
? "J=1 Plane stress: sigz=0":? #1, "J=1 Plane stress: sigz=0" 
CASE 2 
? "J=2 Sosud vuskogo davleniya":? #1, "J=2 Plane strain: epsz=0" 
CASE 3 
? "J=3 Plane strain: epsz=0":? #1, "J=3 Sosud vuskogo davleniya" 
END SELECT 
RETURN

rem Lame Variant 1 Student MPO-18-01 Abdrakhmanov Rinat
cls
defdbl a,b,e,n,r,s:open "result.txt" for output as#1