`clc; clear all; close all;Bx = 720; By = 360; g = gcd(Bx,By);f = factor(g);f_prod = prod(f);`

In this example f = [ 2 2 2 3 3 5 ], and prod(f) = 360, but I just want to receive the product of distinct primes (i.e. 2*3*5). Is there a way to do this?

## Best Answer