error_front = 0;
error_rear = 0;
rear_x = zeros(1,length(T));
front_x = zeros(1,length(T));
delta_fun=@(t) pi/180*sin(2*pi*t)-0.00175;

for i = 1:length(T)

    v = Z_eq(4,i);
    r = Z_eq(5,i);

    % Linear Model

    F_zf_linear = b/L*m*g;
    F_zr_linear = a/L*m*g;

    Ca_r_linear = F_zr_linear * B * C * D;
    Ca_f_linear = F_zf_linear * B * C * D;

    a_f_linear = delta_fun(T(i)) - (v + a * r)/u_0;
    a_r_linear = -(v - b * r)/u_0;

    F_yf_linear = Ca_f_linear * a_f_linear;
    F_yr_linear = Ca_r_linear * a_r_linear;


    % Magic Model
    F_zf_magic = b/L * m*g;
    F_zr_magic = a/L * m*g;

    a_f_magic = delta_fun(T(i)) - atan((v + a * r)/u_0);
    a_r_magic = -atan((v - b * r)/u_0);

    F_yf_magic = F_zf * D * sin(C * atan( B * (1-E) * a_f_magic + E * atan(B * a_f_magic)));
    F_yr_magic = F_zr * D * sin(C * atan( B * (1-E) * a_r_magic + E * atan(B * a_r_magic)));


    % Error calculations
    temp = abs((F_yf_linear - F_yf_magic)/F_yf_magic);
    if temp > error_front
        error_front = temp;
    end

    temp = abs((F_yr_linear - F_yr_magic)/F_yr_magic);
    if temp > error_rear
        error_rear = temp;
    end
end

tireforce_percent_error = 100 * max(error_front,error_rear)