import plotly.offline as py
from plotly.graph_objs import *
import numpy as np
from numpy import pi, cos, sin, exp, log, sqrt
py.init_notebook_mode()

def curve(rfun,tmin=-2,tmax=2,tpts=200,color='black'):
    domain = np.linspace(tmin,tmax,tpts)
    r = [[rfun(t)[i] for t in domain] for i in range(3)]
    trace = Scatter3d(x=r[0],y=r[1],z=r[2],mode='lines',
                      line=Line(color=color,width=3))
    return(trace)

def snake(u,v):
    lu, lv = len(u), len(v)
    path = []

    i, j = 0, 0
    istep, jstep = 1, 1
    while(i < lu):
        while(0 <= j < lv):
            path.append((u[i],v[j]))
            j += jstep
        j -= jstep
        i += istep
        jstep *= -1

    i -= istep
    i -= istep
    istep *= -1
    while(0 <= j < lv):
        while(0 <= i < lu):
            path.append((u[i],v[j]))
            i += istep
        i -= istep
        j += jstep
        istep *= -1
    return path

def surface(rfun,tmin=-2,tmax=2,tpts=20,umin=-2,umax=2,upts=20,color='black'):
    t, u = np.linspace(tmin,tmax,tpts), np.linspace(umin,umax,upts)
    path = snake(t,u)
    r = [[rfun(t,u)[i] for (t,u) in path] for i in range(3)]
    trace = Scatter3d(x=r[0],y=r[1],z=r[2],mode='lines',
                      line=Line(color=color,width=3))
    return(trace)

f = lambda x,y: x*x*x/(x*x+y*y) if x!=0 or y!=0 else 0
x0=0
y0=0
xmin=-1
xmax=1
ymin=-1
ymax=1
s = surface(lambda t,u: (t,u,f(t,u)),color='green',tmin=xmin,tmax=xmax,umin=ymin,umax=ymax) #surface z=f(x,y)
sy = curve(lambda t: (t,y0,f(t,y0)),color='red',tmin=xmin,tmax=xmax) #slice at y=y0
sx = curve(lambda t: (x0,t,f(x0,t)),color='blue',tmin=ymin,tmax=ymax) #slice at x=x0
sd = curve(lambda t: (t,t,f(t,t)),color='black',tmin=ymin,tmax=ymax) #diagonal slice
ss = curve(lambda t: (t,2*t,f(t,2*t)),color='brown',tmin=ymin/2,tmax=ymax/2) #slice at y/x slope 2
sn = curve(lambda t: (t,-t,f(t,-t)),color='purple',tmin=ymin,tmax=ymax) #slice at y/x slope -1
py.iplot(Figure(data=Data([s,sx,sy,sd,ss,sn])))
py.iplot(Figure(data=Data([sx,sy,sd,ss,sn])))

x0=0
y0=0
xmin=-.01
xmax=.01
ymin=-.01
ymax=.01
s = surface(lambda t,u: (t,u,f(t,u)),color='green',tmin=xmin,tmax=xmax,umin=ymin,umax=ymax) #surface z=f(x,y)
sy = curve(lambda t: (t,y0,f(t,y0)),color='red',tmin=xmin,tmax=xmax) #slice at y=y0
sx = curve(lambda t: (x0,t,f(x0,t)),color='blue',tmin=ymin,tmax=ymax) #slice at x=x0
sd = curve(lambda t: (t,t,f(t,t)),color='black',tmin=ymin,tmax=ymax) #diagonal slice
ss = curve(lambda t: (t,2*t,f(t,2*t)),color='brown',tmin=ymin/2,tmax=ymax/2) #slice at y/x slope 2
sn = curve(lambda t: (t,-t,f(t,-t)),color='purple',tmin=ymin,tmax=ymax) #slice at y/x slope -1
py.iplot(Figure(data=Data([s,sx,sy,sd,ss,sn])))
py.iplot(Figure(data=Data([sx,sy,sd,ss,sn])))