AHSAN

40 Reputation

2 Badges

0 years, 353 days

MaplePrimes Activity


These are replies submitted by AHSAN

@Kitonum 

basically I need relation for x in terms of lamda to use in another relation, which is not possible to find from in my case. so I want to generate data points first and then write a polynomial in term of lamda and then use that polynomial to interpolate my original problem.

I appreciate your help but you write in terms of x, but I need a polynomial in terms of lamda, not x means x is a function of lambda. i hope you will understand my issue and solve my problem

@acer 

basically i need relation for x in terms of lamda to use in another relation, which is not possible to find from in my case. so I want to generate data points first and then write a polynomial in term of lamda and then use that polynomial to interpolate my original problem.

@Carl Love 

can you pls upload its maple file because this is not working...

@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. 

@Carl Love 

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

@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

@tomleslie 

sir thanks allot for your a lot help and cooperation

@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

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 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 

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

 

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

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

when i replaceabove command i could get an error

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.

1 2 3 4 5 6 Page 6 of 6