float temp_x = vertices[i].x;

         float temp_y = vertices[i].y;



         vertices[i].x = tstep*(vertices[i].vx * vertices[i].mass + vertices[i].fx)/vertices[i].mass + temp_x;

         vertices[i].y = tstep*(vertices[i].vy * vertices[i].mass + vertices[i].fy)/vertices[i].mass + temp_y;



         vertices[i].vx = (vertices[i].x - temp_x)/tstep;

         vertices[i].vy = (vertices[i].y - temp_y)/tstep;