x=0;
y=0;
tick=5;
temp=0;
lim=101;
plot3(0,0,0,'.');
hold on
while temp<lim
    x=0;
    while x<temp
        y=temp-x;
        z=((x+y)^2+3*x+y)/2;
        plot3(x,y,z,'.');
        x=x+1;
    end
    temp=temp+tick;
end
xlabel('x')
ylabel('y')
zlabel('z')
