clc;
clear all;
[y1,Fs1] = audioread('FM_sinus.wav');

   Fs = 32000;
   dt = 1/Fs;
   time = 10;
   t = (0:dt:time);
   fm=0.5;
   f = 500;
   x = sin(2*pi*f*t);


    int_x = cumsum(x)/Fs;
   y = cos(2*pi*f*t + 2*pi*fm*int_x );
   y =diff(y);
   audiowrite('sinus22.wav',y, Fs);