Differentiation Without Infinitesimals

This shows that infinite techniques are not required for computer calculus.
The program below exactly computes the relative change in x^n with respect to a change in x as a finite digital calculator could do.
The normal formula for derivative is followed except x is scaled up by t so a change in x of 1 is relatively small. x can be scaled up as long as it is not zero.

    If changes(h) in x are continually reduced in relative size in a floating point operation, a size is reached when h+x =x exactly for x<>0, so (f(x+h)-f(x))/h=0/0. So the usual humanly imagined limit process is impossible in a real digital calculator. But 'the scale up' works for non zero x which require relatively few digits compared to calculator precision. Below, the number of digits typed for w must be small compared to the scale t so the didits for various terms of the polynomial expansion are separated by zeros, and the whole string of digits is within the browser's precision.

To stay close to binary logic I use rational numbers: x=w/s and w,s,t,n are integers.
Derivative of x^n ={ (w/s+ 1/(s*t) )^n - (w/s)^n}/(1/(s*t)).
    Factoring 1/(s*t) to outside, { (t*w + 1) +1)^n - (w*t)^n }*  [(1/(s*t))^n] / [1/(s*t)] .
Traditionally we 'find the limit as t is unbounded'. But in an integer calculator we only need t large enough to separate, with zeros, the strings of digits of the products.
To work with only positive integers, I reorder the divisions until they can be applied to sufficiently large numbers. Also t ,s for decimal work are to be a power of ten. This only limits x to finite approximations of infinite decimals. For example: 1/3 is approximated as close as desired well within calculator accuracy. 'well within' is to allow 0's between significant digits in the differentiation process. That is 1/3 is replaced by 333000/1000000=.333000 or if the calculator allows the desired powers of 333300001, use 33330000/100000000.
So we have { (t*w + 1)^n - (t*w)^n }/(s*t)^(n -1) = (Z^n-z^n)/(s*t)^(n-1) for z=t*w, Z=z+1.
This is very general since most physical functions can be generated or approximated by polynomials whose derivative is a sum of derivatives of each term.

scale t= w=. Z=t*w+1= = initial value below. z=t*w
n= Z^n= 
Z^n-z^n = 
(Z^n-z^n)/t^(n-1) = 
/s^(n-1) = traditional derivative of x^n = n*x^(n-1) = 
Here keeping in the integer calculator mode, we only externally represent division by s^(n-1).

Discussion


E-mail author: R Massey
top/home (25k)