Advanced Search

iman

35 Reputation

3 Badges

10 years, 331 days
Contact iman

MaplePrimes Activity


These are questions asked by iman

Parametric solution of polynomial system...

iman 35 Maple 2015
solve parameter
May 31 2016
1 11

hi.

how i can solve two equation with respect to parameter sigma1

SOLL.mw

> 

restart; pprime11 := -16395.36603*q1+5.811117425*q1*sigma1^2+3526.724044*p1-1.250000000*p1*sigma1^2+4.999870968*10^11*p1^3+4.999870970*10^11*p1*q1^2+7.967307034*10^14*p1^2*q1+4.999870966*10^12*sigma1*p1^2*q1-2.655769012*10^14*q1^3+4.999870968*10^12*sigma1*q1^3-17633.62022*q1*sigma1+6.250000000*q1*sigma1^3
> 

qprime11 := 2.655769012*10^14*p1^3-7.967307034*10^14*p1*q1^2+4.999870970*10^11*p1^2*q1-4.999870968*10^12*sigma1*p1^3-4.999870966*10^12*p1*sigma1*q1^2+3526.724044*q1-1.250000000*q1*sigma1^2+16395.36603*p1-5.811117425*p1*sigma1^2+4.999870968*10^11*q1^3+17633.62022*p1*sigma1-6.250000000*p1*sigma1^3:
> 

-50 < sigma1 and sigma1 < 50:
> 

sigma1 <> 53.11665685, -53.11665685:
> 

SOLL := solve({pprime11, qprime11}, real)

Warning,  computation interrupted

  
> 

``

 

Download SOLL.mw

How to differentiate with respect to y?...

iman 35 Maple 2015
differentiation
May 27 2016
1 0

how i can differential with respect to ''y'' variable

thanks...help_(2).mw

T(r, theta) := .5*alpha+sum([gamma*n*(r^n+R^(2*n)+r^(-n))*sin(n*theta)], n = 1 .. infinity)

.5*alpha+[-gamma*r*(r^2-1)*sin(theta)/(4*cos(theta)^2*r^2-4*cos(theta)*r^3+r^4-4*cos(theta)*r+2*r^2+1)-((1/2)*I)*R^2*gamma*(exp(I*theta)*R^4-exp((3*I)*theta)*R^4-exp(I*theta)+exp((3*I)*theta))/((R^2*exp(I*theta)-1)^2*(-R^2+exp(I*theta))^2)+((1/2)*I)*r*gamma*(exp(I*theta)*r^2-exp((3*I)*theta)*r^2-exp(I*theta)+exp((3*I)*theta))/((exp(I*theta)-r)^2*(exp(I*theta)*r-1)^2)]

 (1)

``

y := r, proc (theta) options operator, arrow; r*sin(theta) end proc

r, proc (theta) options operator, arrow; r*sin(theta) end proc

 (2)

 

 

diff(T, y)

Error, invalid input: diff received proc (theta) options operator, arrow; r*sin(theta) end proc, which is not valid for its 3rd argument

  

``

 

Download help_(2).mw

Replica of MATLAB code in Maple?...

iman 35 Maple 2015
May 22 2016
0 0

hi

if possible to convert matlab file in to maple fie program??convert_to_maple_program.txt

thanks

 

Simplification of RootOf...

iman 35 Maple 2015
solve rootof
May 22 2016
2 3

hi...how i can gain result for solve three equations,in which term '' Root of'' dont appear?

thanks

root.mw

restart; Q1 := aa*(y-x)

aa*(y-x)

 (1)

Q2 := -ll*x*z+bb*x

-ll*x*z+bb*x

 (2)

 

-ll*x*z+bb*x

 (3)

Q3 := -cc*z+hh*x*x+kk*y*y

hh*x^2+kk*y^2-cc*z

 (4)

SOLL := solve({Q1, Q2, Q3}, {x, y, z})

{x = 0, y = 0, z = 0}, {x = RootOf((hh*ll+kk*ll)*_Z^2-bb*cc), y = RootOf((hh*ll+kk*ll)*_Z^2-bb*cc), z = bb/ll}

 (5)

``

 

Download root.mw

Incorrect syntax...

iman 35 Maple 2015
syntax programming
May 20 2016
0 1

hi.please help me for remove error'' 

Error, illegal use of an object as a name''

 

thanks

PLATE.mw

   

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

Error, illegal use of an object as a name

"restart:Digits :=15: beta:=10:alpha:=100: xi:=.5: upsilon:=0.2841945289:n:=3: aa:=1:b:=1:N_x:=0.4:N_y:=0.4:N_xy:=0: hl2:=1:mu:=65.8e9:E:=169e9: delta0:=1:delta1:=1: mus:=3:D1:=2;h:=1: lambda:=0.1: D2:=5:A1:=-2:A2:=-2:A3:=-6:A4:=7:A5:=7:A6:=7:A7:=7:A8:=8:A9:=7:A10:=7:A11:=1: A12:=1:tau:=4.730040745:t:=0: g2:=sin(theta):g3:=cos(theta):g1:=cos(theta):a:=0.0:with(Student[Calculus1]): a1:=evalf((A1*ApproximateInt(g3^2,theta=a..1,method=simpson) ) ) : a2:= evalf(A2*ApproximateInt(g3*((&DifferentialD;)^2)/((&DifferentialD;theta)^2)g3,theta=a..1,method=simpson)): a3:=evalf(A3*ApproximateInt(g3*g3,theta=a..1,method=simpson)) : a4:=evalf(A4*ApproximateInt(g3*g3,theta=a..1,method=simpson)) :a5:=evalf(A5*ApproximateInt(g3^2,theta=a..1,method=simpson)) : a6:=evalf(A6*ApproximateInt(((&DifferentialD;)^2)/((&DifferentialD;theta)^2)g3*g3,theta=a..1,method=simpson)) :a7:=evalf(A7*ApproximateInt(g3*g3,theta=a..1,method=simpson)): a8:=evalf(A8*ApproximateInt(g3^2,theta=a..1,method=simpson)):a9:=evalf(ApproximateInt(A9*(g3*g3 ) ,theta=a..1,method=simpson)) :a10:=evalf(A10*ApproximateInt(g3*g3,theta=a..1,method=simpson)):a11:=evalf(ApproximateInt(1,theta=a..1,method=simpson)):a12:=evalf(ApproximateInt(1*(1-1/2 (1)),theta=a..1,method=simpson)):dsys3 := { f3(x)*(a1)+ f3(x)*(a2) +((&DifferentialD;)^2)/((&DifferentialD;x)^2) f3(x)*(a3) +f3(x)*a4+ f3(x)*(a5) +((&DifferentialD;)^2)/((&DifferentialD;x)^2) f3(x)*(a6) +f3(x)*a7= ((&DifferentialD;)^2)/((&DifferentialD;x)^2) f3(x) *(a8) + a9*(&DifferentialD;)/(&DifferentialD;x) f3(x) +f3(x)*a10+ a11+a12 , f3(1) =0,f3(0) =0 , D^(1)(f3)(1) = 0, D^(1)(f3)(0)=0,D^(3)(f3)(1) = 0, D^(3)(f3)(0)=0} :dsol5 := dsolve(dsys3, 'maxmesh'=2024, numeric,abserr=.0001, range=0..1, output=listprocedure):fy3:= eval(f3(x),dsol5):with(CurveFitting):fy33:=PolynomialInterpolation([[0,fy3(0)],[.1,fy3(0.1)],[.2,fy3(0.2)],[0.3,fy3(0.3)],[.4,fy3(0.4)],[.5,fy3(0.5)],[0.6,fy3(0.6)],[0.7,fy3(0.7)],[0.8,fy3(0.8)],[0.9,fy3(0.9)],[1,fy3(1)]],x): d1:=A1*ApproximateInt(((&DifferentialD;)^6)/((&DifferentialD;x)^6)fy33*fy33,x=a..1,method=simpson) :d2:= A2*ApproximateInt(((&DifferentialD;)^4)/((&DifferentialD;x)^4)fy33*fy33 ,x=a..1,method=simpson) :d3:=A3*ApproximateInt(((&DifferentialD;)^2)/((&DifferentialD;x)^2)fy33*fy33,x=a..1,method=simpson): d4:= A4*ApproximateInt(fy33*fy33,x=a..1,method=simpson):d5:=A5*ApproximateInt(((&DifferentialD;)^4)/((&DifferentialD;x)^4)fy33*fy33,x=a..1,method=simpson) : d6:=A6*ApproximateInt(((&DifferentialD;)^2)/((&DifferentialD;x)^2)fy33*fy33,x=a..1,method=simpson) :d7:=A7*ApproximateInt(fy33*fy33,x=a..1,method=simpson) :d8:=A8*ApproximateInt(((&DifferentialD;)^2)/((&DifferentialD;x)^2)fy33*fy33,x=a..1,method=simpson) :d9:=ApproximateInt(A9*(((&DifferentialD;)^1)/((&DifferentialD;x)^1)fy33*fy33 ) ,x=a..1,method=simpson) :d10:=A10*ApproximateInt(fy33*fy33,x=a..1,method=simpson) :d11:=evalf(ApproximateInt(1,theta=a..1,method=simpson)):d12:=evalf(ApproximateInt(1*(1-1/2 (1)),theta=a..1,method=simpson)) : d sys4 := { h3(theta)*(d1)+((&DifferentialD;)^2)/((&DifferentialD;theta)^2) h3(theta)*(d2)+((&DifferentialD;)^4)/((&DifferentialD;theta)^4) h3(theta)*(d3)+ ((&DifferentialD;)^6)/((&DifferentialD;theta)^6) h3(theta)*(d4)+h3(theta) *(d5)+ h3(theta) *(d6) +((&DifferentialD;)^4)/((&DifferentialD;theta)^4) h3(theta)*(d7)= h3(theta)*(d8) +d9*(&DifferentialD;)/(&DifferentialD;theta) h3(theta) +((&DifferentialD;)^2)/((&DifferentialD;theta)^2) h3(theta)*(d10) +d11+d12 ,h3(1) = 0,h3(0) = 0 , D^(1)(h3)(1) = 0, D^(1)(h3)(0)=0,D^(3)(h3)(1) = 0, D^(3)(h3)(0)=0} :dsol6 := dsolve(dsys4, 'maxmesh'=2024, abserr=.0001, range=0..1, numeric, output=listprocedure):g33:= eval(h3(theta),dsol6):with(CurveFitting):g3:=PolynomialInterpolation([[0,g33(0)],[.1,g33(0.1)],[.2,g33(0.2)],[0.3,g33(0.3)],[.4,g33(0.4)],[.5,g33(0.5)],[0.6,g33(0.6)],[0.7,g33(0.7)],[0.8,g33(0.8)],[0.9,g33(0.9)],[1,g33(1)]],theta):"

  

 

``

 

Download PLATE.mw

1 2 3 4 5 6 Page 2 of 6