I was rereading Terence Tao's Analysis 1 and found this question in the section:
Exercise $2.3.4.$ Prove the identity $(a + b)^2 = a^2 + 2ab + b^2$ for all natural numbers a, b.
Prior to this we already have proved:
$1.a\cdot b=b\cdot a\\2.a\cdot b=0\implies a=0\lor b=0\\3.a\cdot (b\cdot c)=(a\cdot b)\cdot c\\4. a\cdot (b+c)=a\cdot b+a\cdot c\\$
So I wrote this down:
$(x+y)^2\\=(x+y)(x+y)\\=x(x+y)+y(x+y)\\=x^2+xy+xy+y^2\\=x^2+2xy+y^2$
Is this a valid proof?
Here is the modified version of your proof with few skipped steps added and properties mentioned.
\begin{align} (x+y)^2 &=(x+y).(x+y) & \\ &=(x(x+y))+(y(x+y))&\text{(using property 4 and 1)} \\ &=(x.x +x.y )+(y.x +y.y) &\text{(using property 4 again)}\\&=x^2+(x.y+y.x)+y^2 &\text{(using associativity of addition)}\\ &=x^2+(x.y+x.y)+y^2 &\text{(since $y.x=x.y$, property 1)}\\&=x^2+(1.xy+1.xy)+ y^2 &\text{(using property: $1.a=a$ for all natural numbers $a$)} \\ &=x^2 + (1+1). xy + y^2 &\text{(using property 4)}\\ &=x^2+2xy+y^2 &\text{(using 1+1=2)}\end{align}