#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); }