function i = mult_inv(a,p)
for(x=1:p)
m = mod(a*x,p);
 if (m==1) i=x; 
 endif;
endfor;
endfunction;