clear all;
clc;
close all;
syms x;
order = 15.74;

den_cmplx = i^order;

den_r = real(den_cmplx); %real part
den_c = imag(den_cmplx); %imaginary part


%To find the magnitude of this transferfunction, we need to substitute s =
%jw and multiply its complex conjugate with it.
t1 = 1/(1 + den_r*x^order + den_c*i*x^order); %just substitution of s = jw
t2 = 1/(1 + den_r*x^order - den_c*i*x^order); % 1/(complex conjugate of denominator of t1)
absolute = t1*t2; %Multiply t1 and t2 to get absolute function
pretty(absolute)
absolute

y=1./((x.^(787/50)*(4133199390247927/4503599627370496 - (7154380369095101*i)/18014398509481984) + 1).*(x.^(787/50)*(4133199390247927/4503599627370496 + (7154380369095101*i)/18014398509481984) + 1)); %Same function obtained from the above code evaluated for matrix x.
loglog(x,y,'-s')