function z=koule(x,y,R,S1,S2,S3)

z=sqrt(R^2-(x-S1)^2-(y-S2)^2)+S3;