HERE IS CODE : clear all

close all

clc

g $=1\mathrm{.}4$;

load U$.$mat;

for i$=1$:size$($U$,1)$

$\%$W$($i$,$:$)=$U$2$W$($U$($i$,$:$),$g$)$;

$\%$lambda$($i$,$:$)=$eigsEuler$($U$($i$,$:$),$g$)$;

$\%$F$($i$,$:$)=$fluxEuler$($U$($i$,$:$),$g$)$;

end

fprintf$(\text{'}$

$---$ W $----\text{'})$;

for i$=1$:size$($U$,1)$

fprintf$(\text{'}$

$\%12\mathrm{.}3$e $\%12\mathrm{.}3$e $\%12\mathrm{.}3$e$\text{'},$W$($i$,1$:$3))$;

end

fprintf$(\text{'}$

$\text{'})$;

fprintf$(\text{'}$

$---$ lambda $----\text{'})$;

for i$=1$:size$($U$,1)$

fprintf$(\text{'}$

$\%12\mathrm{.}3$e $\%12\mathrm{.}3$e $\%12\mathrm{.}3$e$\text{'},$lambda$($i$,1$:$3))$;

end

fprintf$(\text{'}$

$\text{'})$;

fprintf$(\text{'}$

$---$ F $----\text{'})$;

for i$=1$:size$($U$,1)$

fprintf$(\text{'}$

$\%12\mathrm{.}3$e $\%12\mathrm{.}3$e $\%12\mathrm{.}3$e$\text{'},$F$($i$,1$:$3))$;

end

fprintf$(\text{'}$

$\text{'})$;