<![CDATA[
#include <math.h>
#include <chplot.h>
#include <numeric.h>
int main() {
int i, j;
array double s[3][3]={{1,2,1},{2,4,2},{1,2,1}};
array double x[16],y[16],z1[256],z[16][16],Z[18][18],Z1[256];
linspace(x,0,16);
linspace(y,0,16);
for(i=3; i<13; i++)
for(j=3; j<13; j++) {
z1[i*16+j]=1;
z[i][j] = 1;
}
plotxyz(x,y,z1);
filter2(Z,z,s,"full");
for(i = 0; i<16; i++)
for(j=0; j<16; j++)
Z1[i*16+j] = Z[i+1][j+1];
plotxyz(x,y,Z1);
}
]]>
