% 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 0 0 0 0 0 0 0 0.8 
0 1.41836e-06 6.46772e-06 2.80744e-05 0.000121552 0.000526201 0.00227793 0.00986114 0.0426889 0.1848 0.8 
0 9.29285e-06 3.87444e-05 0.000150824 0.000577143 0.00216978 0.00796996 0.0283418 0.096024 0.300244 0.8 
0 3.38227e-05 0.000130442 0.000461806 0.00158929 0.00530785 0.0170599 0.052082 0.147879 0.376426 0.8 
0 9.04175e-05 0.000325718 0.00106093 0.00332988 0.0100547 0.0289317 0.0782004 0.194332 0.429382 0.8 
0 0.000198313 0.00067284 0.002036 0.00589535 0.0163223 0.0427943 0.104755 0.234695 0.467963 0.8 
0 0.000378349 0.00121759 0.00345139 0.00930857 0.0239036 0.0579066 0.130597 0.269443 0.497243 0.8 
0 0.000650591 0.00199812 0.00534382 0.0135341 0.0325419 0.0736637 0.155124 0.299373 0.520245 0.8 
0 0.00103233 0.00304191 0.00772357 0.0184974 0.0419773 0.0896119 0.178074 0.325287 0.538841 0.8 
0 0.00153671 0.00436462 0.0105787 0.0241018 0.0519723 0.10543 0.199382 0.347885 0.554236 0.8 
0 0.00217205 0.00597057 0.0138803 0.0302422 0.0623233 0.120903 0.219093 0.367744 0.567234 0.8 
0 0.00294177 0.00785405 0.017588 0.0368139 0.0728632 0.135894 0.237304 0.385333 0.578389 0.8 
0 0.00384479 0.0100013 0.0216547 0.0437187 0.0834591 0.150325 0.254134 0.401026 0.588092 0.8 
0 0.00487611 0.0123924 0.0260299 0.050868 0.0940084 0.164155 0.269706 0.415124 0.596632 0.8 
0 0.00602767 0.0150034 0.0306631 0.0581839 0.104433 0.177371 0.284142 0.427868 0.60422 0.8 
0 0.00728916 0.0178077 0.0355054 0.0655994 0.114677 0.189978 0.297552 0.439455 0.611022 0.8 
0 0.00864872 0.0207779 0.0405109 0.0730577 0.124698 0.201991 0.310037 0.450045 0.617162 0.8 
0 0.0100937 0.0238865 0.0456372 0.0805114 0.134467 0.213432 0.321688 0.459769 0.622742 0.8 
0 0.0116111 0.0271068 0.0508464 0.0879214 0.143966 0.224328 0.332584 0.468735 0.627841 0.8 
0 0.0131882 0.0304136 0.0561044 0.0952558 0.153183 0.234706 0.342798 0.477036 0.632524 0.8 
0 0.0148125 0.0337835 0.0613815 0.102489 0.162112 0.244595 0.352391 0.484747 0.636844 0.8 
0 0.0164726 0.0371952 0.0666516 0.1096 0.170752 0.254021 0.36142 0.491934 0.640845 0.8 
0 0.0181578 0.0406294 0.0718925 0.116572 0.179103 0.263011 0.369934 0.498651 0.644565 0.8 
0 0.0198581 0.044069 0.0770852 0.123394 0.187168 0.27159 0.377975 0.504946 0.648033 0.8 
0 0.0215648 0.0474989 0.0822134 0.130056 0.194954 0.279782 0.385583 0.51086 0.651278 0.8 
0 0.0232701 0.050906 0.087264 0.136551 0.202466 0.28761 0.392792 0.516429 0.65432 0.8 
0 0.0249672 0.054279 0.0922257 0.142874 0.20971 0.295093 0.399632 0.521682 0.65718 0.8 
0 0.0266501 0.0576083 0.0970897 0.149023 0.216696 0.302251 0.406131 0.526647 0.659875 0.8 
0 0.0283136 0.0608859 0.101849 0.154995 0.223429 0.309102 0.412312 0.531348 0.662418 0.8 
0 0.0299535 0.064105 0.106497 0.16079 0.229919 0.315662 0.418199 0.535806 0.664824 0.8 
0 0.031566 0.0672601 0.111031 0.166409 0.236173 0.321947 0.42381 0.540038 0.667102 0.8 
0 0.0331482 0.0703469 0.115448 0.171854 0.242199 0.327971 0.429164 0.544063 0.669264 0.8 
0 0.0346976 0.0733618 0.119744 0.177126 0.248005 0.333747 0.434277 0.547894 0.671318 0.8 
];
surf ( PLOT_WITH_MATLAB_FD_XVEC, PLOT_WITH_MATLAB_FD_TVEC, PLOT_WITH_MATLAB_FD_ZMAT )
xlabel('x')
ylabel('t')
zlabel('u')
colormap jet
