How to take the product of distinct array values?

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?

Answer 1

[John D'Errico]

prod(unique(f))

[Sean]

Cheers mate.