if youre familiar with MatLab, i jeed your help with my ode45 function code!

  

clc clear close all %Define the system parameters m-7.77; % slugs j-16; %lbs ft*s^2 kf-1100; % lbs/ft kr-1050; % lbs/ft bf=105; %lbs*s/ft br=200; %lbs*s/ft 1=2.5; % ft g=32.2; % ft/s^2 %Time span and initial conditions tspan=[0,2]; x0-[0,0,0,10,-16,0]'; I odefunc=@(t,x)StateEqFunc(t,x,kf, kr, bf, br,j); %Call ode45 [T,X]=ode45(odefunc, tspan,x0); %Extract solution and plot x1=X(:,1); x2=x(:,2); x3=X(:,3); x4=x(:,4); x5=x(:,5); x6-X(:,6); figure