function y=MISE(hh,m,v,n,I,K)
% function y=MISE(hh,m,v,I,K)

if nargin<5
 error('nwmatrix: Too litle input parameters, sorry.');
 return;
end
if nargin==5
 K=K_def('epan');
end
k=K.k;
ms=sym(m);
md=diff(ms,k);
Ims=int(md^2,I(1),I(2));
Im=eval(Ims);
vs=sym(v);
Ivs=int(vs,I(1),I(2));
Iv=eval(Ivs);
bias2=1/k^2*Im*K.beta^2*hh.^(2*k);
var=K.var*Iv/n*1./hh;
y=bias2+var;
pl=plot(hh,bias2,hh,var,hh,y);
set(pl,'Linewidth',2);