f := (x,y) -> sin(sqrt(x))*log(y^2); fx := diff(f(x,y),x); fy := diff(f(x,y),y); D[1](f)(1.,0.5); Dfy := unapply(fy,(x,y)); Dfy(1.,0.5);