gongbeisu.m

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· M 代码 · 共 20 行

function c=gongbeisu(a,b)
a1=a;
b1=b;
if a<b
    t=a;
    a=b;
    b=t;
end
c=b;
r=rem(a,b);
q=fix(a/b);
while r~=0
    a=b;
    b=r;
    c=r;
    r=rem(a,b);
    q=fix(a/b);
end
c=a1*b1/abs(c);