summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Symbdata.sci
blob: d63a34708e7ea9cf81e4c3457e6dd1016c370c5b (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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
// s : 2014.08.25
// 14.09.07

function Out=Symbdata(Moji)
    ML=list(); Sp=[]; CC=[];
    if Moji=="x" then
        M1=[-0.042,0.035;-0.027,0.045;-0.01,0.045;0.004,0.03;0.00,0.00;...
           -0.004,-0.03;0.01,-0.045;0.027,-0.045;0.042,-0.035];
        M2=[0.055,0.035;0.043,0.045;0.027,0.045;0.01,0.03;0.00,0.00;...
           -0.01,-0.03;-0.027,-0.045;-0.043,-0.045;-0.055,-0.035];
        ML=list(M1,M2);
        Sp=[100,100];
        CC=[0,0];
    end;
    if Moji=="y" then
        M1=[-0.06,0.06;...
            -0.045,0.080;...
            -0.03,0.080;...
            -0.022,0.065;...
            -0.032,0.035;...
            -0.035,-0.01;...
             0.007,-0.008;...
             0.037,0.04;...
             0.05,0.09];
        M2=[0.05,0.09;...
            0.037,0.04;...
            0.025,-0.00;...
            0.015,-0.035;...
            -0.000,-0.055;...
            -0.025,-0.062;...
            -0.05,-0.045];
        ML=list(M1,M2);
        Sp=[100,100];
        CC=[0,0];
    end
    if Moji=="z" then
        M1=[-0.035,0.035;-0.023,0.045;-0.01,0.045;...
            0.01,0.04;0.035,0.04;0.042,0.045];
        M2=[0.042,0.045;0.00,0.00;-0.042,-0.045];
        M3=[-0.042,-0.045;-0.035,-0.04;-0.01,-0.04;...
            0.01,-0.045;0.023,-0.045;0.035,-0.035];
        ML=list(M1,M2,M3);
        Sp=[100,0,100];
        CC=[0,0,0];
    end
    Out1=list()
    Out2=list()
    for J=1:length(ML)
        M=ML(J);
        N=Sp(J);
        for K=1:size(M,1)
            P=M(K,:)
            Out1($+1)=P
        end
        if N>0 then
            if CC(J)==0 then
                S=Splinedata(M,"N="+string(N))
            else
                S=Splinedata(M,"N="+string(N),"C")
            end
        else
            S=M
        end;
        Out2($+1)=S
    end
    Out=list(Out1,Out2)
endfunction