MaplePrimes
  • About
  • Help
  • Support
  • Register
  • Login
Advanced Search
  • Questions
  • Posts
  • Products
  • Users
  • Unanswered
  • Maplesoft Blog
  • Badges
  • Recent
  • More
    • Tags
    • RSS Feeds
    • Submit Software Change Request
    • maplesoft.com

Ask a Question

Create a Post

Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer
Unanswered Questions Feed
Edit

Where the error in the following code?...

Asked by: khalil20232023 15
homework differential-equations duplicate-question + Manage Tags
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

Edit

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

Asked by: Rouben Rostamian 7501
ode piecewise dsolve + Manage Tags