% Matlab plot script for ppc_fd1() finite difference solver.
PLOT_WITH_MATLAB_FD_XVEC = [
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 
];
PLOT_WITH_MATLAB_FD_TVEC = [
0 0.015625 0.03125 0.046875 0.0625 0.078125 0.09375 0.109375 0.125 0.140625 0.15625 0.171875 0.1875 0.203125 0.21875 0.234375 0.25 0.265625 0.28125 0.296875 0.3125 0.328125 0.34375 0.359375 0.375 0.390625 0.40625 0.421875 0.4375 0.453125 0.46875 0.484375 0.5 
];
PLOT_WITH_MATLAB_FD_ZMAT = [
0 0 0 1 1 1 1 0 0 0 0 
0 0.0351854 0.25052 0.748515 0.958908 0.95891 0.74853 0.250623 0.0359082 0.00504329 0 
0 0.096542 0.327081 0.666953 0.876828 0.876846 0.667082 0.327831 0.100692 0.0213954 0 
0 0.138477 0.357666 0.624119 0.80126 0.801348 0.624665 0.360296 0.150038 0.0445904 0 
0 0.162859 0.369839 0.591879 0.738163 0.73844 0.593402 0.375983 0.184769 0.0665631 0 
0 0.175538 0.37187 0.563541 0.68588 0.686537 0.566762 0.382973 0.208769 0.0844403 0 
0 0.180781 0.367611 0.537294 0.641925 0.643185 0.542918 0.384582 0.224897 0.0979101 0 
0 0.181327 0.359399 0.512574 0.604238 0.606315 0.521148 0.382586 0.235151 0.107518 0 
0 0.178921 0.348782 0.489178 0.571299 0.574356 0.501029 0.378088 0.24099 0.113988 0 
0 0.174693 0.336798 0.467013 0.542015 0.54615 0.482252 0.371828 0.24351 0.117998 0 
0 0.16938 0.324145 0.446017 0.515606 0.520848 0.464575 0.364323 0.243541 0.120113 0 
0 0.163468 0.311284 0.42614 0.491506 0.497829 0.447812 0.355945 0.241716 0.120788 0 
0 0.157278 0.298514 0.407326 0.469302 0.476639 0.431819 0.346968 0.238516 0.120378 0 
0 0.151018 0.286026 0.389521 0.448687 0.45694 0.416491 0.337595 0.234307 0.119158 0 
0 0.144823 0.273936 0.372666 0.429426 0.438484 0.401749 0.327979 0.229369 0.117338 0 
0 0.138777 0.262311 0.356704 0.411341 0.421084 0.387534 0.318239 0.223914 0.115081 0 
0 0.132932 0.25118 0.341578 0.394291 0.4046 0.373805 0.308464 0.218105 0.11251 0 
0 0.127315 0.240553 0.327235 0.378164 0.388925 0.360532 0.298724 0.212065 0.109719 0 
0 0.121938 0.230426 0.313621 0.362869 0.373975 0.34769 0.289073 0.205891 0.106781 0 
0 0.116806 0.220785 0.30069 0.348332 0.359685 0.335264 0.279551 0.199654 0.10375 0 
0 0.111914 0.211608 0.288394 0.33449 0.346003 0.323238 0.270192 0.19341 0.100669 0 
0 0.107255 0.202876 0.276693 0.32129 0.332886 0.311601 0.261018 0.187204 0.0975704 0 
0 0.102819 0.194563 0.265547 0.308686 0.320298 0.300344 0.252049 0.181066 0.0944789 0 
0 0.0985962 0.186646 0.254922 0.29664 0.30821 0.289456 0.243298 0.175023 0.0914135 0 
0 0.0945743 0.179103 0.244784 0.285116 0.296595 0.27893 0.234774 0.169094 0.0883887 0 
0 0.0907423 0.171911 0.235104 0.274084 0.28543 0.268757 0.226484 0.163292 0.0854155 0 
0 0.0870894 0.165049 0.225854 0.263515 0.274695 0.258927 0.218432 0.157627 0.082502 0 
0 0.083605 0.158499 0.21701 0.253384 0.264369 0.249433 0.210618 0.152107 0.0796545 0 
0 0.0802792 0.152242 0.208548 0.243669 0.254436 0.240267 0.203044 0.146737 0.0768773 0 
0 0.0771029 0.14626 0.200448 0.234348 0.24488 0.231418 0.195708 0.14152 0.0741733 0 
0 0.0740675 0.140539 0.192689 0.225403 0.235684 0.22288 0.188607 0.136457 0.0715446 0 
0 0.0711649 0.135064 0.185254 0.216815 0.226836 0.214642 0.181739 0.131549 0.0689924 0 
0 0.0683878 0.129822 0.178127 0.208568 0.218322 0.206697 0.1751 0.126795 0.0665169 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
