Segundo problema del proyecto Euler - "Resuelto"

El segundo problema dio algo más de trabajo, y teniendo en cuenta en por lo menos en latinoamérica la mayoría de los inscritos no han pasado del segundo, esta bien que ya lo haya resuelto.

El problema radica en codificar un script para que tome todos los términos pares no mayores a 4000000, en la serie de Fibonacci y los sume, a continuación sito el código con el cual resolví esa pregunta.

clc;
f(1)=1;f(2)=1;RTA=0;
n=input('Indique el número de elementos a calcular de la serie de Fibonacci ');
disp(f(1));disp(f(2));
for i=3:n
      f(i)=f(i-1)+f(i-2);
      if f(i)<4000000 && rem(f(i),2)==0 %suma todos los términos pares no mayores a 4000000
          RTA=RTA+f(i); 
      end 
      disp(f(i)) %muestra los elementos de la serie end disp('suma de la serie') 
end
disp('La suma de los términos pares no mayores a 4000000 es')  
disp(RTA) %muestra el resultado de la suma