Numan W

55 Reputation

6 Badges

6 years, 128 days

MaplePrimes Activity


These are replies submitted by Numan W

Thank you so much @dharr 

Thank you @Mariusz Iwaniuk

 

Thank you @Carl Love 

I am exciting to see your answer

Thank you @Preben Alsholm , I like your answer .

But for other equation the singularity will show up , and we can’t change the form of equation , so for that if it's possible  I  need a conditional if to re-solve the equation for example if it's a singularity must return and solve it ( dsolve ) and so that , I hope I have clarified the situation.

Algorithm

Random coefficient

Dsolve

Calculate eng

If eng = singularity

Return and dsolve

Else, end if

 

thank you 

Thank you @J4James ,

but i steal have probleme whene i change the last condition from 70 to 100 

psi(0.174568e-1) = 100, i get the flowing error ;

Error, (in dsolve/numeric/ComputeSolution) Newton iteration is not converging
 

best regards 

thank you Mr.@Preben Alsholm 

i will look for other way 

Best regards!

 

Thank you for your answer Mr.@tomleslie 

i defined the parameters but the problem still not solved 


 

``

restart:with(inttrans):

E:=1;L:=1;k:=1;R:=1;m:=1;

1

 

1

 

1

 

1

 

1

(1)

θl:=arccos(y(t)/R)

arccos(y(t))

(2)

equ := arccos(y(t)/R)*R*L*(diff(y(t), `$`(t, 1)))*abs(diff(y(t), `$`(t, 1)))+diff(y(t), `$`(t, 2))+m*sin(E*t+k*R*sin(`θl`));

arccos(y(t))*(diff(y(t), t))*abs(diff(y(t), t))+diff(diff(y(t), t), t)+sin(t+(1-y(t)^2)^(1/2))

(3)

eq:=fourier(equ,t,w);

fourier(sin(t+(1-y(t)^2)^(1/2)), t, w)+fourier(arccos(y(t))*(diff(y(t), t))*abs(diff(y(t), t)), t, w)-w^2*fourier(y(t), t, w)

(4)

csi := y(0) = 0.2e-1, (D(y))(0) = 0;

y(0) = 0.2e-1, (D(y))(0) = 0

(5)

sol := dsolve({csi, eq}, numeric, maxfun = 1000000000)

Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)

 

Error, (in solve) cannot solve expressions with fourier(arccos(Y[1])*YP[1]*abs(YP[1]), t, w) for YP[1]

 


 

Download Fourier_TRAns_MAPLEprime.mw

Thank you

 

thank you for your answer , but it's not work with dsolve 

Hi

Good idea 

thanks :)

Thank you very much 

how i can plot these solution ?

thank you it's work

thank you 

but theres other tools without LambertW

thank you Preben Alsholm  it's work 

Thank you so much Thomas Richard & thank you for your advice I appreciate that.

 

 

1 2 Page 1 of 2