#!/usr/bin/env python

from math import *
M_PI = pi
true=1
false=0

from epix import *

####################################################################
# Converted from stereo-left.xp
####################################################################

def u1(t):
  return t

def u2(t):
  return t*t

def u3(t):
  return t*t*t

def zero(t):
  return 0

N=60

unitlength("1pt")
bounding_box(P(-1,-1), P(1,1))
picture(P(150,150))

begin()

viewpoint(1, 1.75, 0.5)
camera.range(20)

#  coordinate axes
arrow(P(-1,0,0), P(1,0,0))
arrow(P(0,0,-1), P(0,0,1))

bold()
blue()
plot(u1, zero, u3, -1, 1, N)

t=-1
for i in range(0,21):
  line(P(u1(t), 0, u3(t)), P(u1(t), u2(t), u3(t)))
  t += 0.1

label(P(1,0,1), P(2,2), "$x$-$z$ plane")

red()
plot(zero, u2, u3, -1, 1, N)
t = -1
for i in range(0,21):
  line(P(0, u2(t), u3(t)), P(u1(t), u2(t), u3(t)))
  t += 0.1

label(P(0,1,1), P(2,2), "$y$-$z$ plane")

magenta()
plot(u1, u2, zero, -1, 1, N)

t = -1
for i in range(0,21):
  line(P(u1(t), u2(t), 0), P(u1(t), u2(t), u3(t)))
  t += 0.1

label(P(1,1,0), P(-45,-12), "$x$-$y$ plane")
black()

# The cubic
plot(u1, u2, u3, -1, 1, N)

plain()
arrow(P(0,0,0), P(0,1,0))

end()