need a matlab(R2015a) code for Generating Kuramoto model

I have time series data(amplitude of 2000 values) of two weakly interacting coupled oscillators.could you pleae help me writing a code to generate Kuramoto model to show wheter they are in synchronization or not on unit radius circle.I have stored data of two oscillators separately in dat files which I can't upload here so If you give me a pseudo code I can modify for my data.