% 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 0.390625 0.609375 1 1 0.609375 0.390625 0 0 0 
0 0.152588 0.323486 0.676514 0.847412 0.847412 0.676514 0.323486 0.152588 0 0 
0 0.15974 0.39463 0.60537 0.780655 0.780655 0.60537 0.39463 0.15974 0.0596046 0 
0 0.189096 0.385197 0.59152 0.712184 0.712184 0.59152 0.385197 0.212379 0.0754371 0 
0 0.191832 0.38919 0.558059 0.66505 0.66505 0.558059 0.398285 0.226393 0.0994623 0 
0 0.193991 0.378062 0.533888 0.623257 0.623257 0.537441 0.393551 0.243956 0.110192 0 
0 0.190116 0.367029 0.507928 0.588347 0.589735 0.514756 0.391322 0.25014 0.1194 0 
0 0.184958 0.352961 0.484303 0.557475 0.559904 0.495828 0.384389 0.254219 0.12383 0 
0 0.178335 0.33864 0.46158 0.529841 0.533926 0.477327 0.377072 0.254133 0.126392 0 
0 0.171292 0.324045 0.440221 0.504772 0.510221 0.460274 0.368211 0.252258 0.126919 0 
0 0.16405 0.309757 0.420055 0.481686 0.488582 0.443823 0.358879 0.248592 0.126302 0 
0 0.156885 0.295925 0.401044 0.460305 0.468404 0.428126 0.348979 0.243903 0.124735 0 
0 0.149914 0.282675 0.383131 0.44032 0.449507 0.412943 0.33885 0.238398 0.12256 0 
0 0.143214 0.270056 0.36623 0.421569 0.431635 0.398283 0.328554 0.232388 0.119934 0 
0 0.136819 0.258076 0.350279 0.403884 0.414675 0.384073 0.318227 0.226026 0.117012 0 
0 0.13074 0.246727 0.335202 0.38716 0.398506 0.370306 0.307932 0.219458 0.113888 0 
0 0.124977 0.23598 0.320937 0.371296 0.383058 0.356957 0.297737 0.21278 0.110639 0 
0 0.119518 0.225806 0.307422 0.356219 0.368268 0.34402 0.287683 0.206067 0.107319 0 
0 0.11435 0.216169 0.294602 0.341864 0.354089 0.331485 0.277809 0.199375 0.103971 0 
0 0.109455 0.207034 0.282426 0.328178 0.340484 0.319348 0.268138 0.192746 0.100625 0 
0 0.104816 0.198367 0.270848 0.315113 0.327421 0.3076 0.258692 0.186211 0.097303 0 
0 0.100416 0.190136 0.259826 0.30263 0.314871 0.296238 0.249484 0.179794 0.0940237 0 
0 0.0962379 0.182312 0.249324 0.290691 0.30281 0.285253 0.240525 0.173513 0.0907998 0 
0 0.0922676 0.174866 0.239306 0.279266 0.291218 0.274639 0.23182 0.167379 0.0876408 0 
0 0.0884904 0.167773 0.229743 0.268326 0.280073 0.264389 0.223374 0.161404 0.084554 0 
0 0.0848936 0.16101 0.220607 0.257843 0.269358 0.254494 0.215188 0.155592 0.0815445 0 
0 0.0814652 0.154557 0.211873 0.247796 0.259054 0.244946 0.207262 0.149947 0.0786158 0 
0 0.0781945 0.148394 0.203516 0.238161 0.249145 0.235737 0.199594 0.144472 0.0757702 0 
0 0.0750716 0.142505 0.195517 0.228919 0.239617 0.226856 0.19218 0.139167 0.0730091 0 
0 0.0720878 0.136872 0.187857 0.22005 0.230453 0.218296 0.185017 0.134032 0.070333 0 
0 0.0692347 0.131481 0.180516 0.211538 0.22164 0.210045 0.178101 0.129066 0.0677418 0 
0 0.0665051 0.126321 0.17348 0.203367 0.213165 0.202096 0.171425 0.124265 0.0652348 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
