switchtodesign;
clear;

rise_period = linspace(1e-100, 0.1, 10);
fall_period = rise_period;
tdecq_result = matrix(length(rise_period), 1);

for (ii = 1:length(rise_period)) {
    switchtodesign;
    setnamed("::Root Element" ,"rise period", rise_period(ii));  
    setnamed("::Root Element", "fall period", fall_period(ii));   
    
    run;
    runanalysis;
    
    tdecq_result(ii) = getresult("TDECQ_1", "TDECQ");                
}

plot(rise_period, tdecq_result, "rise oeriod (a.u.)", "TDECQ (dB)", "TDECQ vs rise period");