% 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.0383714 0.234833 0.763974 0.955327 0.955327 0.763974 0.234828 0.038344 0.00610243 0 
0 0.095233 0.318309 0.674985 0.877886 0.877896 0.675059 0.318765 0.0980589 0.0225666 0 
0 0.135829 0.352851 0.627891 0.80473 0.804802 0.628338 0.355059 0.14596 0.0442836 0 
0 0.160308 0.366809 0.594039 0.742377 0.742631 0.595427 0.372431 0.180641 0.0650862 0 
0 0.17343 0.369826 0.56519 0.69009 0.690715 0.568232 0.380281 0.204945 0.0823935 0 
0 0.179126 0.366261 0.538818 0.645853 0.647067 0.544196 0.382444 0.221466 0.0957165 0 
0 0.180063 0.358594 0.514096 0.607811 0.609817 0.522338 0.38087 0.232139 0.105413 0 
0 0.17798 0.34841 0.490724 0.574524 0.577478 0.502154 0.37672 0.2384 0.112085 0 
0 0.174014 0.336768 0.468579 0.544927 0.548925 0.483317 0.370756 0.241325 0.11634 0 
0 0.16891 0.324376 0.447591 0.518246 0.523322 0.465586 0.363507 0.241732 0.118707 0 
0 0.163165 0.311709 0.427707 0.493914 0.500048 0.448775 0.35535 0.240246 0.11962 0 
0 0.157104 0.29908 0.408873 0.471514 0.478645 0.432746 0.346565 0.237347 0.119427 0 
0 0.150944 0.286691 0.391038 0.450731 0.458772 0.417392 0.337359 0.233401 0.118397 0 
0 0.144824 0.274667 0.374146 0.431328 0.440175 0.402632 0.327889 0.228689 0.116744 0 
0 0.138835 0.263081 0.358141 0.413123 0.422661 0.388409 0.318275 0.223428 0.114629 0 
0 0.133031 0.251971 0.342971 0.395969 0.406085 0.374678 0.308611 0.217785 0.112178 0 
0 0.127443 0.241351 0.32858 0.379753 0.390335 0.361407 0.298967 0.211886 0.10949 0 
0 0.122088 0.231221 0.31492 0.364381 0.375325 0.348571 0.289399 0.205831 0.106637 0 
0 0.116969 0.221568 0.301942 0.349777 0.360986 0.336153 0.27995 0.199695 0.103679 0 
0 0.112085 0.212376 0.289602 0.335876 0.347263 0.324136 0.270653 0.193538 0.100658 0 
0 0.10743 0.203624 0.277858 0.322624 0.334113 0.312508 0.261533 0.187403 0.0976106 0 
0 0.102997 0.195291 0.266672 0.309975 0.321497 0.301259 0.25261 0.181327 0.0945617 0 
0 0.0987734 0.187352 0.256008 0.297887 0.309384 0.290379 0.243898 0.175336 0.0915321 0 
0 0.0947502 0.179787 0.245833 0.286325 0.297747 0.279859 0.235407 0.169449 0.0885372 0 
0 0.090916 0.172574 0.236119 0.275258 0.286563 0.26969 0.227144 0.163683 0.0855889 0 
0 0.0872603 0.165692 0.226838 0.264658 0.275809 0.259864 0.219114 0.158048 0.0826962 0 
0 0.083773 0.159122 0.217964 0.254498 0.265466 0.250372 0.211319 0.152553 0.0798659 0 
0 0.0804442 0.152845 0.209474 0.244755 0.255516 0.241205 0.203759 0.147203 0.0771028 0 
0 0.0772649 0.146846 0.201347 0.235409 0.245944 0.232355 0.196434 0.142002 0.0744103 0 
0 0.0742265 0.141108 0.193564 0.226439 0.236733 0.223814 0.189342 0.136952 0.0717909 0 
0 0.0713212 0.135617 0.186105 0.217828 0.22787 0.215572 0.182479 0.132054 0.069246 0 
0 0.0685415 0.13036 0.178956 0.209559 0.21934 0.207621 0.175843 0.127307 0.0667762 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
