wo0olf

20 Reputation

4 Badges

11 years, 15 days

MaplePrimes Activity


These are questions asked by wo0olf

hi,

i cant solve this 

int(1/1-exp(a*x)*erfc(a*x))

hi ,

 

i solve my odes with ode analyzer and obtained h(t) , u(t) and s(t) ... i wana put u(t) and s(t) in another equation to obtain T(t,r) , i did as follow but i have an error :

 

400*(s(t)*(diff(T(t, r), r))+u(t)*(diff(T(t, r), t))/r) = 0.3e-1*((diff(T(t, r), r)+r*(diff(T(t, t), r, r)))/r+(diff(T(t, r), r))/r^2)+0.1e-2*(2*((diff(u(t), t))/r+s(t)/r)^2+((diff(s(t), t))/r-u(t)/r)^2)

 

4.00000000000000*u(t)^2*h(t)*(diff(y(t), t))-8.00000000000000*h(t)^2*x(t)^2+16.0000000000000*u(t)*h(t)*x(t)*(diff(h(t), t))+(1367.18750000000*(-29.8200000000000+u(t)))*h(t)^2*x(t)+363.281250000000*u(t)*h(t)*y(t)-1367.18750000000*h(t)^2*(30-u(t))*x(t)+u(t)^2*h(t)-0.225000000000000e-1+0.225000000000000e-1*(0.749297072727423e-1+1.49578322716831*sin(t-0.75e-1))^2+4.00000000000000*(u(t)*y(t)+h(t)*x(t))^2 = 0

 

diff(s(t), t) = 8.00000000000000*(diff(x(t), t))*h(t)+16.0000000000000*y(t)*x(t)+8.00000000000000*u(t)*(diff(y(t), t))

 

diff(h(t), t) = y(t)

 

diff(u(t), t) = x(t)

 

x(t) = 427.2460938*u(t)+385620.1174/u(t)-25671.38673

 

h(0) = 0.1e-10, s(0) = .6, u(0) = .62, x(0) = 596561.43714, y(0) = .12096, T(0, r) = 0, t = 0 .. .18, r = 0 .. .125

 

pdsolve(1, 2, 3, 4, 5, 6, {7}, numeric);
%;
Error, (in pdsolve) invalid input: `pdsolve/numeric` expects its 2nd argument, IBCs, to be of type {list, set}, but received 2

hi,

i obtained values for h(t) and u(t) and i wana obtaine s(t) , i don't konw how can i do it

 

for h(t) and u(t) ,

 

sol1 := dsolve([4.00000000000000*u(t)^2*h(t)*(diff(diff(h(t), t), t))-8.00000000000000*h(t)^2*(427.2460938*u(t)+385620.1174/u(t)-25671.38673)^2+16.0000000000000*u(t)*h(t)*(427.2460938*u(t)+385620.1174/u(t)-25671.38673)*(diff(h(t), t))+(16406.2500000000*(-29.8200000000000+u(t)))*h(t)^2*(427.2460938*u(t)+385620.1174/u(t)-25671.38673)+7171.87500000000*u(t)*h(t)*(diff(h(t), t))-16406.2500000000*h(t)^2*(30-u(t))*(427.2460938*u(t)+385620.1174/u(t)-25671.38673)+u(t)^2*h(t)-0.225000000000000e-1+0.225000000000000e-1*(.997494986604054+.106105802501554*sin(t-1.5))^2+4.00000000000000*(h(t)*(427.2460938*u(t)+385620.1174/u(t)-25671.38673)+u(t)*(diff(h(t), t)))^2 = 0, diff(u(t), t) = 427.2460938*u(t)+385620.1174/u(t)-25671.38673, h(1.44) = 0.145e-4, (D(h))(1.44) = .12093, u(1.44) = .62], numeric)

 

and s(t) is , s(t)=1/r*d/dt(u(t)*h(t))

 

thanks.

hi,

I want to solve this to data table , i can plot but i need data in table 

 

sol1 := dsolve([diff(u(theta),theta) = 427.2461*u(theta)+385620.123/u(theta)-25671.3871, u(0) = .6], numeric);

plots[odeplot](sol1, 0..0.18, color = red);

 

thanks.

1 2 Page 2 of 2