There were two syntax errors:
- Extra bracket at the end of the last line.
- In one of the terms, you forgot to index into phi(i), using instead just "phi".
More importantly, there is a much better way to code this, using cell arrays. (Notice, in my code below, the curly brackets indexing into the variable "stiffness".) Search the documentation for "cell array" for details. You'll be glad you did.
elements = 3;
nodes = 2;
phi = [pi / 3,pi-pi/3,0];
for i = 1:elements
stiffness{i}=zeros(2)
stiffness{i}=[cos(phi(i))^2,sin(phi(i)) * cos(phi(i));sin(phi(i)) * cos(phi(i)),sin(phi(i))^2]'
end
There is also a good explanation here:
