function h=forma(x)

  global c q m;

h(1)=0.01*x^2+c*sin(x/4);
h(2)=0.01*2*x+c*cos(x/4)/4;
h(3)=0.01*2-c*sin(x/4)/16;

end