function [p r] = longdivision(m,n)
disp(['m=' num2str(m) ', n=' num2str(n)])
p=0;
r=m;
while r>=n
    p=p+1;
    r=m-p*n;
end
disp(['m=' num2str(p) 'n+' num2str(r)])
