分段函数

dy=diff(y);plot(x(2:end),dy,"b");xlim([0,86]);

创建一个f1.m函数：function y = f1(x)if rem(x,1)==0y = 3*factorial(abs(2*x));elsey = 3*4*sin(2*x)* + x^1.8;end创建主函数main.mx = [-4,3,14.37];Y = zeros(1,numel(x));for i = 1:numel(x)Y(i) = f1(x(i));endy = Y(1)+Y(2)-Y(3);两个文件需在同一文件夹下，运行main.m文件即可。