Advanced Search

AHSAN

160 Reputation

6 Badges

5 years, 133 days
Contact AHSAN

MaplePrimes Activity


These are replies submitted by AHSAN

ODE solution...

AHSAN 160
May 01 2020
0 0

@Carl Love

I am interested in any numerical technique to develop interest e.g false position method, newton raphan method, etc rather than in built-in command like fsolve. 

ODE solution...

AHSAN 160
April 30 2020
0 0

@Carl Love 

k is a constant and fixed it valve  k=0.1 then solve with respect to lamda as a variable.

ODE solution...

AHSAN 160
April 28 2020
0 0

@tomleslie 

sir can you pls have a look on it .. I am very beginner to maple so that why face a lot of problems
 

> 

restart
> 

newton := proc (g, X__A, X__B, err) local j, p, X__1, X__2, iters; X__1 := X__A; X__2 := X__B; iters := 0; for j to 50 do iters := iters+1; p := evalf(X__1-(eval(g(X__1)/(diff(g(x), x)), x = X__1))); if abs(g(p)) < err then break else X__1 := p end if end do; return root = p, residual = g(p), numberOfIterations = iters end proc
> 

F := proc (Q) options operator, arrow; 0.2185844281e-2*(428.8960600*sqrt(2)*(0.127600000e-1-(6/5)*Q)*arctan((1/2)*sqrt(2.42704/(2*Q+.99)-2)*sqrt(2))/(2*Q+.99)^3+673.7083558*(0.127600000e-1-(6/5)*Q)*sqrt(2)/(2*Q+.99)^3+60*sqrt(2.42704/(2*Q+.99)-2)*((0.127600000e-1-(6/5)*Q)*(2.42704/(2*Q+.99)-2)^2+(0.68053333e-1-(32/5)*Q)*(2.42704/(2*Q+.99)-2)+1.379488000-8*Q))*(2*Q+.99)^3 end proc

proc (Q) options operator, arrow; 0.2185844281e-2*(428.8960600*sqrt(2)*(0.127600000e-1-(6/5)*Q)*arctan((1/2)*sqrt(2.42704/(2*Q+.99)-2)*sqrt(2))/(2*Q+.99)^3+673.7083558*(0.127600000e-1-(6/5)*Q)*sqrt(2)/(2*Q+.99)^3+60*sqrt(2.42704/(2*Q+.99)-2)*((0.127600000e-1-(6/5)*Q)*(2.42704/(2*Q+.99)-2)^2+(0.68053333e-1-(32/5)*Q)*(2.42704/(2*Q+.99)-2)+1.379488000-8*Q))*(2*Q+.99)^3 end proc

 (1)
> 

``
> 

``


 

Download try1.mw

ODE solution...

AHSAN 160
April 28 2020
0 0

@tomleslie 

sir thanks allot for your a lot help and cooperation

ODE solution...

AHSAN 160
April 27 2020
0 0

@tomleslie 

it's really good for me and thanks a lot sir . can you pls told me the name of this approach means regular false method, newton method, etc

ODE solution...

AHSAN 160
April 27 2020
0 0

many regard sir for your help and cooperation.and the problem you fixed with the student(precalculus)is it possible to solve same problem by using a false position method? @tomleslie 

try.mw   can you please have l...

AHSAN 160
April 26 2020
0 0

try.mw   can you please have look on this file, which I have loaded. I am trying to solve it but could not get the result. and is it possible,you'll solve this problem with(Student[NumericalAnalysis]) false position method @tomleslie 

ODE solution...

AHSAN 160
April 26 2020
0 0

i am trying to solve this equation but could not get the result. need help

> 

restart
> 

with(Student[Precalculus])
> 

k := .1

.1

 (1)
> 

F := (-(288*(lambda+2*k*(1/3)-2/3))*sqrt(2)*lambda^2*arctan((1/2)*sqrt(4*lambda/(1-k)-2)*sqrt(2))/(1-k)^2-(144*(lambda+2*k*(1/3)-2/3))*Pi*lambda^2*sqrt(2)/(1-k)^2-(36*((lambda+2*k*(1/3)-2/3)*(4*lambda/(1-k)-2)+10*lambda*(1/3)+4*k*(1/3)-4/3))*sqrt(4*lambda/(1-k)-2))*(1-k)^2/(256*lambda^2)

0.3164062500e-2*(-355.5555555*(lambda-.6000000000)*2^(1/2)*lambda^2*arctan((1/2)*(4.444444444*lambda-2)^(1/2)*2^(1/2))-558.5053607*(lambda-.6000000000)*lambda^2*2^(1/2)-36*((lambda-.6000000000)*(4.444444444*lambda-2)+(10/3)*lambda-1.200000000)*(4.444444444*lambda-2)^(1/2))/lambda^2

 (2)
> 

`&lambda;__1` := 0

0

 (3)
> 

`&lambda;__2` := 1

1

 (4)
> 

tol := 0.10e-3

0.10e-3

 (5)
> 

for j do p := evalF([Line([`&lambda;__1`, F(`&lambda;__1`)], [`&lambda;__2`, F(`&lambda;__2`)])][5]); if abs(F(p)) < tol then break elif F(p) > 0 then `&lambda;__2` := p else `&lambda;__1` := p end if end do; root = p, residual = F(p)

Error, invalid subscript selector

  

root = p, residual = 0.3164062500e-2*(-355.5555555*(lambda(p)-.6000000000)*2^(1/2)*lambda(p)^2*(arctan((1/2)*(4.444444444*lambda-2)^(1/2)*2^(1/2)))(p)-558.5053607*(lambda(p)-.6000000000)*lambda(p)^2*2^(1/2)-36*((lambda(p)-.6000000000)*(4.444444444*lambda(p)-2)+(10/3)*lambda(p)-1.200000000)*(4.444444444*lambda(p)-2)^(1/2))/lambda(p)^2

 (6)
> 

``
> 

``


 

Download try.mw

 

ODE solution...

AHSAN 160
April 26 2020
0 0

I really appreciated your help, can you pls explain this command 

solve
  ( odetest
    ( genSol, [ode, bc])[2..3],
    [_C1, _C2]
  )[];

why you are taking [2,3] or its meaning

ODE solution...

AHSAN 160
April 26 2020
0 0

cVals:=solve(otest[-2..-1]); 

when i replaceabove command i could get an error

ODE solution...

AHSAN 160
April 25 2020
0 0

dear sir many thanks I have little confusion can you explain the second portion of your answer by talking some second order ode with BCS.

First 16 17 18 Page 18 of 18