Pb integration ds matlab

Bonjour, dans le cadre d'un projet de mathématiques, je suis amené à tracer pour différentes valeurs du paramètre oméga l'intégrale d'une fonction f dépendant d'oméga sur des bornes y1 et y2 dépendant de oméga.

Dans le code suivant, cependant, je n'arrive pas à paramétrer f. Faut il définir plusieurs paramètres à f et intégrer avec une autre fonction que quad? Comment rendre mon "for omega=0.7:0.01:0.8" effectif pour "f" et "trace"?

merci de m'éclairer un peu :)

function Trace

for omega=0.7:0.01:0.8

C=0.5-sqrt(0.5)-omega.*omega.*0.5;

%y2 est la valeur maximale
y2=sqrt(power((-1+(2*sqrt(2)*(1-omega*omega-sqrt(2)/4))/(2*omega*omega-2)),2)-1);

%y4 est la valeur minimale
y4=-sqrt(power((-1+(2*sqrt(2)*(1-omega*omega-sqrt(2)/4))/(2*omega*omega-2)),2)-1);


%Calcul de la periode
T =quad(@f,y4,y2)

% POINTS w,T
plot(omega,abs(T),'r:')
%end

%fonction a integrer
function f=f(x); 
C=0.5-sqrt(0.5)-omega*omega*0.5
f=sqrt(2)./(sqrt(-0.5.*x.*x+0.5.*sqrt(x.*x+1)+0.5.*x.*x.*omega.*omega+C));