magsint.m

%magsint.m
function g=magsint（fnameabnm）
switch m
    case 1
        t=0;
        A=1;
    case 2
        t=[-1/sqrt（3）1/sqrt（3）];
        A=[11];
    case 3
        t=[-sqrt（0.6）0.0sqrt（0.6）];
        A=[5/98/95/9];
    case 4
        t=[-0.86113-0.3399810.3399810.861136];
        A=[0.3478550.