Tamour_Zubair

50 Reputation

One Badge

0 years, 174 days

MaplePrimes Activity


These are replies submitted by Tamour_Zubair

@tomleslie 

Dear Professor;

Hope you are fine. You are absolutly right Professor. Maple builton commands are always better and I always like it. In my situation, when i increase the value of "N", number of nonlinear ODE's increased and further builton comand take a lot of computational time to evaluate the solution. But my target is reduce the computation cost. I have upload the real code. I have two generlized ODE's and further with the help of loops and discritization on different points with value of "N", i obtained system of nonlinear ODE's and further i need to solve with the help of RK-4. For better accuracy, I need to increase "N" which increse the number of nonlinear ODE's. Kindly visit my real code and help me further. I am waiting for your kind response.

Thanks 

@tomleslie 

Dear Professor;

Thanks a lot for this. I cannot use builton command beacuse system of ODE's increase as i will increase the value of "N". Therefore this command will creat heavy computational cost. Is that possible to make code not builton command? I am waiting for your kind response.

Thanks

@vv 

Thanks a lot. It helps me a lot. Can it take first ODE of each variable ?

@Carl Love 

Professor Z_0(t) is extra which i have removed and update the data and file accordingly. Z_1(t) and Z_2(t) is not extra, it contaion odes and part of the system of ODE's. System is nonlinear due this constants. Without this system become linear which can easily handle with mentioed command but i need to convert system of nonlinear ODE's to matrixcs.

@Rouben Rostamian

Yes Professor.

I need matrix formulation of nonlinear ODE's. Is that any solution? 

@Axel Vogt 

Dear Professor;

Thanks a lot for your response. Saperately it is solving but not collectively. I need to find the solution in term of other variable. Not possible to find?

@acer 

Thanks a lot Professor. I have updated the file name according to your instruction. Now it is ok. Thnaks a lot for your response.

I have uploded the file properly but i donot know why maple file download is showing error.

@Carl Love 

Thanks a lot Professor for your kind response. The coding and scheme style you used totally different and may be very hard for me. I have some question about it.

1. How you formulate M and further pi and sys. Can you please explain the steps.

2. What method you used find the error?

I am waiting for your kind response.

Thanks

@tomleslie

Thanks a lot Professor for your kind response. Can we say that our scheme and coding is completely accurate? 

@tomleslie 

Really thanks to you Professor. Its very helpful.

With you given ideas I devolped the code of rk-2, rk-3. But feel that there is no acccurcy difference between rk-2 and rk-4. I checked builton commands that rk-2 and rk-4 has same reults. Is that any issue in scheme? I am waiting for your kind response.

Thanks 

@rcorless 

Sure Professor.

Kindly share the link of that paper. I will explore it. I shall be very greatfull to you. Thanks

@acer 

Professor i want to evaluate this

question11.mw

restart

``

JJx11 := int(int(2.73299999999999999999999999999*10^(-22)*p*cos(.50*p+.50*q)/sqrt(p^2+q^2+7.46928900000000000000000000000*10^(-44)), p), q)

int(int(0.2733000000e-21*p*cos(.50*p+.50*q)/(p^2+q^2+0.7469289000e-43)^(1/2), p), q)

(1)

``

Download question11.mw

@acer 

Professor we the simplest form of my question is given below. I need to find "AA22" in integration evaluated form at p=1 and q=1.

question.mw

restart; printlevel := 4

{--> enter Terminate, args =

<-- exit Terminate (now at top level) = }
{--> enter ModuleUnload, args =
<-- exit ModuleUnload (now at top level) = }

 

with(LinearAlgebra):

``

ff11 := cos(1.0*k*(p+q)):

NULL

JJx11 := int(int(p*ff11/sqrt(p^2+q^2), p), q);

int(int(p*cos(.50*p+.50*q)/(p^2+q^2)^(1/2), p), q)

(1)

AA22 := eval(JJx11, [p = 1, q = 1])

eval(int(eval(int(p*cos(.50*p+.50*q)/(p^2+q^2)^(1/2), p), {p = 1}), q), {q = 1})

(2)
   
 

``

Download question.mw

Page 1 of 1