function buding(a,b,c)
A=gongyuesu(a,b);
if rem(c,A)==0
    a=a/A;
    b=b/A;
    c=c/A;
    B=1;
    if a<0
        a=-a;
        b=-b;
        c=-c;
        if b<0
            b=-b;
            B=-1;
        end
    elseif b<0
           b=-b;
           B=-1;
    end
    k=1;
    if a<b
        
        t=a;
        a=b;
        b=t;
        d=b;
        a2=a;
        b2=b;
        r=rem(a,b);
        q(k)=fix(a/b);
        while r~=0
            a=b;
            b=r;
            d=r;
            r=rem(a,b);
            k=k+1;
            q(k)=fix(a/b);
        end
        p(1)=1;
        p(2)=q(1);
        Q(1)=0;
        Q(2)=1;
        for i=3:k
            p(i)=p(i-1)*q(i-1)+p(i-2);
            Q(i)=Q(i-1)*q(i-1)+Q(i-2);
        end
        a1=(-1)^k*Q(k)*c;
        b1=(-1)^(k-1)*p(k)*c;
        fprintf('x=%d+%d*t\n',b1,a2);
        fprintf('y=%d+(%d)*t',a1,-B*b2);
    else
        a2=a;
        b2=b;
        d=b;
        r=rem(a,b);
        q(k)=fix(a/b);
        while r~=0
            a=b;
            b=r;
            d=r;
            r=rem(a,b);
            k=k+1;
            q(k)=fix(a/b);
        end
        p(1)=1;
        p(2)=q(1);
        Q(1)=0;
        Q(2)=1;
        for i=3:k
            p(i)=p(i-1)*q(i-1)+p(i-2);
            Q(i)=Q(i-1)*q(i-1)+Q(i-2);
        end
        
        a1=(-1)^k*Q(k)*c;
        b1=(-1)^(k-1)*p(k)*c;
        fprintf('x=%d+(%d)*t\n',a1,-B*b2);
        fprintf('y=%d+%d*t',b1,a2);
    end
else
    disp('此不定方程無解')
end