TechnikaObliczeniowa

//a)
function [y]=fun(x)
    if(x<-3)
        y=-6;
    elseif(-3<=x & x<=2)
        y=2*x;
    else y=x^2;
        end
endfunction


x=linspace(-5,5,50);
y = [];

for i = 1:length(x)
    y(i) = fun(x(i));
end
plot2d(x,y);

//b)
x=linspace(-10,10,2000);
y=abs(x);
plot2d(x,y);
plot2d(x,y);

//c)
x=linspace(-10,5,1500);
f1=poly([1,-2,0,3],'x','c');
//f1=x^3-2*x^2+3;
f2=derivat(f1);
f3=derivat(f2);

subplot(1,3,1);
plot2d(x,horner(f1,x),-4);
title("wykres f(x)");
xlabel("x");
ylabel("f(x)");
subplot(1,3,2);
plot2d(x,horner(f2,x),-4);
title("wykres f''(x)");
xlabel("x");
ylabel("f(x)");
subplot(1,3,3);
plot2d(x,horner(f3,x),-4);
title("wykres f''''(x)");
xlabel("x");
ylabel("f(x)");

//d)
r=5;
t=linspace(-%pi,%pi,100);

x=r*cos(t).*(1-2.*sin(t)^2);
y=r*sin(t).*(1-2.*cos(t)^2);

subplot(2,1,1);
plot2d(t,[x' y'],[2,3],leg="x(t)@y(t)");
subplot(2,1,2);
r=10;
x=r*cos(t).*(1-2.*sin(t)^2);
y=r*sin(t).*(1-2.*cos(t)^2);
plot2d(t,[x' y'],[2,3],leg="x(t)@y(t)");

//e)
t= -3*%pi:.01:3*%pi;
a=3;
polarplot(t,4*a.*cos(t./3));

//f)
x=-4:4;
y=-4:4;
z=zeros(length(x),length(y));

for i=1:length(x)
    for j=1:length(y)
        z(i,j)=sin(sqrt(x(i)^2+y(j)^2));
        end
end
plot3d1(x,y,z);