def u(x,y): return x*y; def du_x(x,y): return y; def du_y(x,y): return x; def v(x,y): return 1 + pow(x,2) + pow(y,2); def dv_x(x,y): return 2*x; def dv_y(x,y): return 2*y; def f(x,y): return u(x,y)/float(v(x,y)); def df_x(x,y): return (du_x(x,y)*v(x,y) - u(x,y)*dv_x(x,y))/float(pow(v(x,y),2)); def df_y(x,y): return (du_y(x,y)*v(x,y) - u(x,y)*dv_y(x,y))/float(pow(v(x,y),2)); def h(x,y): return (y - pow(x,2)*y + pow(y,3))/float(pow(1 + pow(x,2) + pow(y,2),2)); x = 1; while(x<15): print h(x,1)-df_x(x,1); x +=1;