#!/usr/bin/python
from numpy import *
import pylab as p
#import matplotlib.axes3d as p3
import mpl_toolkits.mplot3d.axes3d as p3

# Здесь нужно как-то считать данные 

fig=p.figure()
ax = p3.Axes3D(fig)
# x, y, and z are 100x100 arrays
ax.plot_surface(x,y,z)
ax.set_xlabel('x_1')
ax.set_ylabel('x_2')
ax.set_zlabel('x_3')
p.show()

0.0000000000 0.0000000000 1.0000000000
0.1081190184 0.0000000000 0.9941379572
0.1075065026 0.0114923473 0.9941379572
0.0000000000 0.0000000000 1.0000000000


0.0000000000 0.0000000000 1.0000000000
0.1075065026 0.0114923473 0.9941379572
0.1056758951 0.0228544818 0.9941379572
0.0000000000 0.0000000000 1.0000000000


0.0000000000 0.0000000000 1.0000000000
0.1056758951 0.0228544818 0.9941379572
0.1026479375 0.0339576658 0.9941379572
0.0000000000 0.0000000000 1.0000000000
...

[x,y,z] = [[float(i) for i in line.split()] for line in open('surf')]

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

u = np.linspace(0, 2 * np.pi, 100)
v = np.linspace(0, np.pi, 100)

x = 10 * np.outer(np.cos(u), np.sin(v))
y = 10 * np.outer(np.sin(u), np.sin(v))
z = 10 * np.outer(np.ones(np.size(u)), np.cos(v))
ax.plot_surface(x, y, z,  rstride=4, cstride=4, color='b')

x = 10 * np.outer(np.cos(u), np.sin(v))
y = 20 * np.outer(np.sin(u), np.sin(v))
z = 5 * np.outer(np.ones(np.size(u)), np.cos(v))
ax.plot_surface(x, y, z,  rstride=4, cstride=4, color='b')

plt.show()

#!/usr/bin/octave -q
M=load('surf');
x=reshape(M(:,1),4,[]);
y=reshape(M(:,2),4,[]);
z=reshape(M(:,3),4,[]);
patch(x,y,z,'r');
view(3);
pause

M=load('surf');
x=M(1:end,1)';
y=M(1:end,2)';
z=M(1:end,3)';
tri = delaunay(x,y);
#trimesh(tri,x,y,z)
trisurf(tri,x,y,z)
#trimesh(tri,x,y,zeros(size(x)))
pause