x= [18     0.0000000000    5.00000000e+000    47.5878647956;
      23     0.0000000000    4.58758548e-001     3.0559374287;
      28     0.0000000000    4.20918810e-002     0.2175395217;
      33     0.0841837620    8.47058630e-002     0.4802086857;
      38     0.2535954881    1.02581530e-001     0.6539080493;
      43     0.4587585477    1.02581530e-001     0.7252064508;
      48     0.6639216073    8.47058630e-002     0.6450704424;
      53     0.8333333333    4.20918810e-002     0.3340042790;
      58     0.9175170954    9.23207464e-001     8.4659292536;
      63     0.9175170954    8.47058630e-002     0.6955059690;
      68     1.0869288214    1.70462404e-001     1.4756299733;
      73     1.4278536298    2.06435465e-001     1.8824475479;
      78     1.8407245589    2.06435465e-001     1.9486234930;
      83     2.2535954881    1.70462404e-001     1.6408337626;
      88     2.5945202964    8.47058630e-002     0.8228885321;
      93     2.7639320225    1.11803399e+000    11.0656169563;
      98     2.7639320225    1.02581530e-001     1.0016932048;
     103     2.9690950821    2.06435465e-001     2.0288180279;
     108     3.3819660113    2.50000000e-001     2.4728755195;
     113     3.8819660113    2.50000000e-001     2.4836423873;
     118     4.3819660113    2.06435465e-001     2.0558530065;
     123     4.7948369404    1.02581530e-001     1.0227348106;
     128     5.0000000000    1.11803399e+000    11.1682816273;
     133     7.2360679775    9.23207464e-001     9.2308613852;
     138     9.0824829046    4.58758548e-001     4.5874489930];
 xi = x(:,2);
 yi = customized_fun(xi);
 plot(xi,yi,'ro');
 hold on;
 fplot(@customized_fun,[0 10]);
 hold on;
 xi1 =x(:,2)+x(:,3);
 yi1 = customized_fun(xi1);
 plot(xi,yi,'b.');