load wvdata
load weight

global weight tdatatrue tdataret

Tnoise=0.1

tdata=[300 260 340 280 260 260 260 260 260 260 260];

tdatatrue=tdata;

tdatag=[300*ones(1,11)];

Idata=ftoa(tdata,wvdata);



for l=1:20
TBdata=TB(Idata(l),wvdata(l));
TBdatanoise=TBdata+Tnoise*randn(1,1);
Idata_noise(l)=planck(wvdata(l),TBdatanoise);
end

plot(wvdata,Idata,wvdata,Idata_noise)








options = optimset('TolFun',1e-20,'MaxFunEvals',10000,'MaxIter',100);

tdata_ret=LSQCURVEFIT('ftoa',tdatag,wvdata,Idata_noise,[],[],options);


zv=[0:10];

plot(zv,tdata,zv,tdata_ret)