MATLAB: Getting an error “INDEX OUT OF BOUNDS” in the matlab program. Please help.

fem matlab

for i=1:3 fepsn=i; end %%Loop for epsilon values for icounter=1:3 epsn=fepsn(icounter) Attempted to access fepsn(2); index out of bounds because numel(fepsn)=1.(ERROR HERE)
end
%%Shape function for axial effect N1(epsn)=(-1/2)*epsn*(1-epsn); N2(epsn)=(1+epsn)*(1-epsn); N3(epsn)=(1/2)*epsn*(1+epsn);
%%First Derivatives of Shape Function N1=diff(N1(epsn),epsn); N2=diff(N2(epsn),epsn); N3=diff(N3(epsn),epsn);
%%%%Second Derivatives of Shape Function N11=diff(N1(epsn),epsn,2); N22=diff(N2(epsn),epsn,2); N33=diff(N3(epsn),epsn,2);

Best Answer

for i=1:3
fepsn(i)=i; % (i) was added to avoid overwriting of the variable
end