intCellFit

 intCellFit : function used to fit to the shape of the cell.

 INPUT :
   x : value to be fit
   x1 : 
   x2 : 
   x3 : height or radius to be fit to
 OUTPUT :
   y : fitted output

 Copyright (C) 2016 Wiggins Lab 
 Written by Paul Wiggins.
 University of Washington, 2016
 This file is part of SuperSegger.
 
 SuperSegger 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 3 of the License, or
 (at your option) any later version.
 
 SuperSegger 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 for more details.
 
 You should have received a copy of the GNU General Public License
 along with SuperSegger.  If not, see <http://www.gnu.org/licenses/>.

0001 function y = intCellFit( x, x1, x2, h )
0002 % intCellFit : function used to fit to the shape of the cell.
0003 %
0004 % INPUT :
0005 %   x : value to be fit
0006 %   x1 :
0007 %   x2 :
0008 %   x3 : height or radius to be fit to
0009 % OUTPUT :
0010 %   y : fitted output
0011 %
0012 % Copyright (C) 2016 Wiggins Lab
0013 % Written by Paul Wiggins.
0014 % University of Washington, 2016
0015 % This file is part of SuperSegger.
0016 %
0017 % SuperSegger is free software: you can redistribute it and/or modify
0018 % it under the terms of the GNU General Public License as published by
0019 % the Free Software Foundation, either version 3 of the License, or
0020 % (at your option) any later version.
0021 %
0022 % SuperSegger is distributed in the hope that it will be useful,
0023 % but WITHOUT ANY WARRANTY; without even the implied warranty of
0024 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
0025 % GNU General Public License for more details.
0026 %
0027 % You should have received a copy of the GNU General Public License
0028 % along with SuperSegger.  If not, see <http://www.gnu.org/licenses/>.
0029 
0030 
0031 y = 0*x;
0032 
0033 y(and(x >= (h+x1), x <= (-h+x2))) = 2*h;
0034 ind = and(x >= x1, x < (h+x1) );
0035 
0036 y(ind) = 2*h*sqrt(1-((x(ind)-x1-h)/h).^2);
0037 
0038 ind = and(x <= x2, x > (-h+x2));
0039 y(ind) = 2*h*sqrt(1-((x(ind)-x2+h)/h).^2);
0040 
0041 
0042 end