Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
Download

Interpolation basis

1073 views


�eVW�
�@s#ddlmZddd�ZdS)�)�diffFc)Ks�|r�d|kr>|d}td|�tdt|��d|kr�|d\}}}td|�td|�td|�d|kr�|d\}}}}	}
d	|kr|d	\}}}}
td
|||f�td|||f�td|||
f�d
|kr\|d
\}}}td|�td|�td|�d|kr�|d\}}td|�td|�d|kr�|d\}}td|�td|�d|kr|d\}}td|�td|�d|krF|d\}}}}}td|||||f�d|krj|d}td|�d|kr�|d\}}td|�td|�d|kr�|d\}}}td||f�td ||f�d!|kr"|d!\}}}td"||||f�d#|kr[|d#\}}}td$||||f�d%|kr�|d%\}}td|�td|�d&|kr|d&\}}}}td'�td(|�td'�td)�td|�td|�t|�d*|krC|d*\} }!}"}#td+| |!|"|#f�d,|kr�|d,\}}$}%}&}'}(td|�td-|$�td.|%|&f�td/|'|(f�d0S)1z0Work like:
    dbg(**{"string1":value,...})
    zlagrangeInterp|part1zexpr=%sz
diff(expr)=%szlagrangeInterp|part2zlen(c[0][:-1])=%sza=%szd[0]=%szlagrangeInterp|part3zlagrangeInterp|part4zd[%s-1][:-%s]=%szd[%s-2][:-(%s-1)]=%szd[%s-2][:-%s]=%szlagrangeInterp|part5zc=%szd=%szlagrangeInterp|part6zA=%szX=%szlagrangeInterp|part7z
f(a[0])=%sz
exprSum=%szdiffSocket|part1z	points=%szdiffSocket|part21z0b[:-1][%s]=a[:-1][%s]+(x0-points[%s])*b[(%s)]=%szdiffSocket|part22zdiffSocket|part3zb=%szderivations=%szifEqPoints|part1zany c[%s-1]=%sz
d[%s-1][:]=%szifEqPoints|part21zif d[%s-1][:-%s][%s]=%szifEqPoints|part22zelse d[%s-1][:-%s][%s]=%szleastSquares|part1zleastSquares|part2z************zHH=%sz
**********zfindParabs|part1za0=%s|a1=%s|b0=%s|b1=%szfindParabs|part2zinvC=%szp0,p1= %s,%szc0,c1=%s %sN)�printr))�show�part�expr�n�a�d0�i�jZcijZaijZdeltaAijZdij1Zdij2Zdij3�c�d�A�XZfa0�exprSum�points�arg0�arg1�arg2�arg3�arg4�b�deriv�ci�diZdji�HH�a0�a1�b0�b1ZinvC�p0�p1�c0�c1�r$�A/projects/91f280f0-b3af-42e1-a2e0-ed2660b96906/jupyter/src/dbg.py�dbgs�





r&N)�sympyrr&r$r$r$r%�<module>s