% 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.0409178 0.186585 0.809911 0.94661 0.946631 0.810026 0.187089 0.0431012 0.00945201 0 
0 0.0840699 0.278609 0.708728 0.879818 0.879922 0.709249 0.280588 0.0912833 0.0253247 0 
0 0.118337 0.324399 0.647684 0.814694 0.814988 0.649051 0.329005 0.132906 0.0433635 0 
0 0.142328 0.346074 0.605306 0.756053 0.756678 0.608029 0.354367 0.16563 0.0606669 0 
0 0.157601 0.354299 0.572053 0.70468 0.705791 0.57663 0.367085 0.190098 0.0757467 0 
0 0.166236 0.354578 0.543633 0.659931 0.661675 0.550481 0.372348 0.207691 0.0880706 0 
0 0.170074 0.349973 0.518083 0.620788 0.623286 0.527508 0.372919 0.219792 0.0976431 0 
0 0.17055 0.342317 0.494487 0.586248 0.589589 0.506664 0.370377 0.227581 0.104725 0 
0 0.168724 0.332774 0.472394 0.555456 0.559688 0.487377 0.365691 0.232011 0.109673 0 
0 0.165361 0.322112 0.45157 0.527718 0.532857 0.469307 0.359499 0.233837 0.112851 0 
0 0.161004 0.310853 0.431877 0.502489 0.508517 0.452232 0.352237 0.233652 0.114594 0 
0 0.156035 0.299351 0.413224 0.479345 0.486216 0.435998 0.344221 0.231922 0.115194 0 
0 0.150724 0.287852 0.395539 0.457953 0.465605 0.420492 0.335682 0.229012 0.114892 0 
0 0.145257 0.276517 0.378762 0.438055 0.446411 0.405628 0.326795 0.225209 0.113889 0 
0 0.139762 0.265457 0.362836 0.419447 0.428421 0.391342 0.317692 0.22074 0.112345 0 
0 0.134326 0.254738 0.347712 0.401967 0.411471 0.377585 0.308478 0.215781 0.110392 0 
0 0.129009 0.244403 0.33334 0.385484 0.395431 0.364315 0.299231 0.210474 0.108131 0 
0 0.123845 0.234472 0.319674 0.369893 0.380197 0.351504 0.290014 0.204928 0.105645 0 
0 0.118859 0.224951 0.306672 0.355106 0.365688 0.339126 0.280875 0.199231 0.103001 0 
0 0.114061 0.21584 0.294292 0.341053 0.351837 0.327161 0.271853 0.19345 0.100248 0 
0 0.109457 0.207128 0.282497 0.327671 0.338589 0.315591 0.262977 0.18764 0.0974288 0 
0 0.105047 0.198804 0.271251 0.31491 0.325899 0.304404 0.254269 0.181843 0.0945747 0 
0 0.100827 0.190852 0.260522 0.302724 0.313731 0.293586 0.245747 0.176091 0.091711 0 
0 0.0967928 0.183257 0.250278 0.291075 0.302051 0.283126 0.237423 0.170411 0.0888577 0 
0 0.0929367 0.176002 0.240492 0.279929 0.290832 0.273013 0.229307 0.164822 0.0860302 0 
0 0.0892514 0.169068 0.231137 0.269255 0.280049 0.263238 0.221404 0.15934 0.0832406 0 
0 0.0857291 0.162441 0.222188 0.259027 0.269681 0.253792 0.21372 0.153975 0.0804981 0 
0 0.082362 0.156104 0.213624 0.249219 0.259708 0.244664 0.206255 0.148737 0.0778098 0 
0 0.0791424 0.150042 0.205423 0.239811 0.250112 0.235847 0.199012 0.143632 0.0751808 0 
0 0.0760625 0.144241 0.197567 0.23078 0.240877 0.227332 0.191988 0.138662 0.0726152 0 
0 0.0731154 0.138686 0.190037 0.22211 0.231988 0.21911 0.185182 0.133832 0.0701155 0 
0 0.070294 0.133365 0.182816 0.213783 0.223431 0.211172 0.178592 0.129141 0.0676835 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
