SCRIPT S219
function A = s219(a,h,b)
f=@(x) x.*(x+3.*(x+3.*(x+3.*(x+3).^2).^2).^2).^2;
A=[];
j=1;
for i=a:h:b
    A(j,1)=j;
    A(j,2)=i;
    A(j,3)=f(i);
    j=j+1;
end
end
COD S219
a=input('Inceputul intervalului:');
b=input('Sfarsitul intervalului:');
h=input('Pasul intervalului:');
A=s219(a,h,b);
plot(A(:,3),A(:,2),'-d','Color',[0.1 0.3 0],'LineWidth',1)
disp('Matricea este:');
disp(A)