# Standard 3D unit vectors
ivec = vector([1,0,0])
jvec = vector([0,1,0])
kvec = vector([0,0,1])

t=2
p = plot(ivec, thickness=t, color='blue')
p += plot(jvec, thickness=t, color='red')
p += plot(kvec, thickness=t, color='green')

# Vector
v1 = 1
v2 = 0.75
v3 = 0.5
v = vector([v1,v2,v3])
p += plot(v, thickness=t, color='orange')

# Vector components
p += point3d((v1,v2,v3), size=10, color='orange')
p += point3d((v1,0,0), size=10, color='orange')
p += point3d((0,v2,0), size=10, color='orange')
p += point3d((0,0,v3), size=10, color='orange')


# Right triangles
s = 0.2
p += polygon3d([[0,0,0],[v1,v2,v3],[v1,0,0]], opacity=0.1, color='blue')
p += polygon3d([[0,0,0],[s*v1,s*v2,s*v3],[s*v1,0,0]], opacity=0.5, color='blue')
p += polygon3d([[0,0,0],[v1,v2,v3],[0,v2,0]], opacity=0.1, color='red')
p += polygon3d([[0,0,0],[s*v1,s*v2,s*v3],[0,s*v2,0]], opacity=0.5, color='red')
p += polygon3d([[0,0,0],[v1,v2,v3],[0,0,v3]], opacity=0.1, color='green')
p += polygon3d([[0,0,0],[s*v1,s*v2,s*v3],[0,0,s*v3]], opacity=0.5, color='green')


p.show(frame=False)