code

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