SF1672forel20220131matlab.m