Advanced Search

wimvr

10 Reputation

3 Badges

13 years, 50 days
Contact wimvr

MaplePrimes Activity


These are replies submitted by wimvr

My own answer...

wimvr 10
May 14 2012
0 0

Hey!

Thanks for the quick reply Markiyan Hirnyk!

I also found another solution (maybe a bit slower) but also usefull:

> stlsl := {-(40.00000000*(0.4469633000e-3*(diff(y(t), t))^2+0.4469633000e-3*(diff(x(t), t))^2))*(diff(x(t), t))/sqrt((diff(y(t), t))^2+(diff(x(t), t))^2) = diff(x(t), t, t), (40.00000000*(0.4469633000e-3*(diff(y(t), t))^2+0.4469633000e-3*(diff(x(t), t))^2))*(diff(y(t), t))/sqrt((diff(y(t), t))^2+(diff(x(t), t))^2)-9.810000000 = diff(y(t), t, t), x(0) = 0, y(0) = .80, (D(x))(0) = 8.5*cos(0.1744444444e-1*angle), (D(y))(0) = 8.5*sin(0.1744444444e-1*angle)}

>sol := dsolve(eval(stlsl,angle = angeltest), numeric):
>subs(opl(4), x(t));

Using a double for loop, I can do exactly the same ;-)

My own answer...

wimvr 10
May 14 2012

Hey!

Thanks for the quick reply Markiyan Hirnyk!

I also found another solution (maybe a bit slower) but also usefull:

> stlsl := {-(40.00000000*(0.4469633000e-3*(diff(y(t), t))^2+0.4469633000e-3*(diff(x(t), t))^2))*(diff(x(t), t))/sqrt((diff(y(t), t))^2+(diff(x(t), t))^2) = diff(x(t), t, t), (40.00000000*(0.4469633000e-3*(diff(y(t), t))^2+0.4469633000e-3*(diff(x(t), t))^2))*(diff(y(t), t))/sqrt((diff(y(t), t))^2+(diff(x(t), t))^2)-9.810000000 = diff(y(t), t, t), x(0) = 0, y(0) = .80, (D(x))(0) = 8.5*cos(0.1744444444e-1*angle), (D(y))(0) = 8.5*sin(0.1744444444e-1*angle)}

>sol := dsolve(eval(stlsl,angle = angeltest), numeric):
>subs(opl(4), x(t));

Using a double for loop, I can do exactly the same ;-)

Page 1 of 1