i have this matlab code $.$It has a function which with input arguments the diameter d and the flow depth h of the horseshoe

cross section, returns the wetted perimeter P and the wetted cross$-$sectional area A$.$ :function $[$P$,$ A$]=$ horseshoeFlow$($d$,$ h$)$

$\%$ Calculate h$1$

h$1=($d$./2).*(1-(1+(7^(1/2))/4))$;

if h $>=0$ && h $<$ h$1$

$\%$ Case $1$: $0<=$ h $<$ h$1$

A $=($h $-$ d$).*($h $.*(2*$ d $-$ h$)).^(1/2)+$ d$.^2.*(($sin$(($h $-$ d$)./$d$)).^(-1)+$ pi$/2)$;

Wetp $=2*$ d $.*($cos$(1-$ h$./$d$)).^(-1)$;

elseif h$1<=$ h && h $<($d$/2)$

$\%$ Case $2$: h$1<=$ h $<$ diameter$/2$

h$2=($d $/2)-$ h$1$;

C$1=1-(1+7^(1/2))/2$;

C$2=$ C$1/2*(1-(1-($C$1^2)/4).^(1/2))-$ sin$($C$1/2)^(-1)$;

A$1=($h $-$ d$/2).*($h $.*($d $-$ h$)).^(1/2)+($d$^2/4).*$ sin$((2*$ h $-$ d$)/$ d$).^(-1)$;

A$2=$ d$^2.*($C$2+($sin$(2*$ h $-$ d$)/(2*$ d$)).^(-1))-($h $-$ d$/2).*($d $-($d$^2-($h $-$ d$/2)^2).^(1/2))+$ A$1$;

Wetp$1=2*$ d $.*(($cos$(($d $-2*$ h$)/(2*$ d$))).^(-1)-($cos$(-$C$1/2)).^(-1))$;

Wetp$2=2*$ d $.*(($cos$(($d $-2*$ h$)/(2*$ d$))).^(-1)-($cos$(-$C$1/2)).^(-1))+$ Wetp$1$;

A $=$ A$2$;

Wetp $=$ Wetp$2$;

elseif $($d$/2)<=$ h && h $<$ d

$\%$ Case $3$: diameter$/2<=$ h $<$ diameter

h$3=$ d $/2$;

C$1=1-(1+7^(1/2))/2$;

C$2=$ C$1/2*(1-(1-($C$1^2)/4).^(1/2))-$ sin$($C$1/2)^(-1)$;

A$2=($h $-$ d$/2).*($h $.*($d $-$ h$)).^(1/2)+($d$^2/4).*$ sin$((2*$ h $-$ d$)/$ d$).^(-1)$;

A$3=($h $-$ d$/2).*($h $.*($d $-$ h$)).^(1/2)+($d$^2/4).*$ sin$((2*$ h $-$ d$)/$ d$).^(-1)+$ A$2$;

Wetp$2=2*$ d $.*(($cos$(($d $-2*$ h$)/(2*$ d$))).^(-1)-($cos$(-$C$2/2)).^(-1))$;

Wetp$3=$ d $.*(($cos$(1-(2*$ h$)/$ d$)).^(-1)-$ pi$/2)+$ Wetp$2$;

A $=$ A$3$;

Wetp $=$ Wetp$3$;

else

error$(\text{'}$Invalid input for flow depth h$.\text{'})$;

end

$\%$ Calculate hydraulic radius, velocity, and flow rate

Rh $=$ A $/$ Wetp;

V $=(1/$ n$)*$ Rh$^(2/3)*($s$)^(1/2)$;

Q $=$ V $*$ A;

$\%$ Output wetted perimeter and cross$-$sectional area

P $=$ Wetp;

A $=$ A;

end

how can I make a code that reads the diameter,depth$($h$),$manning,slope and discharge from $.$dat files? for exmple this one named mix$1.$dat:MANNING: $0\mathrm{.}016$

DIAMETER: $1\mathrm{.}0$ #m

SLOPE: $0\mathrm{.}5$ # per cent $($actual slope is $0\mathrm{.}005)$

DEPTH$1$: $0\mathrm{.}4$ #m

DISCHARGE$1$: $0\mathrm{.}6$ #m$3/$sec

It has to read other dat files too they are smilar but the have different orders fot the data

thanks! $!!$IN MATLAB CODE PLEASE!!