%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.