Write code to implement the Hermite interpolation algorithm and approximate

the function $f(x)=tanh(x)$ at $x_i=i,i=0,$dots,$20$ using $p_{41}(x)$ such that

$p_{41}(x_i)=f(x_i),p_{41}^{\text{'}}(x_i)=f^{\text{'}}(x_i).$ Make a plot.