program ONEBOX
        real m,nubar,wvl
        integer r
        r0=0.14
        write(6,100)r0
100	format(1H1,1X,'PARTICLE IN A BOX MODEL',//,1X ,'RO= ',f6.3,'NM')
	print*,'GIVE ME THE NUMBER OF DOUBLE BONDS IN THE MOLECULE'
1	read(5,*)r
	if(r.eq.0) stop
	write(6,1001)r
1001	format(///,1X,'THIS CONJUGATED SYSTEM HAS ',I2,' DOUBLE BONDS',/)
	rr=r0*1.0E-9
	h=6.67E-34
	m=9.107E-31
	c=2.997E+8
	xl=(2*r+1)*rr
	nubar=(2*r+1)*h/(8*xl*xl*c*m)
	nubar=nubar*1.0e-2
	wvl=1.0e+7/nubar
	write(6,1002)nubar,wvl
1002	format(1X,'FOR THE LONGEST WAVELENGTH TRANSITION',/,1X,'NUBAR=',
     1	f8.0,' 1/CM',/,1X,'WAVELENGTH =',f6.1,' NM',////)
	go to 1
	end