Here is one way:
T = readtable('Stress_Reactivity.xlsx');
figure
gscatter(T.cort_3,T.diastolic_3,T.condition)
I didn't know exactly what you meant by "scatterplot matrix out of these 4 variables", so I just did two variables. I assumed the binary variable condition was how you wanted to group them.
