SF1672forelasn20220120_26matlab.m