plotting non-functions

plot_nonfun.py

# something that doesn't pass the vertical line test, comes to us unordered.... we want to piecewise-interpolate.

import numpy as np
# partition our data so that each segment DOES pass the vertical line test.
c = np.where(a[:,1] >= 0.5)[0]
d = np.where(a[:,1] < 0.5)[0]

above = a[c,:] # copies of a
below = a[d,:]

print(above)
print(below)

# sort them
ab = np.argsort(above[:,0])
be = np.argsort(below[:,0])

# notice how the sorting works.
print(above[ab])
print(below[be])

plt.scatter(a[:,0], a[:,1])
plt.plot(above[ab,0], above[ab,1])
plt.plot(below[be,0], below[be,1])
plt.show()