% 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.11141e-07 6.80185e-07 4.16273e-06 2.54759e-05 0.000155913 0.000954186 0.00583962 0.0357384 0.218719 0.8 
0 1.23108e-06 6.76805e-06 3.64204e-05 0.000192242 0.000988924 0.00490414 0.0229871 0.0976804 0.33464 0.8 
0 6.84305e-06 3.39087e-05 0.000160994 0.00073674 0.00321336 0.0131187 0.0486727 0.156138 0.404764 0.8 
0 2.56633e-05 0.000115252 0.000485934 0.00194642 0.0073088 0.0252125 0.0776522 0.205811 0.45171 0.8 
0 7.36396e-05 0.000301906 0.00114046 0.00404892 0.0133364 0.0399998 0.106753 0.247183 0.48556 0.8 
0 0.000173707 0.000655598 0.00224164 0.00715105 0.0210538 0.0563154 0.134475 0.281777 0.51132 0.8 
0 0.00035303 0.0012374 0.00387054 0.0112436 0.0300923 0.0732777 0.160241 0.311025 0.531722 0.8 
0 0.000638874 0.00209773 0.00606442 0.0162375 0.0400758 0.090294 0.183929 0.336062 0.548377 0.8 
0 0.00105468 0.00327045 0.00882118 0.0220021 0.0506757 0.106994 0.205616 0.357755 0.562297 0.8 
0 0.00161736 0.00477134 0.0121091 0.0283939 0.0616274 0.123161 0.225464 0.376755 0.574151 0.8 
0 0.0023362 0.00659958 0.0158775 0.0352735 0.0727279 0.138679 0.243653 0.393558 0.584402 0.8 
0 0.00321298 0.00874089 0.0200651 0.0425151 0.0838268 0.153496 0.260361 0.408548 0.593379 0.8 
0 0.00424299 0.0111712 0.0246073 0.0500095 0.0948153 0.167599 0.275753 0.42202 0.601326 0.8 
0 0.00541648 0.0138601 0.02944 0.0576646 0.105617 0.181 0.289975 0.434211 0.608424 0.8 
0 0.00672015 0.0167736 0.0345026 0.0654044 0.116178 0.193724 0.303156 0.445308 0.614813 0.8 
0 0.00813852 0.0198767 0.0397396 0.0731667 0.126465 0.205802 0.315407 0.455464 0.620605 0.8 
0 0.00965504 0.0231349 0.0451011 0.0809013 0.136455 0.217269 0.326827 0.464802 0.625887 0.8 
0 0.0112531 0.0265151 0.0505431 0.0885681 0.146135 0.228163 0.3375 0.473425 0.630729 0.8 
0 0.0129164 0.0299869 0.0560272 0.0961353 0.155501 0.238517 0.3475 0.481419 0.635188 0.8 
0 0.0146298 0.0335224 0.0615205 0.103578 0.164551 0.248367 0.356892 0.488856 0.639312 0.8 
0 0.0163793 0.0370968 0.0669948 0.110877 0.173288 0.257744 0.36573 0.495796 0.643141 0.8 
0 0.0181522 0.040688 0.0724264 0.118017 0.181716 0.266677 0.374065 0.502292 0.646708 0.8 
0 0.0199371 0.0442771 0.0777955 0.124987 0.189843 0.275196 0.38194 0.508387 0.650042 0.8 
0 0.0217243 0.0478476 0.0830857 0.131779 0.197676 0.283325 0.389393 0.51412 0.653165 0.8 
0 0.023505 0.0513854 0.0882837 0.138389 0.205223 0.291088 0.396457 0.519524 0.656099 0.8 
0 0.0252721 0.0548789 0.0933787 0.144811 0.212493 0.298506 0.403163 0.524628 0.658862 0.8 
0 0.0270192 0.0583182 0.0983622 0.151046 0.219496 0.305599 0.409536 0.529457 0.661469 0.8 
0 0.0287411 0.0616952 0.103228 0.157092 0.226239 0.312387 0.415601 0.534034 0.663934 0.8 
0 0.0304335 0.0650037 0.10797 0.16295 0.232733 0.318884 0.421379 0.538379 0.666267 0.8 
0 0.032093 0.0682384 0.112587 0.168621 0.238985 0.325108 0.426889 0.542508 0.668481 0.8 
0 0.0337168 0.0713954 0.117074 0.174109 0.245004 0.331073 0.432149 0.546437 0.670583 0.8 
0 0.0353026 0.0744718 0.121433 0.179416 0.250799 0.336791 0.437173 0.550179 0.672582 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
