mpowers

8 Reputation

2 Badges

18 years, 236 days

MaplePrimes Activity


These are replies submitted by mpowers

I'm running Maple 10. The problem is that the function "inversion(z)" can't be expressed separately... I have to cut and past the expression in literally for every reference to inversion (z) Here's what I want my code to be: dsys := { > > diff(P_p_co(z),z)=(alpha_p+g_p)*inversion(z)*P_p_co(z)+g_p*inversion(z)*m*h*nu_p*nu_bw_p-(alpha_p+bgLoss_p)*P_p_co(z), > > diff(P_p_counter(z),z)=-(alpha_p+g_p)*inversion(z)*P_p_counter(z)-g_p*inversion(z)*m*h*nu_p*nu_bw_p+(alpha_p+bgLoss_p)*P_p_counter(z), > > inversion(z) = (((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta))), > > diff(P_s(z),z)=(alpha_s+g_s)*inversion(z)*P_s(z)+g_s*inversion(z)*m*h*nu_s*nu_bw_s-(alpha_s+bgLoss_s)*P_s(z), > > P_s(0) = 0.050e-3, P_p_co(0) = 280e-3, P_p_counter(25) = 0 > > }; --------------------------------------------------------------- Here's how it actually is after cutting and pasting: dsys := { > > diff(P_p_co(z),z)=(alpha_p+g_p)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_p_co(z)+g_p*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_p*nu_bw_p-(alpha_p+bgLoss_p)*P_p_co(z), > > diff(P_p_counter(z),z)=-(alpha_p+g_p)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_p_counter(z)-g_p*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_p*nu_bw_p+(alpha_p+bgLoss_p)*P_p_counter(z), > > > diff(P_s(z),z)=(alpha_s+g_s)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_s(z)+g_s*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_s*nu_bw_s-(alpha_s+bgLoss_s)*P_s(z), > > P_s(0) = .001, P_p_counter(L) = .1, P_p_co(0) = .1};
I'm running Maple 10. The problem is that the function "inversion(z)" can't be expressed separately... I have to cut and past the expression in literally for every reference to inversion (z) Here's what I want my code to be: dsys := { > > diff(P_p_co(z),z)=(alpha_p+g_p)*inversion(z)*P_p_co(z)+g_p*inversion(z)*m*h*nu_p*nu_bw_p-(alpha_p+bgLoss_p)*P_p_co(z), > > diff(P_p_counter(z),z)=-(alpha_p+g_p)*inversion(z)*P_p_counter(z)-g_p*inversion(z)*m*h*nu_p*nu_bw_p+(alpha_p+bgLoss_p)*P_p_counter(z), > > inversion(z) = (((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta))), > > diff(P_s(z),z)=(alpha_s+g_s)*inversion(z)*P_s(z)+g_s*inversion(z)*m*h*nu_s*nu_bw_s-(alpha_s+bgLoss_s)*P_s(z), > > P_s(0) = 0.050e-3, P_p_co(0) = 280e-3, P_p_counter(25) = 0 > > }; --------------------------------------------------------------- Here's how it actually is after cutting and pasting: dsys := { > > diff(P_p_co(z),z)=(alpha_p+g_p)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_p_co(z)+g_p*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_p*nu_bw_p-(alpha_p+bgLoss_p)*P_p_co(z), > > diff(P_p_counter(z),z)=-(alpha_p+g_p)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_p_counter(z)-g_p*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_p*nu_bw_p+(alpha_p+bgLoss_p)*P_p_counter(z), > > > diff(P_s(z),z)=(alpha_s+g_s)*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*P_s(z)+g_s*(((P_s(z)*alpha_s)/(h*nu_s*zeta)) + ((P_p_co(z)*alpha_p)/(h*nu_p*zeta)) + ((P_p_counter(z)*alpha_p)/(h*nu_p*zeta)))/(1+((P_s(z)*(alpha_s+g_s))/(h*nu_s*zeta))+((P_p_co(z)*(alpha_p+g_p))/(h*nu_p*zeta))+((P_p_counter(z)*(alpha_p+g_p))/(h*nu_p*zeta)))*m*h*nu_s*nu_bw_s-(alpha_s+bgLoss_s)*P_s(z), > > P_s(0) = .001, P_p_counter(L) = .1, P_p_co(0) = .1};
Page 1 of 1