wxMaxima script - 4 keer integreren en oplossen

[K.S.]

Het lukt mij niet om een script in wxMaxima kloppend te krijgen. Ik heb ook chatGPT gebruikt om een deel te schrijven. Wellicht dat hier iemand is die mij kan helpen.

In wxMaxima wil ik een formule 4-maal integreren. Vervolgens geef ik 4 randvoorwaarden om de 4 onbekenden op te lossen. Uiteindelijk wil ik dus 4 formules hebben die zijn opgelost.

De basis is alsvolgt:

/* Definieer de functie f(x) */
f(x) := q; /* Vervang dit door de functie die je wilt integreren */

/* Integreer de functie vier keer */
I1(x) := integrate(f(x), x);
I2(x) := integrate(I1(x), x);
I3(x) := integrate(I2(x), x)/EI;
I4(x) := integrate(I3(x), x);

/* Definieer de randvoorwaarden */
v0: subst(x=0, I2) = 0;
v1: subst(x=L, I2) = 0;
v2: subst(x=0, I4) = 0;
v3: subst(x=L, I4) = 0;

/* Los de onbekenden (integratieconstanten) op */
eq1: I1(x0) = v3;
eq2: I2(x0) = v2;
eq3: I3(x0) = v1;
eq4: I4(x0) = v0;

stelsel: solve([eq1, eq2, eq3, eq4], [C, C1, C2, C3]); /* Los het stelsel van vergelijkingen op */

/* Vervang de gevonden constanten in de geintegreerde functies */
I1_final(x) := subst(stelsel, I1(x));
I2_final(x) := subst(stelsel, I2(x));
I3_final(x) := subst(stelsel, I3(x));
I4_final(x) := subst(stelsel, I4(x));

/* Toon de resultaten */
I1_final(x);
I2_final(x);
I3_final(x);
I4_final(x);


De functies worden niet goed opgelost en de gekregen formules kloppen niet.
Om een betere beeldvorming te krijgen: het gaat om het bepalen van functievoorschriften voor een ligger. De dwarskracht, momentenlijn, hoekvedraaiing en zakking wil ik weten. Als ik een basis script heb dan kan ik daarmee de 1e formule opgeven (de belastingslijn) en door de randvoorwaarden aan te passen, pas ik de manier van opleggen aan (vrij, rol, ingeklemd)

Wie O wie kan mij helpen?

[DJMaze]

Kan dit helpen?
https://github.com/wxMaxima-developers/wxmaxima/discussions

[K.S.]

Daar had ik ook zitten kijken. Dat forum is echter niet zeer actief.