Табулирование функций y(x) на отрезке [a,b] с последующим построением графика

eps=0.0000001; 
xn=-5*pi+eps; 
xk=10-eps; 
n=100; 
x=xn: (xk-xn)/n:xk; 
Nk=length(x); 
f=zeros(Nk,1); 
for k=1:Nk 
if x(k)<-3*pi 
f(k)=(cos(x(k))+sin(2*x(k))); 
if -3*pi<=x(k)<0
    f(k)=(cos(x(k))+cos(2*x(k)));
if 0<=x(k)<1
    f(k)=((1-2*x(k))^3);
    if x(k)>=1
        f(k)=(abs(abs(x(k)-6)-4)-2);
else 
f(k)=(-1); 
end 
end 
end
end
end
[ymin, kmin]=min(f); 
[ymax, kmax]=max(f); 
plot(x,f,'b',x(kmin),f(kmin),'rx',x(kmax),f(kmax),'rx'); 
text(x(kmin),f(kmin),'min'); 
text(x(kmax),f(kmax),'max'); 
fmax=max(f); 
fmin=min(f);

clear all;
clc;
close all;
xn=-5*pi; 
xk=10; 
n=100; 
x=xn: (xk-xn)/n:xk; 
Nk=length(x); 
f=zeros(Nk,1); 
for k=1:Nk 
if x(k)<-3*pi 
    f(k)=(cos(x(k))+sin(2*x(k))); 
    else
        if -3*pi<=x(k) && x(k)<0
            f(k)=(cos(x(k))+cos(2*x(k)));
        else
            if 0<=x(k) && x(k)<1
                f(k)=((1-2*x(k))^3);
            else
                if x(k)>=1
                    f(k)=(abs(abs(x(k)-6)-4)-2);
                else 
                    f(k)=(-1); 
                end 
            end 
        end
    end
end
[ymin, kmin]=min(f); 
[ymax, kmax]=max(f); 
plot(x,f,'b',x(kmin),f(kmin),'rx',x(kmax),f(kmax),'rx'); 
text(x(kmin),f(kmin),'min'); 
text(x(kmax),f(kmax),'max'); 
fmax=max(f); 
fmin=min(f);

clear all; 
clc; 
close all; 
xn=-3*pi; 
xk=3*pi; 
n=100; 
x=xn:  (xk-xn)/n:xk; 
Nk=length(x); 
f=zeros(Nk,1); 
for k=1:Nk 
if x(k)>=0 
f(k)=(sin(sqrt(x(k)))); 
else 
if x(k)<0 
f(k)=(cos(x(k)^1/3)); 
else 
f(k)=(-1); 
end 
end 
end 
[ymin, kmin]=min(f); 
[ymax, kmax]=max(f); 
plot(x,f,'b',x(kmin),f(kmin),'rx',x(kmax),f(kmax),'rx');grid on 
text(x(kmin),f(kmin),'min'); 
text(x(kmax),f(kmax),'max'); 
fmax=max(f); 
fmin=min(f); 
[max_val,max_ind] = max(f); 
fprintf('ymax = %f, Xmax = %f\n',max_val,x(max_ind)); 
[min_val,min_ind] = min(f); 
fprintf('Ymin = %f, Xmin = %f\n',min_val,x(min_ind));