Pretty please help me fix my MATLAB code, I've been struggling for quite some time. I hope I've included enough context information to help with the MATLAB portion of the project.

Notes

utilize MATLAB function "ode$45"$ to solve the system equations

Each plot should include $2$ curves

The first curve should show the pedal force $($Fp $).$ The second curve should show the CVT shift angle $(\backslash$theta $)$ or theta.

Both curves should be plotted as a function of time. The simulation should start at time t $=0$ and continue at least until the pedal force falls within $2\%$ of the desired FD $($the response time$)\%$ script that represents an accurate model of the automatic bicycle transmission

$\%$ Initial shift angle $($deg$)$: $-15\mathrm{.}0$

$\%$

$\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%$

$\%$

$\%$ author $[335]$

$\%$ version $[1\mathrm{.}0]$

$\%$ date $[2024\mathrm{.}09\mathrm{.}27]$

$\%$

$\%$

$\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%$

$\%$ DESCRIPTION OF VARIABLES

$\%$

$\%$

$\%$ START OF EXECUTABLE CODE

$\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%$

$\%$ List of parameters in structure form

close; clear; clc;

params.wi $=(2*$pi$)/(0\mathrm{.}121)$; $\%$rad$/$s

params.c$3=384$;

params.ri $=20\mathrm{.}7$; $\%$mm

params.c$2=3\mathrm{.}176$;

params.c$1=0\mathrm{.}121$;

params.rp $=170$; $\%$mm

params.rv $=25\mathrm{.}9$; $\%$mm

params.k $=2\mathrm{.}24$;$\%$N$/$mm

params.T$\_$level $=(8\mathrm{.}6*1000)$; $\%$N$-$mm

params.theta$0=-15*($pi$/180)$; $\%$rad

tspan $=$ linspace$(0,15,1000)$; $\%$s

Tin $=$ zeros$($length$($tspan$),1)$;

Fp $=$ zeros$($length$($tspan$),1)$;

$\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%$

$[$t$,$ theta$]=$ ode$45($@$($t$,$ theta$)$ eq$1($t$,$ theta, params$),$ tspan, params.theta$0)$;

plot$($t$,$ theta$*180/$pi$,\text{'}$r$\text{'})$;

hold on;

$\%[$t$,$ Fp$]=$ ode$45($@$($t$,$ Fp$)$ eq$5($t$,$ theta, params$),$ tspan, params.theta$0)$;

$\%$hold off;

xlabel$(\text{'}$Time $($s$)\text{'})$;

ylabel$(\text{'}$Theta and Fp$\text{'})$;

for i$=1$:length$($t$)$

Tin$($i$)=$ eq$6($params$,$ t$($i$),$ theta$($i$))$;

Fp$($i$)=$ Tin$($i$)/($

params.cl $*$params$.$rp$)$;

end

plot$($t$,$ Fp$,\text{'}$b$\text{'})$

$\%$ Equation $1$: de$/$dt

function dtheta$\_$dt $=$ eq$1($t$,$ theta, params$)$

x $=$ eq$2($params$,$ t$,$ theta$)$;

dtheta$\_$dt $=($params$.$wi $/($params$.$c$3*$ params.ri$))*$ x;

end

$\%$ Equation $2$: x$($t$)$

function x $=$ eq$2($params$,\}$t$\backslash$mathrm$\{,$ theta$)$

F$\_$set $=$ eq$4($params$)$;

$\%$Fp$=$ eq$5($params$)$;

Fc$\}=\backslash$mathrm$\{$ eq$3($params$,$ t$,$ theta$)$;

x $=(1/$params$.$k$)*(2*$Fc $-$ F$\_$set$)$;

end

$\%$ Equation $3$: Fc$($t$)$

function FC $=$ eq$3($params$,$ t$,$ theta$)$

Fp $=$ eq$5($params$,$ t$,$ theta$)$;

FC $=(($params$.$c$1*$ params.rp$)/$ params.rv$)*$ Fp;

end

$\%$ Equation $4$: F$\_$set$($t$)$

function F$\_$set $=$ eq$4($params$)$

F$\_$set $=(($params$.$c$1*$ params.rp$)/$ params.rv$)*156$;

end

$\%$ Equation $5$: Fp$($t$)$

function Fp $=$ eq$5($params$,$ t$,$ theta$)$

Tin eq$6($params$,$ t$,$ theta$)$;

Fp $=$ Tin $/($params$.$c$1*$ params.rp$)$;

end

$\%$ Equation $6$: Tin$($t$)$

function Tin $=$ eq$6($params$,$ t$,$ theta$)$

TRN $=$ eq$7($params$,$ t$)$;

Tin $=($TRW $/$ params.c$2)*(1-$tan$($theta$))/(1+$tan$($theta$))$;

end

$\%$ Equation $7$: TRW

function TRW $=$ eq$7($params$,$ t$)$

end