% Matlab script for porous medium solver
PLOT_WITH_MATLAB_FD_XVEC = [
-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
];
PLOT_WITH_MATLAB_FD_TVEC = [
0 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 0.3 0.325 0.35 0.375 0.4 0.425 0.45 0.475 0.5 0.525 0.55 0.575 0.6 0.625 0.65 0.675 0.7 0.725 0.75 0.775 0.8 0.825 0.85 0.875 0.9 0.925 0.95 0.975 1 
];
PLOT_WITH_MATLAB_FD_ZMAT = [
0 0 0 0 0 -0.5 -0.5 -0.5 -0.5 -0.5 0.5 0.5 0.5 0.5 0.5 0.5 0 0 0 0 0 
0 -1.98603e-16 -1.32254e-06 -0.0101898 -0.201266 -0.371259 -0.4382 -0.451229 -0.429082 -0.325847 0.336699 0.434688 0.46443 0.469407 0.454413 0.400952 0.164934 0.00155004 1.16381e-09 -1.98603e-16 0 
0 -7.98681e-13 -8.61199e-05 -0.040571 -0.253221 -0.347003 -0.394846 -0.404927 -0.375709 -0.27439 0.293805 0.389608 0.427264 0.434658 0.416498 0.363246 0.241294 0.0243763 5.45853e-06 -1.98603e-16 0 
0 -8.12479e-10 -0.000865663 -0.0822711 -0.262993 -0.331107 -0.365189 -0.369434 -0.338985 -0.243653 0.27039 0.35894 0.396116 0.404657 0.388857 0.345591 0.264292 0.0654207 0.000234251 4.30392e-12 0 
0 -6.98788e-08 -0.00379381 -0.121322 -0.259995 -0.316427 -0.342472 -0.342633 -0.311995 -0.221049 0.254022 0.335808 0.371162 0.38042 0.367904 0.332607 0.267252 0.108596 0.00191673 3.11434e-09 0 
0 -1.51141e-06 -0.0103185 -0.150892 -0.254759 -0.302981 -0.323976 -0.321563 -0.29102 -0.20323 0.241355 0.317392 0.350982 0.360742 0.350909 0.321013 0.264856 0.144385 0.00710639 2.32721e-07 0 
0 -1.40324e-05 -0.0207376 -0.170034 -0.249977 -0.291209 -0.308507 -0.304369 -0.274053 -0.188655 0.231021 0.302257 0.334312 0.344398 0.336586 0.31067 0.26182 0.169509 0.0169763 4.6006e-06 0 
0 -7.50405e-05 -0.0341677 -0.180856 -0.245729 -0.281069 -0.295384 -0.289967 -0.259927 -0.176423 0.222318 0.289521 0.320251 0.330543 0.324313 0.301627 0.259019 0.184942 0.0310259 3.92329e-05 0 
0 -0.000274175 -0.0492445 -0.18612 -0.241637 -0.272267 -0.284112 -0.277677 -0.247914 -0.165954 0.214827 0.278602 0.308185 0.318613 0.313684 0.293713 0.256226 0.193395 0.0477618 0.000192645 0 
0 -0.000761521 -0.0646697 -0.188065 -0.237498 -0.264482 -0.274301 -0.267033 -0.237533 -0.156853 0.208277 0.269101 0.297683 0.30821 0.304373 0.28667 0.253206 0.197491 0.0655338 0.000649696 0 
0 -0.00172562 -0.0794137 -0.18821 -0.233304 -0.25747 -0.265643 -0.257694 -0.228443 -0.148837 0.202476 0.260732 0.288435 0.299033 0.296112 0.280278 0.249937 0.199137 0.0829204 0.00167957 0 
0 -0.00334906 -0.0927541 -0.187487 -0.229133 -0.25107 -0.257907 -0.249408 -0.220396 -0.141703 0.197286 0.253284 0.280207 0.290849 0.288693 0.274392 0.246522 0.199555 0.0988484 0.00357159 0 
0 -0.00576299 -0.104267 -0.186429 -0.225077 -0.245182 -0.250922 -0.24198 -0.213204 -0.135294 0.192602 0.246597 0.272818 0.283482 0.28196 0.26893 0.243094 0.199446 0.112637 0.00655108 0 
0 -0.00901981 -0.113795 -0.185298 -0.221203 -0.23974 -0.244563 -0.235264 -0.206721 -0.12949 0.188344 0.240546 0.26613 0.276793 0.275802 0.263846 0.239762 0.199158 0.123994 0.0107195 0 
0 -0.0130909 -0.121392 -0.184196 -0.21755 -0.234699 -0.238736 -0.229148 -0.200836 -0.124196 0.184447 0.235033 0.260033 0.270677 0.270136 0.259112 0.236589 0.198818 0.132953 0.0160422 0 
0 -0.0178835 -0.127248 -0.183133 -0.214129 -0.230022 -0.23337 -0.223542 -0.195458 -0.119337 0.180861 0.229979 0.254441 0.265053 0.264901 0.254704 0.233599 0.198433 0.139763 0.0223763 0 
0 -0.0232649 -0.13162 -0.182085 -0.210929 -0.225674 -0.228407 -0.218378 -0.190516 -0.114851 0.177545 0.225321 0.249284 0.259857 0.260047 0.250596 0.230788 0.197969 0.144776 0.029513 0 
0 -0.029085 -0.134783 -0.181018 -0.207932 -0.221625 -0.223801 -0.213599 -0.185954 -0.11069 0.174463 0.221007 0.244507 0.255036 0.255531 0.246759 0.228136 0.197392 0.148362 0.0372157 0 
0 -0.0351916 -0.136992 -0.17991 -0.205111 -0.217844 -0.219512 -0.209159 -0.181723 -0.106812 0.171588 0.216996 0.240063 0.250546 0.251318 0.243167 0.225618 0.196684 0.150862 0.0452455 0 
0 -0.0414402 -0.138469 -0.178748 -0.202443 -0.214303 -0.215507 -0.20502 -0.177786 -0.103184 0.168896 0.213253 0.235916 0.246352 0.247375 0.23979 0.22321 0.195845 0.152563 0.0533752 0 
0 -0.047699 -0.139398 -0.17753 -0.199906 -0.210974 -0.211755 -0.20115 -0.174109 -0.0997763 0.166367 0.209747 0.232033 0.242422 0.243673 0.236604 0.220894 0.194889 0.153696 0.0613989 0 
0 -0.0538515 -0.13993 -0.176265 -0.197483 -0.207836 -0.20823 -0.19752 -0.170666 -0.0965662 0.163985 0.206455 0.228386 0.238728 0.240187 0.233588 0.218657 0.193841 0.154444 0.0691374 0 
0 -0.0597978 -0.140181 -0.174962 -0.195159 -0.204867 -0.20491 -0.194107 -0.167432 -0.0935329 0.161735 0.203354 0.224952 0.235246 0.236894 0.230721 0.216491 0.192727 0.15494 0.076444 0 
0 -0.0654559 -0.140241 -0.173636 -0.192925 -0.202051 -0.201773 -0.190889 -0.164386 -0.0906588 0.159604 0.200426 0.22171 0.231957 0.233775 0.227989 0.214392 0.191574 0.155281 0.0832078 0 
0 -0.0707617 -0.140178 -0.172301 -0.190773 -0.199371 -0.198803 -0.187848 -0.161513 -0.0879288 0.157582 0.197655 0.218642 0.228841 0.230813 0.22538 0.212357 0.190403 0.155532 0.0893549 0 
0 -0.0756699 -0.14004 -0.170969 -0.188698 -0.196817 -0.195984 -0.184968 -0.158794 -0.0853297 0.155659 0.195028 0.215733 0.225882 0.227994 0.222883 0.210386 0.189234 0.155731 0.0948475 0 
0 -0.0801524 -0.139861 -0.169651 -0.186695 -0.194376 -0.193303 -0.182235 -0.156218 -0.0828497 0.153827 0.19253 0.212967 0.223067 0.225305 0.22049 0.20848 0.18808 0.155899 0.0996795 0 
0 -0.0841973 -0.13966 -0.168355 -0.184762 -0.192039 -0.190747 -0.179635 -0.153772 -0.0804787 0.152079 0.190152 0.210334 0.220384 0.222736 0.218194 0.206638 0.18695 0.156042 0.103871 0 
0 -0.0878064 -0.139451 -0.167088 -0.182897 -0.1898 -0.188307 -0.177159 -0.151445 -0.0782077 0.150407 0.187884 0.207822 0.217822 0.220278 0.21599 0.20486 0.18585 0.156159 0.10746 0 
0 -0.0909929 -0.139238 -0.165855 -0.181097 -0.187651 -0.185973 -0.174795 -0.149228 -0.0760285 0.148806 0.185717 0.205421 0.215371 0.217922 0.213872 0.203145 0.18478 0.156244 0.110497 0 
0 -0.0937785 -0.139022 -0.164657 -0.179359 -0.185587 -0.183738 -0.172536 -0.147112 -0.073934 0.147271 0.183643 0.203123 0.213023 0.215663 0.211835 0.20149 0.183741 0.156289 0.113038 0 
0 -0.0961907 -0.138801 -0.163496 -0.177683 -0.183602 -0.181594 -0.170373 -0.14509 -0.0719177 0.145796 0.181655 0.200921 0.210772 0.213493 0.209875 0.199893 0.182731 0.156289 0.115141 0 
0 -0.0982602 -0.138572 -0.162371 -0.176065 -0.181692 -0.179536 -0.1683 -0.143155 -0.0699739 0.144378 0.179747 0.198808 0.208609 0.211406 0.207988 0.198351 0.181745 0.156238 0.116862 0 
0 -0.100019 -0.13833 -0.16128 -0.174502 -0.179853 -0.177557 -0.166311 -0.1413 -0.0680973 0.143012 0.177914 0.196777 0.20653 0.209398 0.20617 0.196861 0.180781 0.156133 0.118253 0 
0 -0.1015 -0.138071 -0.160222 -0.172994 -0.178079 -0.175653 -0.164399 -0.139521 -0.0662832 0.141695 0.176151 0.194824 0.204529 0.207464 0.204416 0.19542 0.179835 0.155973 0.119359 0 
0 -0.102733 -0.137794 -0.159194 -0.171536 -0.176369 -0.17382 -0.16256 -0.137812 -0.0645273 0.140425 0.174453 0.192943 0.202601 0.205599 0.202722 0.194023 0.178904 0.15576 0.120222 0 
0 -0.103748 -0.137495 -0.158194 -0.170126 -0.174718 -0.172052 -0.16079 -0.136168 -0.0628258 0.139197 0.172815 0.191129 0.200742 0.2038 0.201085 0.192668 0.177984 0.155494 0.120878 0 
0 -0.104571 -0.137173 -0.15722 -0.168761 -0.173124 -0.170347 -0.159084 -0.134587 -0.0611751 0.13801 0.171235 0.189379 0.198948 0.202062 0.199502 0.191351 0.177075 0.155178 0.121358 0 
0 -0.105226 -0.136827 -0.156268 -0.167439 -0.171582 -0.1687 -0.157438 -0.133063 -0.0595721 0.136861 0.169708 0.187689 0.197215 0.200381 0.197968 0.190069 0.176174 0.154817 0.121689 0 
0 -0.105734 -0.136458 -0.155338 -0.166158 -0.170091 -0.167109 -0.155849 -0.131595 -0.0580137 0.135748 0.168232 0.186056 0.195539 0.198755 0.196482 0.18882 0.175279 0.154415 0.121892 0 
0 -0.106115 -0.136066 -0.154426 -0.164914 -0.168647 -0.16557 -0.154314 -0.130177 -0.0564974 0.134668 0.166804 0.184475 0.193916 0.19718 0.195039 0.187601 0.17439 0.153974 0.121987 0 
];
surf(PLOT_WITH_MATLAB_FD_XVEC, PLOT_WITH_MATLAB_FD_TVEC, PLOT_WITH_MATLAB_FD_ZMAT)
xlabel('x')
ylabel('t')
zlabel('u')
colormap jet
