Problema al modelar un filtro de paso bajo en MATLAB

Estoy intentando ejecutar una simulación en MATLAB para modelar dos filtros de paso bajo. Es un filtro de paso bajo activo que se muestra en el siguiente esquema:

^{simular este circuito : esquema creado usando CircuitLab}

La única parte en la que estoy enfocado ahora es el circuito de retroalimentación con el condensador y la resistencia. Estoy tratando de ver cómo cambiar el condensador ajustaría el filtro de paso bajo mientras mantiene la resistencia igual.

Solo para esta parte, he logrado obtener la siguiente ecuación para el filtro de paso bajo, y lo simplifiqué enfocándome en la parte del filtro, asumiendo que se puede usar R como parte de la ganancia proporcionada por el amplificador:

\ $ \ frac {R} {1 + j * 2 * pi * f * R * C} \ $ --- > \ $ \ frac {1} {1 + j. * 2. * pi. * f. * R. * C} \ $

Usando el modelo simplificado, utilicé MATLAB para escribir diagramas para los diagramas de magnitud y bode de los filtros con sus respectivas capacidades, pero algo parece extraño con el gráfico de bode. La frecuencia en el punto de -3 dB no coincide con el punto 1 / sqrt (2) en el gráfico de magnitud. Sé que esta es una pregunta muy simple, pero parece que no puedo entender por qué la trama de Bode no coincide. ¿Podría haber algo que no esté tomando en cuenta? Mi código se muestra a continuación:

clc; clear all; close all;  f = 1:1000:(1e9);

Ra = 470;             Ca = 3.*(10.^-12);
AA = 1./(1+j.*2.*pi.*f.*Ra.*Ca);
magAA = abs(AA);      bodeAA = 0-20.*log(magAB);

Rb = 470;               Cb = 6.*(10.^-12);
AB = 1./(1+j.*2.*pi.*f.*Rb.*Cb);
magAB = abs(AB);        bodeAB = 0-20.*log(magAB);

%Graph Lines for -3 dB and 1/sqrt(2)
line3 = ones(size(f)); line3 = line3 .*-3;   
line2 = ones(size(f)); line2 = line2 .*(1./sqrt(2));

figure(1); %Plot the filter Magnitude
plot(f,magAA);      title('MagPlot 3pF');   hold on; plot(f,line2,'--r');

figure(2); %Plot the filter magnitude
plot(f,magAB);      title('MagPlot 6pF');   hold on; plot(f,line2,'--r');

figure(3);  %Plot the filter bode plot
plot(f,bodeAA);     title('Bode Plot 3pF'); hold on; plot(f,line3,'--r');

figure(4);  %Plot the filter bode plot
plot(f,bodeAB);     title('Bode Plot 6pF'); hold on; plot(f,line3,'--r');