compute circular using linear matlab code