for i = 1:1001,
    x(i)=-6 + (0.01*(i-1));
    y(i)= fungsi(x(i));
end
plot(x,y), grid on, xlabel('x'), ylabel('y = f(x)')