matlab code

for rap