% 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 2.45138e-08 0.661438 0.866025 0.968246 1 0.968246 0.866025 0.661438 0 0 0 0 0 0 0 
0 -1.98603e-16 -1.98603e-16 8.83548e-10 0.000890786 0.0893212 0.416343 0.566436 0.731071 0.824234 0.855478 0.828626 0.741139 0.584481 0.335827 0.0175708 1.69717e-06 -1.98603e-16 -1.98603e-16 -1.98603e-16 0 
0 -1.98603e-16 5.66018e-15 1.67256e-05 0.0237397 0.260019 0.452908 0.566408 0.6625 0.72887 0.752545 0.730834 0.66438 0.558526 0.420925 0.167567 0.00217955 3.2432e-09 -1.98603e-16 -1.98603e-16 0 
0 -1.98603e-16 2.01698e-09 0.00117291 0.0969608 0.347909 0.463413 0.550168 0.615985 0.661116 0.677299 0.660703 0.612719 0.538736 0.439787 0.288688 0.0367423 1.85174e-05 1.68812e-15 -1.98603e-16 0 
0 -1.98603e-16 1.91344e-06 0.0115167 0.194562 0.369146 0.463917 0.531638 0.580808 0.612476 0.623238 0.610644 0.575531 0.520419 0.444179 0.336392 0.115783 0.00116883 5.23105e-10 -1.98603e-16 0 
0 1.39866e-12 0.000103822 0.0430886 0.262034 0.378347 0.457196 0.513192 0.552367 0.576192 0.583704 0.573815 0.546747 0.50306 0.440993 0.352821 0.196382 0.0113766 6.17825e-07 -1.98603e-16 0 
0 2.59734e-09 0.00127561 0.0950021 0.292515 0.384066 0.449314 0.496473 0.528964 0.548056 0.553716 0.545565 0.523497 0.487205 0.434875 0.361042 0.249081 0.0407268 4.5031e-05 2.94925e-14 0 
0 3.81409e-07 0.00670156 0.151593 0.304728 0.384805 0.441086 0.481705 0.509463 0.525481 0.53004 0.523045 0.50428 0.473163 0.428104 0.365318 0.275352 0.0858637 0.000689658 1.23955e-10 0 
0 1.17533e-05 0.0206429 0.197307 0.311778 0.382663 0.432681 0.468567 0.492915 0.50682 0.510678 0.504509 0.488076 0.460755 0.421145 0.366167 0.287771 0.134779 0.00415089 3.55691e-08 0 
0 0.000133855 0.0446282 0.226758 0.31746 0.379723 0.424608 0.45685 0.478642 0.491013 0.494392 0.488846 0.474146 0.449685 0.414181 0.364909 0.295049 0.176472 0.0139183 1.88017e-06 0 
0 0.000794498 0.0760377 0.243309 0.321373 0.376629 0.417149 0.446408 0.466178 0.47737 0.480397 0.475336 0.461976 0.439745 0.407476 0.362818 0.300538 0.20597 0.031881 3.19744e-05 0 
0 0.00300625 0.109967 0.252468 0.323227 0.37332 0.410278 0.437063 0.455178 0.465421 0.468169 0.4635 0.451217 0.430796 0.401202 0.36047 0.30465 0.224406 0.05683 0.000251535 0 
0 0.00820507 0.141448 0.258361 0.323555 0.36976 0.403872 0.428621 0.445366 0.454826 0.457348 0.453004 0.441616 0.422709 0.395375 0.35795 0.307214 0.235415 0.0856084 0.00116522 0 
0 0.0175666 0.16713 0.263131 0.323088 0.366075 0.397858 0.420923 0.436525 0.445333 0.447668 0.443602 0.432972 0.415351 0.389927 0.355232 0.308408 0.242349 0.11456 0.00371279 0 
0 0.0313739 0.186043 0.267355 0.322322 0.362429 0.392216 0.413853 0.42849 0.436748 0.43893 0.435104 0.425125 0.408601 0.38479 0.352355 0.308676 0.24742 0.140539 0.00899644 0 
0 0.0489448 0.199054 0.270879 0.321436 0.358928 0.386941 0.407335 0.421139 0.428925 0.430978 0.427362 0.417945 0.402364 0.379932 0.349416 0.308449 0.251701 0.161679 0.0177642 0 
0 0.0689596 0.207765 0.273507 0.320427 0.355593 0.382019 0.401307 0.414377 0.421752 0.423692 0.420261 0.411334 0.396573 0.375343 0.346513 0.307998 0.255467 0.177615 0.0300788 0 
0 0.0898214 0.213741 0.275244 0.319248 0.352402 0.377414 0.395714 0.40813 0.415138 0.416979 0.41371 0.405216 0.391181 0.371017 0.343705 0.307424 0.258637 0.189042 0.045369 0 
0 0.109954 0.218194 0.276265 0.317889 0.349319 0.373085 0.390503 0.402332 0.409011 0.410762 0.407638 0.39953 0.386146 0.366944 0.341006 0.306731 0.261107 0.197081 0.062669 0 
0 0.128081 0.221912 0.276801 0.316386 0.346327 0.368995 0.385626 0.396928 0.403311 0.404981 0.401988 0.394229 0.381432 0.363102 0.338402 0.305895 0.262881 0.202838 0.0808353 0 
0 0.143434 0.225285 0.277056 0.314797 0.343419 0.365114 0.381042 0.391871 0.397987 0.399584 0.396709 0.389268 0.377006 0.359465 0.335873 0.30491 0.264066 0.207207 0.0987152 0 
0.0507852 0.155842 0.228413 0.277165 0.313181 0.340605 0.361424 0.376719 0.387121 0.392996 0.394528 0.391762 0.38461 0.372834 0.356007 0.333406 0.303795 0.264818 0.210821 0.115465 0.0507852 
0.0769231 0.165999 0.231281 0.277193 0.311578 0.337893 0.357911 0.372631 0.382646 0.388301 0.389775 0.38711 0.380222 0.36889 0.352708 0.330996 0.30259 0.265293 0.214081 0.130712 0.0769231 
0.0943341 0.174354 0.23387 0.277168 0.310014 0.33529 0.354565 0.368757 0.378417 0.383874 0.385295 0.382722 0.376077 0.365148 0.349553 0.328648 0.30134 0.265617 0.2172 0.144252 0.0943341 
0.107575 0.18125 0.236164 0.277095 0.308499 0.332798 0.35138 0.36508 0.374413 0.379686 0.38106 0.378572 0.37215 0.361592 0.346534 0.326371 0.300087 0.265881 0.220244 0.155949 0.107575 
0.118227 0.18697 0.238163 0.276973 0.307034 0.330414 0.348345 0.361586 0.370614 0.375717 0.377046 0.374637 0.36842 0.358205 0.343646 0.324176 0.298865 0.266133 0.223193 0.165843 0.118227 
0.127075 0.191749 0.239876 0.276795 0.305614 0.328134 0.345452 0.358262 0.367004 0.371948 0.373234 0.370899 0.364873 0.354976 0.340884 0.32207 0.297692 0.266389 0.225985 0.174107 0.127075 
0.134582 0.195773 0.241324 0.27656 0.304236 0.32595 0.342691 0.355095 0.363568 0.368361 0.369608 0.36734 0.361493 0.351896 0.338246 0.320058 0.296576 0.266641 0.228563 0.180972 0.134582 
0.141049 0.199188 0.242533 0.276267 0.302892 0.323854 0.340053 0.352073 0.360292 0.364943 0.366152 0.363948 0.35827 0.348957 0.335727 0.318137 0.295517 0.266876 0.230887 0.186679 0.141049 
0.146683 0.202107 0.243531 0.275917 0.301578 0.321839 0.337528 0.349186 0.357164 0.361682 0.362854 0.36071 0.355192 0.34615 0.333322 0.316306 0.294506 0.267076 0.232943 0.191445 0.146683 
0.151636 0.204619 0.244344 0.275514 0.300289 0.319897 0.335108 0.346424 0.354175 0.358565 0.359703 0.357617 0.352251 0.343468 0.331025 0.314558 0.293535 0.267226 0.234733 0.195451 0.151636 
0.156022 0.206793 0.244997 0.275064 0.299022 0.318022 0.332784 0.343778 0.351313 0.355582 0.356689 0.354658 0.349439 0.340903 0.328829 0.312885 0.292595 0.267318 0.236273 0.198842 0.156022 
0.159931 0.208683 0.245511 0.274571 0.297776 0.316209 0.330549 0.341239 0.34857 0.352725 0.353802 0.351824 0.346746 0.338448 0.326727 0.311279 0.291677 0.267347 0.237586 0.201734 0.159931 
0.163431 0.210332 0.245906 0.274042 0.296548 0.314452 0.328396 0.338799 0.345938 0.349985 0.351034 0.349108 0.344166 0.336096 0.324711 0.309733 0.290774 0.267313 0.238695 0.204215 0.163431 
0.16658 0.211777 0.246198 0.27348 0.295339 0.312748 0.32632 0.336453 0.343409 0.347354 0.348378 0.346502 0.34169 0.333838 0.322774 0.308241 0.289882 0.267216 0.239625 0.206358 0.16658 
0.169423 0.213045 0.246401 0.272892 0.294147 0.311093 0.324315 0.334192 0.340977 0.344826 0.345826 0.344 0.339313 0.33167 0.320911 0.306797 0.288997 0.267062 0.240396 0.208216 0.169423 
0.171999 0.21416 0.246527 0.272282 0.292972 0.309484 0.322376 0.332013 0.338635 0.342393 0.343372 0.341594 0.337027 0.329584 0.319114 0.305395 0.288117 0.266855 0.24103 0.209835 0.171999 
0.174339 0.215142 0.246587 0.271654 0.291815 0.307919 0.3205 0.329909 0.336377 0.340051 0.34101 0.339278 0.334827 0.327574 0.317378 0.304032 0.287241 0.2666 0.241544 0.21125 0.174339 
0.176471 0.216006 0.246588 0.27101 0.290674 0.306394 0.318682 0.327876 0.334199 0.337792 0.338734 0.337047 0.332707 0.325635 0.3157 0.302703 0.286367 0.266301 0.241953 0.21249 0.176471 
0.178417 0.216767 0.246539 0.270355 0.289551 0.304908 0.316919 0.32591 0.332095 0.335613 0.336539 0.334896 0.330661 0.323762 0.314073 0.301407 0.285496 0.265965 0.24227 0.213579 0.178417 
0.180197 0.217437 0.246447 0.269691 0.288446 0.303459 0.315208 0.324006 0.330062 0.333508 0.334419 0.332819 0.328685 0.321951 0.312496 0.300141 0.284629 0.265595 0.242508 0.214537 0.180197 
];
surf(PLOT_WITH_MATLAB_FD_XVEC, PLOT_WITH_MATLAB_FD_TVEC, PLOT_WITH_MATLAB_FD_ZMAT)
xlabel('x')
ylabel('t')
zlabel('u')
colormap jet
