def Lmodtestm(lam,i, p):
x0=np.array([56.,100.,760.,35.,9888.])
dx=np.array([5.6,10.,76.,3.5,988.8])
dxlam=lam*dx
x0h=x0+dxlam
mxdx=DAm.nldalec(x0h,i)[1][i-1]
mx=DAm.nldalec(x0,i)[1][i-1]
#Mlist=DAm.ldalec(x0,i)[0]
#Mdx=DAC3.Mfac(Mlist,i)*np.matrix(dxlam).T
Mdx=DAm.ldalec2(x0,dxlam,i)[i-1]
if p=='all':
err=(np.linalg.norm(mxdx-mx))/(np.linalg.norm(Mdx))
else:
err=(np.linalg.norm(mxdx[p]-mx[p]))/(np.linalg.norm(Mdx[p]))
return err
