fuction c=func1(x)
c=(x./1.6).*(x>0.6)+(x.^2+0.3).*(x<=0.6);