[Theta,tri]=cubesurf([1 10],50);
Pos=forward(Theta,[10 10 10 10]);
trisurf(tri,Pos(:,1),Pos(:,2),Pos(:,3))