diff --git a/Demo.m b/Demo.m
index 63640af522c4d4dc9b458836ae034494d0454c5c..33118783767f013632723995310dd1d1ff8cb091 100644
--- a/Demo.m
+++ b/Demo.m
@@ -1,3 +1,5 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
+%% Demo
close all
clear
diff --git a/Export2inp.m b/Export2inp.m
index a24a2763a3e9ff77481e2179faab2ed7b1c86936..bbccd7fc496e6bee7f8cd0bcabf9e49155989ade 100644
--- a/Export2inp.m
+++ b/Export2inp.m
@@ -1,5 +1,5 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% export points to inp
-
samples = {'out_out_CD','out_out_MD','in_in_CD','in_in_MD'};
degree = {'0','45','90','180'};
diff --git a/Functions/functions_auxiliary/ginputWhite.m b/Functions/functions_auxiliary/ginputWhite.m
index c9c1b220a25f861a2aac1791e128a814eb429343..1b9a64c81579596bc6600766f271ac729e61e1b6 100644
--- a/Functions/functions_auxiliary/ginputWhite.m
+++ b/Functions/functions_auxiliary/ginputWhite.m
@@ -1,4 +1,5 @@
function [out1,out2,out3] = ginputWhite(arg1)
+% The code is modified by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%GINPUT Graphical input from mouse.
% [X,Y] = GINPUT(N) gets N points from the current axes and returns
% the X- and Y-coordinates in length N vectors X and Y. The cursor
diff --git a/Functions/functions_cost/smoothsurface_cost.m b/Functions/functions_cost/smoothsurface_cost.m
index eb93a6736624c64e9af9d86b8d0d53564461e4ff..6b30586aa8e2b2950fb1bf50bb191f7c0dfc7329 100644
--- a/Functions/functions_cost/smoothsurface_cost.m
+++ b/Functions/functions_cost/smoothsurface_cost.m
@@ -1,4 +1,5 @@
function cost = smoothsurface_cost(V,mu,std,s1,s2)
+% The code is modifed by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
% abda@dtu.dk
if nargin == 1
diff --git a/Functions/functions_mesh/elements_mesh.m b/Functions/functions_mesh/elements_mesh.m
index f62f26a20d0c53f42c160fd2bf29251b8c16900e..858420a2a76c02a2b4501a1301543afc80ce2f82 100644
--- a/Functions/functions_mesh/elements_mesh.m
+++ b/Functions/functions_mesh/elements_mesh.m
@@ -1,5 +1,6 @@
function ele = elements_mesh(mesh)
-
+% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version
+% 2020.04.10
sz0 = size(mesh);
ny = sz0(1);
nz = sz0(2);
diff --git a/Functions/functions_mesh/write_mesh.m b/Functions/functions_mesh/write_mesh.m
index 9495f455914a59a7cdee8f15784c05d714fa7c90..9f9b915cfef9ddbe9726629d814a4f61b62bca6a 100644
--- a/Functions/functions_mesh/write_mesh.m
+++ b/Functions/functions_mesh/write_mesh.m
@@ -1,5 +1,6 @@
function write_mesh(nodes,elements,filename)
-
+% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version
+% 2020.04.10
fid = fopen([filename,'.inp'],'w');
[~,name,~] = fileparts(filename);
fprintf(fid,'*Part, name=%s \n',name);
diff --git a/Functions/functions_visualization/suptitle_modified.m b/Functions/functions_visualization/suptitle_modified.m
index 3fafe385b8903535b811ec06fa7cd647510c3af7..b1e93e65d4e1f523732a3af3a2fcb813f73db771 100644
--- a/Functions/functions_visualization/suptitle_modified.m
+++ b/Functions/functions_visualization/suptitle_modified.m
@@ -1,4 +1,6 @@
function hout=suptitle_modified(str,interpreter)
+% The code is modified by Behnaz Pirzamanbein, bepi@dtu.dk last version
+% 2020.04.10
%SUPTITLE puts a title above all subplots.
%
% SUPTITLE('text') adds text to the top of the figure
diff --git a/LayerDetection.m b/LayerDetection.m
index 81252ca57cce20e8d22b08c552a60c8358323c36..54e3dbfd0adfa570de0f2cb31910d07177fa5dde 100644
--- a/LayerDetection.m
+++ b/LayerDetection.m
@@ -1,3 +1,4 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% detecting layers
% create result folder
save_folder = 'Results\layer\';
diff --git a/Measure.m b/Measure.m
index a7d0e357d8034803db88cddc9fd632d789df39e3..81b5bb75075d90f62802a9cdc8107420f5bd70d8 100644
--- a/Measure.m
+++ b/Measure.m
@@ -1,3 +1,4 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% fix the meaure distance plot same xlim and ylim and also caxes
samples = {'out_out_CD','out_out_MD','in_in_CD','in_in_MD'};
degree = {'0','45','90','180'};
diff --git a/PointSelection.m b/PointSelection.m
index 7eb0e6fc980b731f18b7e1a4e7846dbc4ea80864..18fc38aab6e65d8f8a94724842f311bb492d9188 100644
--- a/PointSelection.m
+++ b/PointSelection.m
@@ -1,5 +1,5 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% move from pipeline to here to select the points
-
save_folder = 'Results\point\';
if ~exist(save_folder,'dir')
mkdir(save_folder)
diff --git a/PrepMeasure.m b/PrepMeasure.m
index dbd7bad1e48ad2c54f8ae030713f2378fc93f6c3..801e9bf4465d825921642e551f708a918b510691 100644
--- a/PrepMeasure.m
+++ b/PrepMeasure.m
@@ -1,3 +1,4 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% measuring the distance betwee two sruface and length of each surface
addpath(genpath('.\Functions'))
diff --git a/PrepRegistration.m b/PrepRegistration.m
index 0ae1eeb3f3572c52eaeafd6c466bfc27b7f8dd8a..1d259a326a95efcfd910fb71a638fc1e531605b6 100644
--- a/PrepRegistration.m
+++ b/PrepRegistration.m
@@ -1,3 +1,4 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% summing over 10 pixels along normal direction of the smoothed surface
addpath(genpath('.\Functions'))
diff --git a/Registration.m b/Registration.m
index 7460ebd479aac485f2c1a5a6fb3e47762cd0c603..b0b043a4a4e2b6e35483b8dbd3b9e116acb2ac1f 100644
--- a/Registration.m
+++ b/Registration.m
@@ -1,3 +1,4 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
%% registration
samples = {'out_out_CD','out_out_MD','in_in_CD','in_in_MD'};
diff --git a/Transformation.m b/Transformation.m
index 0779d1c5f018665a510435f2ad04311c1895f2e0..5d1da56d970f39bfad6628bf33f2df1e72eee4b6 100644
--- a/Transformation.m
+++ b/Transformation.m
@@ -1,3 +1,5 @@
+%% The code is written by Behnaz Pirzamanbein, bepi@dtu.dk last version 2020.08.04
+%% Transformation
addpath(genpath('.\Functions'))
samples = {'out_out_CD','out_out_MD','in_in_CD','in_in_MD'};