summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Kukannozoku.sci
blob: 0e38dae789e550c83da03825d992541e978ff755 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
// 08.07.13
// 08.10.16

function Res=Kukannozoku(Jokyo,KukanL)
  Eps=10^(-6);
  N=size(KukanL,1);
  T1=Jokyo(1); T2=Jokyo(2);
  Tmp=KukanL(1,:);
  T1=max(T1,Tmp(1));
  Tmp=KukanL(N,:);
  T2=min(T2,Tmp(2));
  Res=[];
  Flg=0;
  for I=1:N
    Ku=KukanL(I,:);
    if Flg==0
      if Ku(2)<T1
        Res=[Res;Ku];
      else
        Flg=1;
        if Ku(1)<T1-Eps 
          Tmp=[Ku(1),T1];
          Res=[Res;Tmp];
        end;
        if Ku(2)>T2+Eps
          Tmp=[T2,Ku(2)];
          Res=[Res;Tmp];
        end;
      end;
    elseif Flg==1
      if Ku(2)<T2
        continue;
      else
        Flg=2;
        if Ku(1)<T2-Eps
          Ku=[T2,Ku(2)];
        end;
        Res=[Res;Ku];
      end;
    else
      Res=[Res;Ku];
    end;
  end;
endfunction;