function solving1(u)

options.abstol=1e-5;



[z,info,msg]=fsolve('f',1)


	
endfunction k=f(x)
k=x^2-u
function