306
307
308
309
310
311
312
313 REAL X, Y
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333 REAL ONE, ZERO
334 parameter( one = 1.0e+0, zero = 0.0e+0 )
335
336
337 REAL W, Z
338
339
341
342
343
346
347 IF( z.EQ.zero ) THEN
348 x = w
349 ELSE
350 x = w*sqrt( one+( z / w )**2 )
351 END IF
352
353 RETURN
354
355
356