% script for doing problem 5 on practice quiz 1
x=[1:sqrt(10000000)].^2;
y=x(find(mod(x,10)==1));
sum(y)

