function [x,y,STDY]=xydata(k_noise)
%xydata.m		產生仿真數據
x=[0:0.2:4]';
yo=3*exp(-0.4*x)+12*exp(-3.2*x);
rand('seed',234)
y_noise=k_noise*(rand(size(x))-0.5);
y=yo+y_noise;
STDY=std(y_noise);