mx=[...
   (1-vrs(4))*exp(vrs(1))    vrs(4)*vrs(5)       vrs(4)*vrs(5);
    vrs(4)*vrs(5)       (1-vrs(4))*exp(vrs(2))   vrs(4)*vrs(5);
    vrs(4)*vrs(5)       vrs(4)*vrs(5)       (1-vrs(4))*exp(vrs(3))];