Hnx

20 Reputation

One Badge

5 years, 154 days

MaplePrimes Activity


These are questions asked by Hnx

I have a non-linear deq of 2nd order which I want to solve numerically. The integration should stop if the integrated variable exceeds a certain value, i.e. phi(t)>phi_ end or equivalently phi(t)-phi_end>0. Maple doesn't accept this type of event and the test phi(t)=phi_end doesn't catch (obviously). How to work around?

NULL

Programm zur Brechnung der Kreisteilungspolynome_2022-04-05 Ki

 

NULL

restart; with(Algebraic); with(NumberTheory)

n := 6

6

(1)

 

cyclo_poly := Vector[row](1 .. n, 0)

Vector[row](%id = 36893489967244295332)

(2)

i := 1; cyclo_poly[1] := X-1

1

 

X-1

(3)

NULL
c_poly := proc (i, kt_poly) local j, hz1, hz2, X; j := i+1; hz1 := X^j-1; hz2 := kt_poly; kt_poly := Quotient(hz1, hz2, X) end proc

cyclo_poly[2] := c_poly(i, cyclo_poly[i])

Error, (in c_poly) illegal use of a formal parameter

 

i := 1; kt_poly := X-1

1

 

X-1

(4)

j := i+1; hz1 := X^j-1; hz2 := kt_poly; kt_poly := Quotient(hz1, hz2, X)

2

 

X^2-1

 

X-1

 

X+1

(5)

NULL

Download Programm_zur_Brechnung_der_Kreisteilungspolynome_2022-04-05_Ki.mw

Good morning

Why doesn't Quotient(x,y,z) not work within a procedure whereas it works well standing alone? See example.

Greetings  Heinz

Is there a possibility to test the result of a numerically solved deq?

1. The attached file is a small test sample out of a larger program. I do not understand what is wrong with the code.

2. What does 'global variables are deprecated' mean when I want to solve a deq numerically? Cannot acces link indicated for further explanation.

3. Is there a means to use 'odetest' when a deq was solved numerically?

I have the numeric solution of differential equation and interpolate the data points with a spline of third degree. Is there a way of differentiating the spline as a whole? Differentiating every single polynom seems awkward.

1 2 Page 1 of 2