Question: pde system solve

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

Please Wait...