function dekartovlist
title(‘Декартов лист’)
hold on
t=-1000:0.01:1000;
plot(3*t.*(1+t.^3).^(-1),3*t.^2.*(1+t.^3).^(-1),’b’)
axis([-3 3 -3 3])
plot(t,-t-1,’r’)
legend(‘кривата’, ‘асимптотата’)
grid on
end

function graphcos
x=-pi:0.001:pi;
plot(x,x.^2.*(2+cos(x.^(-1))),’r’)

function sigmaproc
x=-1:0.01:1;
y=-1:0.01:1;
[X,Y]=meshgrid(x,y);
Z=X.*Y.*(X.^2+Y.^2).^(-1);
surfl(X,Y,Z), shading interp
colormap(‘gray’)
end

function vitline
hold on
t=0:0.1:4*pi; x=-10:0.1:10;
y=-10:0.1:10; [X,Y]=meshgrid(x,y);
plot3(t.*cos(2t),t.*sin(2*t),t)
surfl(X,Y,sqrt(X.^2+Y.^2))
shading interp
grid on
end

16.10

function sequence
nrf=400;
grid on
hold on
axis([-2 2 -2 2])
function anime
nrf=100;
hold on
axis([-1 2 -1 2])
for i=1:nrf/2
fimplicit(@(x,y) x.^3+y.^3-3x.y-(nrf/2-i)/10,[-2 2 -2 2])
getframe;
end for i=nrf/2:nrf
fimplicit(@(x,y) x.^3+y.^3-3x.*y+(i-nrf/2)/10,[-2 2 -2 2])
getframe;
end

function anime2
t=-100:0.1:100;
x=3*t./(1+t.^3); y=3*t.^2./(1+t.^3);
grid on
axis([-2 2 -2 2])
for k=1:length(t)
xx(k)=x(k); yy(k)=y(k);
plot(xx,yy,’g’);
getframe;
end
end