matlabSkryptyWykresy

% zad 1.
x=-6:0.1:6;
y=x./(x.^2-4);
plot(x,y)
grid on

% zad 2.
x=0:0.1:5;
y=sqrt(x.^2 +1)./(x+1);
plot(x,y)
grid on

% zad 3.
x=0:0.1:2*pi;
y=((cos(x)).^2 + 1)./(x.^2 + 1).^2;
plot(x,y)
grid on

% zad 4.
x=0:0.001:2;
for k=1:5
    y=sin(k.*x*pi);
    subplot(5,3,k)
    plot(x,y)
end
grid on

% zad 5.
a=2;
b=4;
for x=1:0.01:10
    if x>b+a
        y=1;
    elseif (b-a<=x) && (x<=b+a)
        y=0.5*(1+sin(pi*((x-b)/(2*a))));
    elseif x<b-a
        y=0;
    end
    plot(x,y)
    hold on
end
grid on
hold off

% zad 6.
x=0:0.01:10*pi;
sum=0;
for k=1:10
    sum=sum+k*(sin(x)).^k;
end
plot(x,sum)
grid on
hold off

% zad 7.
vx=-1:0.01:1;
vy=vx;
% surf, meshc, surfc, surfl, contour, contour3 - do wytestowania
[x,y]=meshgrid(vx,vy);
z=cos(x.*y).*exp((x.^2+y.^2)./(10));
contour3(x,y,z)
grid on

% zad 8.
vx=-1:0.01:1;
vy=vx;
[x,y]=meshgrid(vx,vy);
z=abs(sqrt(x.^2+y.^2)-1.5);
mesh(x,y,z)
grid on

% zad 9.
hold on
vx=-2:0.01:2;
vy=vx;
[x,y]=meshgrid(vx,vy);
n = length(vx);
for i=1:n
    for j=1:n
        if vx(i)+vy(j)>=0
            z(i,j)=1./(1+exp(-(vx(i)+vy(j))));
        else
            z(i,j)=0.5;
        end
    end
end
mesh(x,y,z)
grid on
hold off

% zad 10.
hold on
x=-2:0.01:2;
plot(x,-x+1,x,x.^2-1,x,sin(x))
for x=-2:0.01:2
    for y=-2:0.01:2
        if (y<-x+1) && (y>x^2-1) && (y<sin(x))
            plot(x,y)
        end
    end
end
grid on
hold off

% zad 11.
hold on
for k=1:5
    x=-k:0.1:k;
    y=(k^2-x.^2).^(1/2);
    plot(x,y)
    y=-(k^2-x.^2).^(1/2);
    plot(x,y)
end
grid on
hold off