Advanced Search

Syed Asadullah Shah

15 Reputation

One Badge

5 years, 180 days
Contact Syed Asadullah Shah

MaplePrimes Activity


These are questions asked by Syed Asadullah Shah

What is the error in my code (DTM Pade ...

Syed Asadullah Shah 15 Maple 2019
March 06 2023
1 0

HI, I attach  a document (Pade_Approximants_for_alpha.mw) extract showing the errors (Error, (in Engine:-Dispatch) badly formed input to solve: not fully algebraic
).

Can any one help us please?

Thanks. Syed Asadullah Shah

Please help with this DTM-Pade CODE...

Syed Asadullah Shah 15 Maple 2016
solve syntax homework
October 11 2022
1 7

DTM_Pade_for_A_and_B.mw

> 

restart; Digits := 10; m := 2; F[0] := 0; M := 1; F[1] := 1; F[2] := (1/2)*A; T[0] := 0; T[1] := B; for k from 0 to m+2 do F[k+3] := (M*(k+1)*F[k+1]+sum((k-r+1)*(r+1)*F[r+1]*F[k-r+1], r = 0 .. k)-(sum((k-r+1)*(k-r+2)*F[r]*F[k-r+2], r = 0 .. k)))*factorial(k)/factorial(k+3); T[k+1] := (2*(sum((k-r+1)*F[k-r+1], r = 0 .. k))+(2*(k+2))*(k+1)*F(k+2))*factorial(k)/factorial(k+1) end do; f := 0; t := 0; for k from 0 to m do f := f+F[k]*x^k; t := t+T[k]*x^k end do; print(f); print(t); with(numapprox); pade(t, x, [2, 2]); pade(diff(f, x), x, [2, 2]); solve({limit(pade*(t, x, [2, 2]), x = infinity) = 0., limit(pade(diff(f, x), x, [2, 2]), x = infinity) = 0.}, [A, B])

x+(1/2)*A*x^2

  

(2+4*(table( [( 0 ) = 0, ( 1 ) = 1, ( 2 ) = (1/2)*A, ( 3 ) = 1/3, ( 4 ) = (1/12)*A, ( 5 ) = 1/60+(1/120)*A^2, ( 6 ) = (1/180)*A, ( 7 ) = 1/840+(1/5040)*A^2 ] ))(k+2))*x+(1+A+6*(table( [( 0 ) = 0, ( 1 ) = 1, ( 2 ) = (1/2)*A, ( 3 ) = 1/3, ( 4 ) = (1/12)*A, ( 5 ) = 1/60+(1/120)*A^2, ( 6 ) = (1/180)*A, ( 7 ) = 1/840+(1/5040)*A^2 ] ))(k+2))*x^2

  

Error, (in testeq) invalid arguments

  
> 

NULL

Download DTM_Pade_for_A_and_B.mw

Differential Transform Method (DTM) ...

Syed Asadullah Shah 15 Maple
ode homework numeric
September 23 2022
1 3

Hi

I have an equation which i want solve it with DTM and sketch the graph of f(η) along y axis and η along x axis, f ‘( η) and g( η)

 

 

M:=0..5

F[0]:=0;F[1]:=1;F[2]:=a/2;G(0)=0

That's my pleasure If there is somebody help me to code and solve this problem.

Thanks for every one very very much.

Page 1 of 1