code: def check_angle_solution(ang,n_vec,toy_aa,timeseries):
      
      f,a=plt.subplots(3,1)
      for i in range(3):
          a[i].plot(toy_aa.T[i+3],'.')
          size = len(ang[6:])/3
          AA = np.array([np.sum(ang[6+i*size:6+(i+1)*size]*np.sin(np.sum(n_vec*K,axis=1))) for K in toy_aa.T[3:].T])
          a[i].plot((ang[i]+ang[i+3]*timeseries-2.*AA) % (2.*np.pi),'.')
          a[i].set_ylabel(r'$\theta$'+str(i+1))
      a[2].set_xlabel(r'$t$')
      plt.show()
docs: Plots the toy angle solution against the toy angles ---
          Takes true angles and frequencies ang,
          the Fourier vectors n_vec,
          the toy action-angles toy_aa
          and the timeseries