Changeset 2878 for production/sydney_2006/find_max.py

Timestamp:

May 16, 2006, 2:06:46 PM (18 years ago)

Author:

sexton

Message:

intro updates

File:

• production/sydney_2006/find_max.py (modified) (3 diffs)

production/sydney_2006/find_max.py

@@ -6 +6 @@
 from math import exp, cosh
 import Numeric
+
 
 # Grilli and Watts 2005 example
@@ -33 +34 @@
 tol = 0.001
 c = 0
-direction = 'positive'
+direction = 'negative'
 while c < 100000 and deriv > 0:
-    deriv = dfuncdx(x,1.0)
+    deriv = dfuncdx(x,0.83)
     if deriv < 0: xstar = x
     if direction == 'positive': x += step
@@ -50 +51 @@
 test = 0.0
 for xi in x:
-    testi = func(xi,0.0,1.0)
+    testi = func(xi,0.0,0.83)
     if direction == 'positive':
         if testi > test: