i have these matlab codes :function $[$P$,$ A$]=$ horseshoeFlow$($d$,$ h$,$ n$,$ s$)$

$\%$ Constants

$\%$ 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

function data $=$ readData$($filename$)$

filename$=$'mix$1.$dat'

fid $=$ fopen$($filename$,\text{'}$r$\text{'})$;

if fid $==-1$

error$([\text{'}$Error opening file $\text{'}$ filename$])$;

end

data $=$ struct$()$;

while ~feof$($fid$)$

line $=$ fgetl$($fid$)$;

if isempty$($line$)$

continue;

end

$\%$ Split the line into parameter and value

$[$parameter$,$ value$]=$ strtok$($line$,\text{'}$:$\text{'})$;

$\%$ Remove the colon and whitespace from the value

value $=$ strtrim$($strrep$($value$,\text{'}$:$\text{'},\text{'}\text{'}))$;

$\%$ Remove comments $($lines starting with $\text{'}$#$\text{'})$

value $=$ strtok$($value$,\text{'}$#$\text{'})$;

$\%$ Convert the value to a number if it is not empty

if ~isempty$($value$)$

numValue $=$ str$2$double$($value$)$;

$\%$ Check if the conversion was successful

if ~isnan$($numValue$)$

data.$($strtrim$($parameter$))=$ numValue;

else

warning$([\text{'}$Could not convert value for parameter $\text{'}$ parameter $\text{'}$ to a number.$\text{'}])$;

end

else

warning$([\text{'}$Empty value for parameter $\text{'}$ parameter $\text{'}.\text{'}])$;

end

end

fclose$($fid$)$;

end

$\%$ Call readData function to get values from the $.$dat file

data $=$ readData$(\text{'}$mix$1.$dat'$)$;

$\%$ Extract values of d$,$ h$,$ and n from the structure

d $=$ data.DIAMETER; $\%$ Assuming 'DIAMETER' is the parameter for diameter

h $=$ data.DEPTH$1$; $\%$ Assuming 'DEPTH$1\text{'}$ is the parameter for flow depth

n $=$ data.MANNING;

s $=$ data.SLOPE; $\%$ Assuming 'MANNING' is the parameter for Manning's n

$\%$ Call the horseshoeFlow function with the obtained values

$[$P$,$ A$]=$ horseshoeFlow$($d$,$ h$,$ n$,$ s$)$;

$\%$ Display the results

fprintf$(\text{'}$Wetted Perimeter $($P$)$: $\%\mathrm{.}4$f

$\text{'},$ P$)$;

fprintf$(\text{'}$Cross$-$sectional Area $($A$)$: $\%\mathrm{.}4$f

$\text{'},$ A$)$;

which reads the mix$1.$dat files and gives me A nd P

How can I modify thse codes so that it can read multiple $.$dat files with almost the same structure and with the 'input' command I can write which $.$dat file I want it to take the data from and calculate A and P

for example I have 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 and mix$8.$dat:#Petaloid section

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

MANNING: $0\mathrm{.}016$

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

when I run the programm: enter file prefix:mix$8.$dat

and it gives me A and P for mix$8.$dat

thanks!MATLAB CODE ONLY!!