function [y1] = reuna_arvo(dy0)
%  syote: derivaatta-arvo
%  tulos: y(1) 
[t y]=ode45(@shoot,[0 0.5 1],[0;dy0]);
y1=y(end,1);
end

