function res = check_orthogonality(f,g)

eq = f*g;
res = int(eq,0,1);
end