function[Area] = find_AreaCell(Cell_Number,X_Coordinate,Y_Coordinate,subtract_area)

num_cells = size(Cell_Number,1);
num_nodes = size(X_Coordinate,1);
time_steps = size(X_Coordinate,2);

for i=1:1:time_steps
    for j=1:1:num_cells
        mat1_node1 = Cell_Number(j,1);
        mat1_node2 = Cell_Number(j,2);
        mat1_node3 = Cell_Number(j,3);
        
        mat2_node1 = Cell_Number(j,1);
        mat2_node2 = Cell_Number(j,3);
        mat2_node3 = Cell_Number(j,4);
        
        mat3_node1 = Cell_Number(j,1);
        mat3_node2 = Cell_Number(j,6);
        mat3_node3 = Cell_Number(j,4);
        
        mat4_node1 = Cell_Number(j,4);
        mat4_node2 = Cell_Number(j,5);
        mat4_node3 = Cell_Number(j,6);
        
        Mat1 = [1 X_Coordinate(mat1_node1,i) Y_Coordinate(mat1_node1,i);...
            1 X_Coordinate(mat1_node2,i) Y_Coordinate(mat1_node2,i);...
            1 X_Coordinate(mat1_node3,i) Y_Coordinate(mat1_node3,i)];
        
        Mat2 = [1 X_Coordinate(mat2_node1,i) Y_Coordinate(mat2_node1,i);...
            1 X_Coordinate(mat2_node2,i) Y_Coordinate(mat2_node2,i);...
            1 X_Coordinate(mat2_node3,i) Y_Coordinate(mat2_node3,i)];
        
        Mat3 = [1 X_Coordinate(mat3_node1,i) Y_Coordinate(mat3_node1,i);...
            1 X_Coordinate(mat3_node2,i) Y_Coordinate(mat3_node2,i);...
            1 X_Coordinate(mat3_node3,i) Y_Coordinate(mat3_node3,i)];
        
        Mat4 = [1 X_Coordinate(mat4_node1,i) Y_Coordinate(mat4_node1,i);...
            1 X_Coordinate(mat4_node2,i) Y_Coordinate(mat4_node2,i);...
            1 X_Coordinate(mat4_node3,i) Y_Coordinate(mat4_node3,i)];
        
        Area(j,i) = (0.5*(abs(det(Mat1)) + abs(det(Mat2)) + abs(det(Mat3)) + abs(det(Mat4))) - subtract_area)*10^-9;  
    end
end