function yinit=guess(x)
global lamida
y1=exp(-lamida*x)
y2=0.*y1
yinit=[y1
       y2];