diff --git a/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.cpp b/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.cpp
new file mode 100644
index 00000000..3c4855c5
--- /dev/null
+++ b/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.cpp
@@ -0,0 +1,99 @@
+/****************************************************************************
+* VCGLib                                                            o o     *
+* Visual and Computer Graphics Library                            o     o   *
+*                                                                _   O  _   *
+* Copyright(C) 2004-2016                                           \/)\/    *
+* Visual Computing Lab                                            /\/|      *
+* ISTI - Italian National Research Council                           |      *
+*                                                                    \      *
+* All rights reserved.                                                      *
+*                                                                           *
+* This program is free software; you can redistribute it and/or modify      *
+* it under the terms of the GNU General Public License as published by      *
+* the Free Software Foundation; either version 2 of the License, or         *
+* (at your option) any later version.                                       *
+*                                                                           *
+* This program is distributed in the hope that it will be useful,           *
+* but WITHOUT ANY WARRANTY; without even the implied warranty of            *
+* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the             *
+* GNU General Public License (http://www.gnu.org/licenses/gpl.txt)          *
+* for more details.                                                         *
+*                                                                           *
+****************************************************************************/
+/*! \file trimesh_intersection_plane.cpp
+\ingroup code_sample
+
+\brief An example of computing the intersection of a mesh with a plane, 
+saving this polyline as a well projected 2D SVG 
+and splicing the mesh in the two components below and over the plane
+
+*/
+
+#include <vcg/complex/complex.h>
+#include <vcg/complex/algorithms/clean.h>
+#include <vcg/complex/algorithms/intersection.h>
+#include <vcg/complex/algorithms/refine.h>
+#include <vcg/complex/algorithms/Update/color.h>
+#include <vcg/complex/algorithms/create/platonic.h>
+#include<vcg/complex/algorithms/isotropic_remeshing.h>
+
+#include <wrap/io_trimesh/export_off.h>
+
+using namespace std;
+using namespace vcg;
+
+class MyFace;
+class MyVertex;
+
+struct MyUsedTypes : public UsedTypes<	
+    Use<MyVertex>::AsVertexType,
+    Use<MyFace>  ::AsFaceType>{};
+
+
+class MyVertex  : public Vertex< MyUsedTypes, vertex::Coord3f, vertex::VFAdj, vertex::BitFlags, vertex::Normal3f, vertex::Qualityf, vertex::Mark>{};
+class MyFace    : public Face  < MyUsedTypes, face::VertexRef, face::VFAdj, face::FFAdj, face::Color4b, face::BitFlags, face::Mark, face::Normal3f> {};
+
+class MyMesh:     public tri::TriMesh< vector<MyVertex>, vector<MyFace > >{};
+
+int main(int ,char **)
+{  
+  MyMesh m1, m2;
+  Sphere(m1);
+  Hexahedron(m2);
+  printf("Created a sphere mesh of %i vertices %i edges\n",m1.VN(),m1.FN());
+  printf("Created a cube mesh of %i vertices %i edges\n",m2.VN(),m2.FN());
+  
+  math::MarsenneTwisterRNG rnd;
+  Point3f direction = vcg::math::GeneratePointOnUnitSphereUniform<float,math::MarsenneTwisterRNG>(rnd);
+  tri::UpdatePosition<MyMesh>::Translate(m1,direction);
+  for(int i=0;i<10;++i)
+  {
+  tri::Clean<MyMesh>::SelectIntersectingFaces(m1,m2);
+  tri::UpdateSelection<MyMesh>::FaceDilate(m1);
+  tri::Clean<MyMesh>::SelectIntersectingFaces(m2,m1);
+  tri::UpdateSelection<MyMesh>::FaceDilate(m2);
+  IsotropicRemeshing<MyMesh>::Params params;
+
+  float len = (tri::Stat<MyMesh>::ComputeFaceEdgeLengthAverage(m1,true) + tri::Stat<MyMesh>::ComputeFaceEdgeLengthAverage(m1,true));
+  params.SetTargetLen(len);
+  params.SetFeatureAngleDeg(10);
+  params.iter=3;
+  params.selectedOnly=true;
+  printf(" Input mesh %8i v %8i f\n",m1.VN(),m1.FN());
+  IsotropicRemeshing<MyMesh>::Do(m1, params);
+  IsotropicRemeshing<MyMesh>::Do(m2, params);
+  printf(" Input mesh %8i v %8i f\n",m1.VN(),m1.FN());
+  }
+  tri::Clean<MyMesh>::SelectIntersectingFaces(m1,m2);
+  tri::Clean<MyMesh>::SelectIntersectingFaces(m2,m1);  
+  int selCnt= tri::UpdateColor<MyMesh>::PerFaceConstant(m1,Color4b::Red,true);
+  printf("Intersected %i faces on sphere\n",selCnt);
+  selCnt= tri::UpdateColor<MyMesh>::PerFaceConstant(m2,Color4b::Red,true);
+  printf("Intersected %i faces on cube\n",selCnt);
+  
+  tri::io::ExporterOFF<MyMesh>::Save(m1,"sphere.off",tri::io::Mask::IOM_FACECOLOR);
+  tri::io::ExporterOFF<MyMesh>::Save(m2,"cube.off",tri::io::Mask::IOM_FACECOLOR);
+  
+  return 0;
+}
+
diff --git a/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.pro b/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.pro
new file mode 100644
index 00000000..9d676c3c
--- /dev/null
+++ b/apps/sample/trimesh_intersection_mesh/trimesh_intersection_mesh.pro
@@ -0,0 +1,3 @@
+include(../common.pri)
+TARGET = trimesh_intersection_mesh
+SOURCES += trimesh_intersection_mesh.cpp