mardiamv25 function

Simulated data