Hi, I have a ‘slight’ problem (you will probably recognize it Joe! :-) ). It concerns the values of Tau, omega in my worksheet (see below). If I set Tau=0.7, omega=0.7*m*Eta everything is rosy, and works fine. If I start tweaking these values (which I have to) things go a bit pair shaped. I either get an error message after the first call to dsolve (e.g. when tau=0.5, omega=0.7*m*Eta) : "Error, (in dsolve/numeric/checksing) ode system has a removable singularity at r=1. Initial data is restricted to {Phi(r) = .20650095602297*diff(Phi(r),r)+.82088920025557e-1*I*diff(Phi(r),r)}"

Hi, I am trying to solve a 2nd order ODE with linear symmetries (according to odeadvisor!) subject to two boundary conditions. The ODE itself is not particularly pretty (It can't be solved analytically), and it involves quite a lot of other functions that I have to define before hand - sorry if this looks ugly: N:=2: Xi:=1: mu:=0: m:=4: l:=2: sigma:=1: Tau:=0.7: M:=1/(2*(1-(Eta)^2)): a:=Eta: Eta:=Tau*sqrt(2/3): omega:=0.9*m*Eta: g:=r->sqrt((1-2*M/r^(2*N)+2*M*a^2/r^(2*N+2))^(-1)): h:=r->sqrt(r^2*(1+2*M*a^2/r^(2*N+2))): Omega:=r->2*M*a/(r^(2*N+2)+2*M*a^2): f:=r->r/(g(r)*h(r)):

Hi, I'm trying to get maple to output the numerical answer, from a call to dsolve, to an external txt file; but I'm having a bit of trouble. I want to do this so that I may then plot the results in an external program like easy-plot, or just analyse them. I'm using a "read" command in Maple to read in a procedure from a txt file (shown below), which I then execute in Maple: #Utilise prime notation for derivatives wrt t: with(PDEtools, declare): with(plots): declare(y(t), prime=t, quiet): #Also shorten output notation: alias(y=y(t),y0=y(0),yp0=D(y)(0)): #define the ODE: ODE:=diff(y,[t$2])-(1-y^2)*diff(y,t)+y=0;