ans=[4145306.61         512963.80         18.96 ;

4145365.17         511893.37         19.28 ;

4145272.03         512127.81         22.11 ;

4145362.48         512051.28        20.83 ;

4145420.53         511982.28         19.30 ;

4145258.77         512237.06        23.42 ;

4145299.82         512187.76        22.63 ;

4145341.06         512138.86        21.78 ;

4145138.16         512474.12        25.10 ;

4145200.09         512399.98        23.99 ;

4145224.39         512371.23        23.83 ;

4145269.97         512316.88        23.59 ;

4145340.64         512232.51        22.85 ;

4145386.51         512178.34        21.97 ;

4145443.66         512109.8                        21.05 ;

4145275.35         512388.29        23.52 ;

4145300.38         512350.25        23.22 ;

4145335.95         512315.68        22.89;

4145309.43         512457.3                         23.44;

4145313.50         512419.28        23.31 ;

4145346.68         512382.49        22.92 ;

4145367.85         512387.46        23.00 ;

4145400.56         512316.84        22.55 ;

4145452.73         512254.4                          21.78 ;

4145510.78         512185.52        20.91 ;

4145567.15         512118.33        20.16 ;

4145337.82         512607.48        23.63 ;

4145401.81         512533.76        23.32 ;

4145392.04         512481.98        20.67 ;

4145444.98         512419.39        22.50 ;

4145493.02         512361.46        21.69 ;

4145556.40         512443.45         21.78 ;]

X=ans(:,1)

Y=ans(:,2)

Z=ans(:,3)

X1=min(X)

X2=max(X)

X3=X1:10:X2;

Y1=min(Y);

Y2=max(Y);

Y3=Y1:10:Y2;

[X4,Y4]=meshgrid(X3,Y3)

Z4=griddata(X, Y, Z, X4, Y4)

mesh(X4, Y4 , Z4)