280 do 100 mrs = 1, nrs 281 C 282 do 10 mq = 1, nq 283 do 10 mi = 1, morb 284 10 xrsiq(mi, mq) = zero 285 do 40 mp = 1, np 286 do 30 mq = 1, mp 287 mrspq = mrspq+1 288 val = xrspq(mrspq) 289 if (val .eq. zero) 289 goto 30 290 do 20 mi = 1, morb 291 xrsiq(mi, mq) = xrsiq(mi, mq)+val*v(mp, mi) 292 xrsiq(mi, mp) = xrsiq(mi, mp)+val*v(mq, mi) 293 20 continue 294 30 continue 295 40 continue 296 C 297 mrsij = mrsij0298 do 90 mi = 1, morb 299 do 50 mj = 1, mi 300 50 xij(mj) = zero 301 do 70 mq = 1, nq 302 val = xrsiq(mi, mq) 303 if (val .eq. zero) 303 goto 70 304 do 60 mj = 1, mi 305 60 xij(mj) = xij(mj)+val*v(mq, mj) 306 70 continue 307 do 80 mj = 1, mi 308 mrsij = mrsij+1 309 80 xrsij(mrsij) = xij(mj) 310 90 continue 311 C 312 mrsij0 = mrsij0+nrs 313 100 continue