Unanswered Questions

Ask a Question

Create a Post

This page lists MaplePrimes questions that have not yet received an answer

Unanswered Questions Feed

Where the error in the following code?...

Asked by: khalil20232023 15

April 15 2023

1 0

$alpha1 := .8; alpha2 := .6; eq1 := fracdiff(x(t), t, alpha1)-t^2*x(t)+y(t)-2*t = 0; eq2 := fracdiff(y(t), t, alpha2)-2*z(t)+2*(t+1) = 0; eq3 := z(t)-y(t)-2*t*x(t)+t^4-t-1 = 0; N := 4; x0 := 0; y0 := 0; z0 := 1; XT1 := D[xt](eq1); XT2 := D[xt](eq2); XT3 := D[xt](eq3); XT1[0] := x0; XT2[0] := y0; XT3[0] := z0; for n to N do XT1[n] := unapply(simplify(rhs(subs({seq(D[xt](x(t)), j = 0 .. n-1), seq(D[xt](y(t)), j = 0 .. n-1), seq(D[xt](z(t)), j = 0 .. n-1)}, XT1))), t); XT2[n] := unapply(simplify(rhs(subs({seq(D[xt](x(t)), j = 0 .. n-1), seq(D[xt](y(t)), j = 0 .. n-1), seq(D[xt](z(t)), j = 0 .. n-1)}, XT2))), t); XT3[n] := unapply(simplify(rhs(subs({seq(D[xt](x(t)), j = 0 .. n-1), seq(D[xt](y(t)), j = 0 .. n-1), seq(D[xt](z(t)), j = 0 .. n-1)}, XT3))), t) end do; xseries := add(XT1[n](t)*t^(n*alpha1)/factorial(n*alpha1), n = 0 .. N); yseries := add(XT2[n](t)*t^(n*alpha2)/factorial(n*alpha2), n = 0 .. N); zseries := add(XT3[n](t)*t^(n*alpha2)/factorial(n*alpha2), n = 0 .. N); x_sol := unapply(xseries, t); y_sol := unapply(yseries, t); z_sol := unapply(zseries, t); print(x_sol); print(y_sol); print(z_sol)$

proc (t) options operator, arrow; 1.000000000+1.119174954*XT3[1](t)*t^.6+.9076036841*XT3[2](t)*t^1.2+.5964840410*XT3[3](t)*t^1.8+.3354346720*XT3[4](t)*t^2.4 end proc

(1)

Download Ex2_FDAEs.mw

How to find a symbolic solution for this piecewise...

Asked by: Rouben Rostamian 7501