diff options
author | Santosh G. Vattam | 2009-10-08 22:48:59 +0530 |
---|---|---|
committer | Santosh G. Vattam | 2009-10-08 22:48:59 +0530 |
commit | 01824f31456fb6e6068ca72163421fca19d6f9ce (patch) | |
tree | 44f6ad4ad4a23246ff176d73ecbf19ca1aa02ab9 /day1/exercise/roots.py | |
parent | 1c50afc1e22012a9bdb3d2b67f8d369ce2fea6f4 (diff) | |
parent | 9936d8647f15e656351fc2551b08bfb91b38590f (diff) | |
download | workshops-more-scipy-01824f31456fb6e6068ca72163421fca19d6f9ce.tar.gz workshops-more-scipy-01824f31456fb6e6068ca72163421fca19d6f9ce.tar.bz2 workshops-more-scipy-01824f31456fb6e6068ca72163421fca19d6f9ce.zip |
Branches merged.
Diffstat (limited to 'day1/exercise/roots.py')
-rw-r--r-- | day1/exercise/roots.py | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/day1/exercise/roots.py b/day1/exercise/roots.py new file mode 100644 index 0000000..c6bf43e --- /dev/null +++ b/day1/exercise/roots.py @@ -0,0 +1,28 @@ +import math + +def linspace(a, b, N): + lns = [] + step = (float(b) - float(a)) / float(N - 1) + print step + for i in range(N): + lns.append(a + i*step) + + return lns + +def sinsin_func(): + x = linspace(0, 5, 11) + sin_list = [] + for i in x: + sin_list.append(math.sin(i) + math.sin(10*i)) + + return sin_list + +def find_root_range(): + sin_list = sinsin_func() + for i, sins in enumerate(sin_list): + if (sin_list[i] > 0 and sin_list[i+1] < 0) or (sin_list[i] > 0 and sin_list[i+1] < 0): + print "Roots lie between: %f and %f" % (sin_list[i], sin_list[i+1]) + if sin_list[i] == 0: + print "%f is the root" % sin_list[i] + +find_root_range() |