x^2=4*p*(y-k);

solve (%,y);

y(x):=(x^2+4*k*p)/(4*p);

fp:diff(y(x),x);

fpp:diff(y(x),x,2);

k:abs(fpp/(1+fp^2)^(3/2));

rho:1/k;