一个函数与其导数的图象绘制

clc
clear
close
format long
syms x y
t0=1.001;
t=-1:0.001:5;
x1=t;
y1=t.*(t-1).^(2/5);
z1=t*0;
plot(x1,y1);
hold on
Rx0=t0;
Ry0=t0*(t0-1)^(2/5);
plot(Rx0,Ry0,'*');

grid on
axis equal

df=diff(x);
dg=diff(x*(x-1)^(2/5));
df2=diff(df);
dg2=diff(dg);
k=abs(df*dg2-dg*df2) / (df^2+dg^2)^(3/2);
% pretty(k)
kt0=subs(k,t0);
df0=subs(df,t0);
dg0=subs(dg,t0);
Nx0=dg0/sqrt(dg0^2+df0^2);
Ny0=-df0/sqrt(dg0^2+df0^2);

Cx0=Rx0 + (1/kt0) *Nx0;
Cy0=Ry0 + (1/kt0) *Ny0;
plot(Cx0,Cy0,'*')

x2=Cx0+cos(t)*(1/kt0);
y2=Cy0+sin(t)*(1/kt0);
plot(x2,y2);
double(kt0)
%n=-dg i+df j
%n=dg i- df j

% eval(solve('3*t/sqrt(1+t^2)=2.281'))