(%i1) rhs(first(solve(x^2+d*y*4+d^2*(-4),y)))$
(%i2) [%,subst(x=-x,diff(%,x)), diff(%,x)]$
Define a function larry that returns a list [y(x,d), dy/dx (-x,d), dy/dx(x,d)]
(%i3) larry(x,d) := ''%$
Check one value
(%i4) larry(10,5);
(%o4) [0,1,-1]
Generate a table (matrix) of values:
(%i5) funmake('matrix, map(lambda([x], append(x, float(apply(larry,x)))), [[10,5], [10,10],[50,25], [50,50],[100,50],[100,100]]));
(%o5) matrix([10,5,0.0,1.0,-1.0],[10,10,7.5,0.5,-0.5],[50,25,0.0,1.0,-1.0],[50,50,37.5,0.5,-0.5],[100,50,0.0,1.0,-1.0], [100,100,75.0,0.5,-0.5])
--bw