259
260
261
262
263
264
265
266 DOUBLE PRECISION V1( 2 ), V2( 2 )
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287 DOUBLE PRECISION ZERO
288 parameter( zero = 0.0d+0 )
289
290
291
292 IF( v1( 1 ).GE.v2( 1 ) ) THEN
293 IF( v1( 1 ).NE.zero )
294 $ v1( 2 ) = v1( 2 ) + ( v2( 1 ) / v1( 1 ) )**2 * v2( 2 )
295 ELSE
296 v1( 2 ) = v2( 2 ) + ( v1( 1 ) / v2( 1 ) )**2 * v1( 2 )
297 v1( 1 ) = v2( 1 )
298 END IF
299
300 RETURN
301
302
303