Items tagged with equation

Feed

hi guys , i have this warning for solving a complicated equation with 7 parameters. how can i overcome to this warning ?

 


odesys := {(1/4)*(-4*r^(2+p+a)*p*a-11*r^(2+a+c)*a*c-4*r^(2+p+c)*p*c+22*r^(2+a+n)*a*n+8*r^(2+p+n)*p*n-22*r^(2+a+c)*a+8*r^(2+p+a)*p-8*r^(2+p+c)*p+22*r^(2+b)*b+32*r^(2+p)*p^2+32*r^(2+p)*p+22*r^(2+b)*b^2+22*r^(2+2*a)*a+65*r^(2+2*a)*a^2-8*r^(2+p+c)*p^2+8*r^(2+p+a)*p^2-22*r^(2+a+c)*a^2)/r^4+(1/4)*(4*r^(a+n)*n^2-2*r^(n+c)*n*c-4*r^(n+c)*n^2+3*r^(2*n)*n^2-4*r^(a+c)*c+8*r^(a+n)*n+4*r^(2*c)*c-8*r^(n+c)*n+4*r^m*m^2-4*r^d*d+8*r^m*m+4*r^m-4*r^d)/r^2}

{(1/4)*(-4*r^(2+p+a)*p*a-11*r^(2+a+c)*a*c-4*r^(2+p+c)*p*c+22*r^(2+a+n)*a*n+8*r^(2+p+n)*p*n-22*r^(2+a+c)*a+8*r^(2+p+a)*p-8*r^(2+p+c)*p+22*r^(2+b)*b+32*r^(2+p)*p^2+32*r^(2+p)*p+22*r^(2+b)*b^2+22*r^(2+2*a)*a+65*r^(2+2*a)*a^2-8*r^(2+p+c)*p^2+8*r^(2+p+a)*p^2-22*r^(2+a+c)*a^2)/r^4+(1/4)*(4*r^(a+n)*n^2-2*r^(n+c)*n*c-4*r^(n+c)*n^2+3*r^(2*n)*n^2-4*r^(a+c)*c+8*r^(a+n)*n+4*r^(2*c)*c-8*r^(n+c)*n+4*r^m*m^2-4*r^d*d+8*r^m*m+4*r^m-4*r^d)/r^2}

(1)

res := op(odesys);

(1/4)*(-4*r^(2+p+a)*p*a-11*r^(2+a+c)*a*c-4*r^(2+p+c)*p*c+22*r^(2+a+n)*a*n+8*r^(2+p+n)*p*n-22*r^(2+a+c)*a+8*r^(2+p+a)*p-8*r^(2+p+c)*p+22*r^(2+b)*b+32*r^(2+p)*p^2+32*r^(2+p)*p+22*r^(2+b)*b^2+22*r^(2+2*a)*a+65*r^(2+2*a)*a^2-8*r^(2+p+c)*p^2+8*r^(2+p+a)*p^2-22*r^(2+a+c)*a^2)/r^4+(1/4)*(4*r^(a+n)*n^2-2*r^(n+c)*n*c-4*r^(n+c)*n^2+3*r^(2*n)*n^2-4*r^(a+c)*c+8*r^(a+n)*n+4*r^(2*c)*c-8*r^(n+c)*n+4*r^m*m^2-4*r^d*d+8*r^m*m+4*r^m-4*r^d)/r^2

(2)

SOL1 := solve(identity(res = 0, r), {a, b, c, d, m, n, p})

Warning, solutions may have been lost

 

``


Download sol.mw

thanks

guys, i have this equation which a, n, c, phi[0] are parameters and r is variable. maple solved this equation with n=0,c=c,phi[0]=5/4-1/4 c but i obtained another solution for this equation : a = 1, c = 5, n = 1, phi[0] = 1 ( you can check). 

 

restart

odesys := {5*r^a+4*n-c*r^n-4*phi[0]}

{5*r^a+4*n-c*r^n-4*phi[0]}

(1)

res := op(odesys);

5*r^a+4*n-c*r^n-4*phi[0]

(2)

 

SOL := solve(identity(res = 0, r), [a, n, c, phi[0]]);

[[a = 0, n = 0, c = c, phi[0] = 5/4-(1/4)*c]]

(3)

NULL

eval(5*r^a+4*n-c*r^n-4*phi[0], [a = 1, c = 5, n = 1, phi[0] = 1])

0

(4)

NULL


how can i get all solutions for a equation like this ?

Download eq000.mw

Is it possible to somehow extract a derivative from numeric solution of partial differential equation?

I know there is a command that does it for dsolve but i couldn't find the same thing for pdsolve.

The actual problem i have is that i have to take a numeric solution, calculate a derivative from it and later use it somewhere else, but the solution that i have is just a set of numbers an array of some sort and i can't really do that because obviously i will get a zero each time.

Perhaps there is a way to interpolate this numeric solution somehow?

I found that someone asked a similar question earlier but i couldn't find an answer for it.

Hello everyoene, please i have a problem solving this delay differential equation:

y'(x)=cos(x)+y(y(x)-2)    0<x<=10

y(x)=1      x<=0

tau=x-y(x)+2

please i need the solution urgently

guys ,

in a differential equation i want to expand its variable, but i  have some problem with it for exponential term :jadid.mw

 

 

thanks guys

 

 

Hi, I'm trying to solve a system of equation and I keep getting this error. Could anyone help me figure out what I'm doing wrong?

My problem is:

> alpha := .3; G := 3.5; L := 6; f := 1.1;

for i to 50 do

I0 := x(z)+y[i](z); ICon := x(0) = 1, y[i](0) = 0;

for j to 50 do

i <> j;

d1 := diff(x(z), z) = -G*x(z)*y[i](z)/IC-alpha*x(z);

d2 := diff(y[i](z), z) = G*y[i](z)*y[j](z)/IC-alpha*y[i](z);

dsys := {d1, d2};

F := dsolve({ICon, op(dsys)}, [x(z), y[i](z)], numeric);

end do;

end do;
Error, (in dsolve/numeric/process_input) unknown y[2] present in ODE system is not a specified dependent variable or evaluatable procedure


 

can anybody help me? i want to check the consistency of my scheme. My equation is too long if i check manually, so i used maple 13 to simplify my equation. But it cannot simplify it because of length of output exceed limit 1000000

restart

eqn1 := u+(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2)+(u-(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2))-u-(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)-u+(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3))*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3));

(1-exp(-m))^2*u[tt]-(1-exp(-m))^2*u[xx]+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-16*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+4*(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3))

(1)

a := simplify(eqn1);

`[Length of output exceeds limit of 1000000]`

(2)

``


Download consistency_expmle_4.mw

.

i want to find the stability of this equation, but there is seem to have some problems..can somebody help me..

 

y := A*(1/x+x*exp(-2*sqrt(-1)*b))+4*(exp(h)-1)^2*(2*exp(-sqrt(-1)*b)-3*(exp(h)-1)^2*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1))+3*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1)))/(3*(1-r))-exp(-2*sqrt(-1)*b)/x-x;

subs(A = (1+r)/(1-r), %);

subs(r = (1/3)*(exp(h)-1)^2, %);

subs(b = m*Pi*(exp(h)-1), %);

subs(m = 1, %);

subs(h = 0.5e-1, %);

> ans := solve(%, x);
Warning, solutions may have been lost
ans:=

> r1 := ans[1];
Error, invalid subscript selector
> r2 := ans[2];
Error, invalid subscript selector

Any suggestions (or perhaps related examples?) illustrating how I might numerically solve for f(t) in the following non-linear integral equation?  In Fortran, I would start with a guess f(t)=T0, and then search in the neighborhood for a minimum (in the error), but I am not familiar with numerical searches and methods in Maple.  Thank you for any suggestions or leads.

(a,b,... etc are all real)


T__0 := 298.

`&Delta;T` := 25.

0 < beta and beta <= 1

``

f*t = T[0]+`&Delta;T`*[1-exp(-a(int(exp(-b/f(y)), y = y[1] .. t))^beta)]

NULL


Download Integraleqn.mw

 

 

 

hi, the equations read as 

eq:=[2*x-0.2e-1*y-2.04*sqrt(-v^2+1)*v, 2*y-0.2e-1*x-2.16*sqrt(-u^2+1)*u, 2*u+2.16*u^2*y/sqrt(-u^2+1)-2.16*sqrt(-u^2+1)*y, 2*v+2.04*v^2*x/sqrt(-v^2+1)-2.04*sqrt(-v^2+1)*x] ;

i do as follows using DirectSearch package v.2

i find the solutions not the same,some time the results not much difference,but another,sols1 have one solution,sols2 have three solutions.in some time,some solutions are lost,the result show  me  random.may i have run the command serveral times? regards.

Can I use Maple to solve equation like |a-b|+\sqrt{2b+c}+c^2-c+1/4=0 for a,b,c?

 

a,b,c are real numbers and I need to solve it in real domain.

Hey,

I've been trying to divide an equation into matrixes, with not a lot of luck.

Simple example:

A=2b remade into [A]=[2]*[b]

This is a basic overview of the problem, i need it to work with bigger matrixes.

Thank you in advance.

Hello!
Please help solve this equation!

restart;

Dd:=4.9*10^(-10):
n:=4.5:
m:=3.5:
ny:=20:
B:=1.47*10^(-11):
H1:=70:
H2:=20:
H:=100:
b1:=40:
b2:=80:
delt_y:=H/ny:
A_max:=b1*H1+b2*H2+10*b1:
M:=1000:




b:=y->piecewise(`and`(y>=0,y<=H1),b1,`and`(y>H1,y<H1+H2),b2,`and`(y>=H1+H2,y<=H),b1):
y:=0: eq:=0:

for j from 1 to ny do
y:=(j-1)*delt_y+delt_y/2;

eq:=evalf(eq+b(y)*sign(y-y0)*abs((y-y0))^(1/n)*delt_y);
end do;


#eq=0, y0-?
sol1:=fsolve({eq=0},{y0});

Hi all,

 

It's been a while since I have used Maple. To be honest I haven't used it for over six years.

 

I am trying to solve simple differential equations, however I have many issues.

 

I am trying to simulate what author of this paper did 06421188.pdf

 

My file looks like this (Pendulum.mw)

 

Can someone help me to simulate this system? I simply can't remember how to do it.

 

Cheers,

Bart

Hello there

I'm quite an amature so please don't judge.  I'm trying to use fsolve to solve a system of non-linear equations but Maple is just "spitting" on me the equations with no intention to solve them:

> delta5 := P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+zeq^2)*(sqrt(x^2+zeq^2)*x))+x*zeq/sqrt(x^2+zeq^2)^3)/(2*Pi*E5);
print(`output redirected...`); # input placeholder
> shrinkage := P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+Zb^2)*(sqrt(x^2+Zb^2)*x))+x*Zb/sqrt(x^2+Zb^2)^3)/(2*Pi*E5)-P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+Za^2)*(sqrt(x^2+Za^2)*x))+x*Za/sqrt(x^2+Za^2)^3)/(2*Pi*E5);
> eq10 := subs(x = 1800, delta5)+subs(x = 1800, Zb = z2, Za = z1, shrinkage)+subs(x = 1800, Zb = z3, Za = z2, shrinkage)+subs(x = 1800, Zb = z4, Za = z3, shrinkage)+subs(x = 1800, Zb = z5, Za = z4, shrinkage) = 36.7*10^(-3);
print(`output redirected...`); # input placeholder
> eq9 := subs(x = 1500, delta5)+subs(x = 1500, Zb = z2, Za = z1, shrinkage)+subs(x = 1500, Zb = z3, Za = z2, shrinkage)+subs(x = 1500, Zb = z4, Za = z3, shrinkage)+subs(x = 1500, Zb = z5, Za = z4, shrinkage) = 47.2*10^(-3);
print(`output redirected...`); # input placeholder
> eq8 := subs(x = 1200, delta5)+subs(x = 1200, Zb = z2, Za = z1, shrinkage)+subs(x = 1200, Zb = z3, Za = z2, shrinkage)+subs(x = 1200, Zb = z4, Za = z3, shrinkage)+subs(x = 1200, Zb = z5, Za = z4, shrinkage) = 63.8*10^(-3);
> eq7 := subs(x = 900, delta5)+subs(x = 900, Zb = z2, Za = z1, shrinkage)+subs(x = 900, Zb = z3, Za = z2, shrinkage)+subs(x = 900, Zb = z4, Za = z3, shrinkage)+subs(x = 900, Zb = z5, Za = z4, shrinkage) = 91.1*10^(-3);
print(`output redirected...`); # input placeholder
> eq6 := subs(x = 600, delta5)+subs(x = 600, Zb = z2, Za = z1, shrinkage)+subs(x = 600, Zb = z3, Za = z2, shrinkage)+subs(x = 600, Zb = z4, Za = z3, shrinkage)+subs(x = 600, Zb = z5, Za = z4, shrinkage) = 137.9*10^(-3);
> eq5 := subs(x = 450, delta5)+subs(x = 450, Zb = z2, Za = z1, shrinkage)+subs(x = 450, Zb = z3, Za = z2, shrinkage)+subs(x = 450, Zb = z4, Za = z3, shrinkage)+subs(x = 450, Zb = z5, Za = z4, shrinkage) = 175.2*10^(-3);
> eq4 := subs(x = 300, delta5)+subs(x = 300, Zb = z2, Za = z1, shrinkage)+subs(x = 300, Zb = z3, Za = z2, shrinkage)+subs(x = 300, Zb = z4, Za = z3, shrinkage)+subs(x = 300, Zb = z5, Za = z4, shrinkage) = 230.9*10^(-3);
print(`output redirected...`); # input placeholder
> sys := {eq10, eq5, eq6, eq7, eq8, eq9};
print(`output redirected...`); # input placeholder
> fsolve(sys, {E1 = 1000 .. 2000, E2 = 0 .. 2000, E3 = 0 .. 2000, E4 = 0 .. 2000, E5 = 0 .. 2000, h4 = 100 .. 400});

and this is what Maple gives after the fsolve

 

fsolve({(3937.500000*(.2/(202500+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(450*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(202500+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.3888888889e-2/E5+(3937.500000*(.2/(202500+(650+h4)^2)+(450*(650+h4))/(202500+(650+h4)^2)^(3/2)))/E5 = .1752000000, (3937.500000*(.2/(360000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(600*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(360000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.2187500000e-2/E5+(3937.500000*(.2/(360000+(650+h4)^2)+(600*(650+h4))/(360000+(650+h4)^2)^(3/2)))/E5 = .1379000000, (3937.500000*(.2/(810000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(900*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(810000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.9722222220e-3/E5+(3937.500000*(.2/(810000+(650+h4)^2)+(900*(650+h4))/(810000+(650+h4)^2)^(3/2)))/E5 = 0.9110000000e-1, (3937.500000*(.2/(1440000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1200*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(1440000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.5468750000e-3/E5+(3937.500000*(.2/(1440000+(650+h4)^2)+(1200*(650+h4))/(1440000+(650+h4)^2)^(3/2)))/E5 = 0.6380000000e-1, (3937.500000*(.2/(2250000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1500*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(2250000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.3500000000e-3/E5+(3937.500000*(.2/(2250000+(650+h4)^2)+(1500*(650+h4))/(2250000+(650+h4)^2)^(3/2)))/E5 = 0.4720000000e-1, (3937.500000*(.2/(3240000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1800*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(3240000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.2430555555e-3/E5+(3937.500000*(.2/(3240000+(650+h4)^2)+(1800*(650+h4))/(3240000+(650+h4)^2)^(3/2)))/E5 = 0.3670000000e-1}, {E1, E2, E3, E4, E5, h4}, {E1 = 1000 .. 2000, E2 = 0 .. 2000, E3 = 0 .. 2000, E4 = 0 .. 2000, E5 = 0 .. 2000, h4 = 100 .. 400})

3 4 5 6 7 8 9 Last Page 5 of 21