I've got my TI Nspire CX CAS today and I'm really surprised in a positive way how easy it is to write programs for it.

However, I've got a problem here. I'm writing some tools that should make solving mechanical problems a bit easier. One of these mechanical problems regards Euler-Bernoulli beam theory. The dynamics of a beam are described by a special function and its differentials…. well I'm not going to go to deep into details here. I'm just looking for a way to create a program that asks for the constraints like this:

Enter order of differentiation n=…

Enter position of constraint x=…

Enter value of constraint f^{(n)}(x)=…

Looping through this requests until n < 0 (I'm not gonna take account of cases where n is indeed < 0) and recording every given constraint.

With these given constraints, solving the system of equations (the given f(x) and its n differentials) for all required coefficients should be possible. At least manually using solve().

**Now, my problem is**, I don't have the slightest clue on how to get from the constraint matrix that might look like this

ca:=[[0, 0, 0][1, 0, 0][0, 500, 0]]

to a TI Nspire Basic command that somehow looks like this

solve({ f^{(0)}(0)=0, f^{(1)}(0)=0, f^{(0)}(500)=0}, {a, b, c, d})

where f(x) is for example

f(x):=ax^{3}+bx^{2}+cx+d

Since the system consists of three equations for four unknowns I can define *a* for example to be *1* (a valid method to get the coefficients in case of a Euler-Bernoulli beam).

So how do I convert the matrix to a solve command?

Cheers,

Hendrik