Items tagged with syntax

Feed

yep the errors recieved using some packages are very very specific for maple, for example, the one i got today using the ODE package was profoundly helpful:

 

Error, (in ODEtools/info) unable to handle derivatives as {diff(1/(ln(X)-Psi(1-f(X))-Psi(f(X))), [`$`(X, n-k[1]-k[2])]), diff(1/f(X), [`$`(X, n-k[1]-k[2])]), diff(Psi(1-f(X)), [`$`(X, k[2])]), diff(Psi(f(X)), [`$`(X, k[2])]), diff(f(X), [`$`(X, k[1])])} while solving w.r.t f(X)

Hi,

1st post. I'm trying to integrate the following function:

h:=t->(2*t-1)*cos*sqrt(3*(2*t-1)^2+6)/(sqrt(3*(2*t-1)^2+6));

h:=t->(2*t-1)*cos*(sqrt(3*(2*t-1)^2+6))/(sqrt(3*(2*t-1)^2+6));

h:=t->(2*t-1)*(cos*(sqrt(3*(2*t-1)^2+6)))/(sqrt(3*(2*t-1)^2+6));

h:=t->((2*t-1)*(cos*(sqrt(3*(2*t-1)^2+6))))/(sqrt(3*(2*t-1)^2+6));

int(h(t),t); ** Integration command. I've also replaced the "h(t)" with the entire function.

I've tried the following:

1. Changed the "t" to "x" throughout function.

2. Added parens around sqrt portion.

3. Added parens to include "cos" and then added to include the beginning (2*t-1).

4. I've added brackets around the numerator but this just causes Maple to reprint the function with the inegration sign in front of the function.

5. I've also tried using the Integration tutor. It returns that maple is unable to calculate.

6. Repeat all the above in Maple 2015, same answer.

I always get cos(t^2-t).

The math book claims the answer is 1/6*sin*sqrt(3*(2*t-1)^2+6). When I perform the inegration on paper I get the same answer.

Any suggestions or corrections would be great.

Thank you,

Jay.

 

pls help review this code, its doesnt return a solution

 

 

restart;
Digits := 16;
M := .5; lambda := .5; Pr := .72; beta := 1; L[w] := 0; m := 1; R := 1; Ec := 1;
N := 7;
for j from 0 to N do J[j] := sum(f[k](t)*(diff(f[j-k](t), `$`(t, 2))), k = 0 .. j) end do;
for i from 0 to N do K[i] := sum((diff(f[k](t), t))*(diff(f[i-k](t), t)), k = 0 .. i) end do;
for j from 0 to N do G[j] := sum(f[k](t)*(diff(theta[j-k](t), t)), k = 0 .. j) end do;
for j from 0 to N do H[j] := sum((diff(f[k](t), t))*theta[j-k](t), k = 0 .. j) end do;
for i from 0 to N do P[i] := sum((diff(f[k](t), t, t))*(diff(f[i-k](t), t)), k = 0 .. i) end do;
epsilon := 1; delta := 0;
f[0] := proc (t) options operator, arrow; L[w]+epsilon+delta*A*t+(1/2)*A*t^2 end proc;
1 2
t -> L[w] + epsilon + delta A t + - A t
2
theta[0] := proc (t) options operator, arrow; 1+B*t end proc;
t -> 1 + B t
NULL;
;
NULL;
NULL;
NULL;
NULL;
for i to N do f[i] := simplify(-((m+1)*(1/2))*(int(int(int(J[i-1], t = 0 .. eta), t = 0 .. eta), t = 0 .. eta))+m*(int(int(int(1-K[i-1], t = 0 .. eta), t = 0 .. eta), t = 0 .. eta))-M*(int(int(int(diff(f[i-1](t), t)-1, t = 0 .. eta), t = 0 .. eta), t = 0 .. eta))-lambda*(int(int(int(theta[i-1](t), t = 0 .. eta), t = 0 .. eta), t = 0 .. eta))); f[i] := unapply(f[i], eta); theta[i] := simplify(-3*Pr*R*(((m+1)*(1/2))*(int(int(G[i-1], t = 0 .. eta), t = 0 .. eta))-(2*m-1)*(int(int(H[i-1], t = 0 .. eta), t = 0 .. eta))+Ec*(int(int(P[i-1], t = 0 .. eta), t = 0 .. eta)))/(4+3*R)); theta[i] := unapply(theta[i], eta) end do;
NULL;
F(eta):=collect((∑)f[z](eta),eta):
Theta(eta):=collect((∑)theta[z](eta),eta):
with(numapprox);
for k from 2 to 5 do W[k] := pade(diff(F(eta), eta), eta, [k, k]); Q[k] := pade(Theta(eta), eta, [k, k]); SOLL1[k] := expand(coeff(numer(W[k]), eta^k)) = 1; SOLL2[k] := expand(coeff(numer(Q[k]), eta^k)) = 0; SOL[k] := solve({SOLL1[k], SOLL2[k]}, {A, B}); print([k] = SOL[k]) end do;
Warning, computation interrupted

 

 

 

 

 

 

I'm sorry for asking too many similar technical questions, but I just can't help this.

I'm solving a system of differential equations with dsolve and getting this type of error:

Error, (in f) unable to store 'HFloat(0.10664489706950975)+HFloat(1.1891638418458722e-5)*sin^2-HFloat(1.6095871822513048e-6)*sin' when datatype=float[8]

I checked the dsolve syntax, checked that all the constants are defined, checked that the number of equations matches the number of unknown functions (and there is no misprints of them in the code). In brief, I checked all the stuff that I usually have mistakes in, and yet the error stands. 

Here's the file: 1.mw

cannot find the error in loop

 

 

> restart; u[0] := (4/3)*c^2*cos((1/4)*x)^2; alpha := 2;
 
> iteration := 3;
> for k from 0 while k <= Iteration do u[s] := eval(u[k], t = xi); u[k+1] := simplify(u[k]-(int(diff(u[s], [`$`(xi, alpha)])+diff(u[s]*u[s], x)+diff(u[s]*u[s], x, x, x), xi = 0 .. t))) end do;

how can solve using assume option ?

restart:

F(z):=m*z^4-4*m*z^3+(3*m+3)*z^2-6*u;

m*z^4-4*m*z^3+(3*m+3)*z^2-6*u

(1)

solve(F(z),z assuming -0.5<m<0.5, 0<u<0.5)

Error, `<` unexpected

 

 

 

Download assume.mw

Hi,

When I type sin(pi/2) the result is sin(pi/2) (not "1.0)

What should I do to get "1.0" instead of replicating the sin(pi/2) ?!

I mean why the expression contating "pi" can not being simplified?

Is there any solution to this problem?

 

Thanks

how are iterated functions represented in maple? as in f(f(f(x))) is f^3(x)(x)  in conventional notation where by the reader knows it is refering to the iteration conducted 3 times on the argument x, but what does maple use to differentiate between iteration ,exponentiation and differentiation?

when writing by error the following

restart;
f:=proc()
dsolve(diff(y(x),x)+y(x)=0,y(x));
end proc();

Maple did not complain and returned  f := y(x) = _C1*exp(-x)

Then I noticed I needed to change "end proc();" to "end proc;" 

The question is, why did Maple not generate a syntax error? If "end proc();" is valid Maple code, then what does it mean? or Did Maple ignore it? what is the parsing steps used to make Maple generate the above output? 

Hello every one,please help me for plotting this equation(T).L=1,t=10,a=0.1,x=0..1

I have two sets 

f:={1,2,3,4};

h:={1,2,4,5}

L=seq(i,i=1..4):

I want to program , if it is true that f[i] = h[i], then it prints f[i]. So the output should be {1,2}.

for i in L do
if evalb(f[i]=h[i]) then
print(f[i]);
end if;
end do;

I get no output. https://i.imgur.com/qA5hU3i.png

I tried changing the set f to list, f:= [1,2,3,4], still no output.

hi .please help me for solve this equations.

bbb2.mw

restart; d[11] := 1; mu[11] := 1; q[311] := 1; d[33] := 1; mu[33] := 1; a[11] := 1; e[311] := 1; a[33] := 1; A := 1; g[111111] := 1; c[1111] := 1; g[113113] := 1; f[3113] := 1; beta[11] := 1; `&Delta;T` := 1; II := 1; L := 1

J := d[11]*(diff(Phi(x, z), x, x))+mu[11]*(diff(psi(x, z), x, x))+q[311]*(diff(w(x), x, x))+d[33]*(diff(Phi(x, z), z, z))+mu[33]*(diff(psi(x, z), z, z));

diff(diff(Phi(x, z), x), x)+diff(diff(psi(x, z), x), x)+diff(diff(w(x), x), x)+diff(diff(Phi(x, z), z), z)+diff(diff(psi(x, z), z), z)

(1)

B := a[11]*(diff(Phi(x, z), x, x))+d[11]*(diff(psi(x, z), x, x))+e[311]*(diff(w(x), x, x))+a[33]*(diff(Phi(x, z), z, z))+d[33]*(diff(psi(x, z), z, z));

diff(diff(Phi(x, z), x), x)+diff(diff(psi(x, z), x), x)+diff(diff(w(x), x), x)+diff(diff(Phi(x, z), z), z)+diff(diff(psi(x, z), z), z)

(2)

R := A*(g[111111]*(diff(u[0](x), x, x, x, x))-c[1111]*(diff(u[0](x), x, x)+(1/2)*(diff((diff(w(x), x))^2, x)))+e[311]*(diff(diff(Phi(x, z), z), x))+q[311]*(diff(diff(psi(x, z), z), x)));

diff(diff(diff(diff(u[0](x), x), x), x), x)-(diff(diff(u[0](x), x), x))-(diff(w(x), x))*(diff(diff(w(x), x), x))+diff(diff(Phi(x, z), x), z)+diff(diff(psi(x, z), x), z)

(3)

S := -II*g[111111]*(diff(w(x), x, x, x, x, x, x))-II*c[1111]*(diff(w(x), x, x, x, x))+A*g[113113]*(diff(w(x), x, x, x, x))-A*f[3113]*(diff(diff(Phi(x, z), z), x, x))-A*(c[1111]*(diff(u[0](x), x, x)+(1/2)*(diff((diff(w(x), x))^2, x)))+e[311]*(diff(diff(Phi(x, z), z), x))+q[311]*(diff(diff(psi(x, z), z), x)))*(diff(w(x), x))-A*(diff(w(x), x, x))*(c[1111]*(diff(u[0](x), x)+(1/2)*(diff(w(x), x))^2)+e[311]*(diff(Phi(x, z), z))+q[311]*(diff(psi(x, z), z))-beta[11]*`&Delta;T`);

-(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))-(diff(diff(diff(Phi(x, z), x), x), z))-(diff(diff(u[0](x), x), x)+(diff(w(x), x))*(diff(diff(w(x), x), x))+diff(diff(Phi(x, z), x), z)+diff(diff(psi(x, z), x), z))*(diff(w(x), x))-(diff(diff(w(x), x), x))*(diff(u[0](x), x)+(1/2)*(diff(w(x), x))^2+diff(Phi(x, z), z)+diff(psi(x, z), z)-1)

(4)

dsys := {B, J, R, S}; BCS := {D@@2*w(0) = 0, D@@2*w(L) = 0, Phi(x = 0) = 0, Phi(x = L) = 0, Phi(z = -(1/2)*h) = 0, Phi(z = (1/2)*h) = 0, psi(x = 0) = 0, psi(x = L) = 0, psi(z = -(1/2)*h) = 0, psi(z = (1/2)*h) = 0, w(x = 0) = 0, w(x = L) = 0, u[0](x = 0) = 0, u[0](x = L) = 0, (D(w))(0) = 0, (D(w))(L) = 0, (D(u[0]))(0) = 0, (D(u[0]))(L) = 0}

{D@@2*w(0) = 0, D@@2*w(L) = 0, Phi(x = 0) = 0, Phi(x = L) = 0, Phi(z = -(1/2)*h) = 0, Phi(z = (1/2)*h) = 0, psi(x = 0) = 0, psi(x = L) = 0, psi(z = -(1/2)*h) = 0, psi(z = (1/2)*h) = 0, w(x = 0) = 0, w(x = L) = 0, u[0](x = 0) = 0, u[0](x = L) = 0, (D(w))(0) = 0, (D(w))(L) = 0, (D(u[0]))(0) = 0, (D(u[0]))(L) = 0}

(5)

dsol5 := dsolve(dsys, numeric)

Error, (in dsolve/numeric/process_input) missing differential equations and initial or boundary conditions in the first argument: dsys

 

NULL

NULL

NULL

if former equations are not solvable , please help me for another way, in which at first two equation solve..in this way in equation [J and B] assume that q[311]=e[311]=0 and dsolve perform to find Φ and  ψ

after by finding Φ and  ψ is use for detemine w and u0

please see attached file below[bbb2_2.mw]

bbb2_2.mw

Download bbb2.mw

 

hi all,

How could I write long commands in Maple?

For instance: a long vector in Matlab could be written in this way:

A=[ 3, 4, 5, 6, 6, 45, 37...

5, 4, 67, 39, -967 ];

But what is that in Maple??? 


Trying to build a block matrix. Having a problem getting the syntax correct. Can't add in the predefind matrices to the lower band.

restart

with(LinearAlgebra):

interface(displayprecision = 5)

5

(1)

interface(rtablesize = 81)

10

(2)

S := 2

2

(3)

dmax := 7

7

(4)

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

(5)

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

(6)

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

(7)

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

(8)

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

(9)

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

(10)

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

(11)

S*dmax

14

(12)

M1 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), Matrix(S, S*dmax, [seq(-CCnew || n, n = 0 .. dmax-1)])])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -336750255587/3769550688757, (1, 2) = 14853552191797/1696297809940650, (1, 3) = -49655436033349/56543260331355, (1, 4) = 75647656451147/1413581508283875, (1, 5) = -299962512141959/80776086187650, (1, 6) = 1231816081155781/8481489049703250, (1, 7) = -50445725001719/5769720441975, (1, 8) = 9065291388901/40388043093825, (1, 9) = -142685068141037/11539440883950, (1, 10) = 116560067351321/565432603313550, (1, 11) = -60560690824604/5769720441975, (1, 12) = 88774498025543/848148904970325, (1, 13) = -20398649489879/4121228887125, (1, 14) = 5446073753219/242328258562950, ...

 

M2 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity)])

M2 := Matrix(14, 14, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 1, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (1, 10) = 0, (1, 11) = 0, (1, 12) = 0, (1, 13) = 0, (1, 14) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 1, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (2, 10) = 0, (2, 11) = 0, (2, 12) = 0, (2, 13) = 0, (2, 14) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (3, 5) = 1, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (3, 10) = 0, (3, 11) = 0, (3, 12) = 0, (3, 13) = 0, (3, 14) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 0, (4, 5) = 0, (4, 6) = 1, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (4, 10) = 0, (4, 11) = 0, (4, 12) = 0, (4, 13) = 0, (4, 14) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (5, 7) = 1, (5, 8) = 0, (5, 9) = 0, (5, 10) = 0, (5, 11) = 0, (5, 12) = 0, (5, 13) = 0, (5, 14) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0, (6, 7) = 0, (6, 8) = 1, (6, 9) = 0, (6, 10) = 0, (6, 11) = 0, (6, 12) = 0, (6, 13) = 0, (6, 14) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 0, (7, 5) = 0, (7, 6) = 0, (7, 7) = 0, (7, 8) = 0, (7, 9) = 1, (7, 10) = 0, (7, 11) = 0, (7, 12) = 0, (7, 13) = 0, (7, 14) = 0, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 0, (8, 5) = 0, (8, 6) = 0, (8, 7) = 0, (8, 8) = 0, (8, 9) = 0, (8, 10) = 1, (8, 11) = 0, (8, 12) = 0, (8, 13) = 0, (8, 14) = 0, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = 0, (9, 6) = 0, (9, 7) = 0, (9, 8) = 0, (9, 9) = 0, (9, 10) = 0, (9, 11) = 1, (9, 12) = 0, (9, 13) = 0, (9, 14) = 0, (10, 1) = 0, (10, 2) = 0, (10, 3) = 0, (10, 4) = 0, (10, 5) = 0, (10, 6) = 0, (10, 7) = 0, (10, 8) = 0, (10, 9) = 0, (10, 10) = 0, (10, 11) = 0, (10, 12) = 1, (10, 13) = 0, (10, 14) = 0, (11, 1) = 0, (11, 2) = 0, (11, 3) = 0, (11, 4) = 0, (11, 5) = 0, (11, 6) = 0, (11, 7) = 0, (11, 8) = 0, (11, 9) = 0, (11, 10) = 0, (11, 11) = 0, (11, 12) = 0, (11, 13) = 1, (11, 14) = 0, (12, 1) = 0, (12, 2) = 0, (12, 3) = 0, (12, 4) = 0, (12, 5) = 0, (12, 6) = 0, (12, 7) = 0, (12, 8) = 0, (12, 9) = 0, (12, 10) = 0, (12, 11) = 0, (12, 12) = 0, (12, 13) = 0, (12, 14) = 1, (13, 1) = 0, (13, 2) = 0, (13, 3) = 0, (13, 4) = 0, (13, 5) = 0, (13, 6) = 0, (13, 7) = 0, (13, 8) = 0, (13, 9) = 0, (13, 10) = 0, (13, 11) = 0, (13, 12) = 0, (13, 13) = 0, (13, 14) = 0, (14, 1) = 0, (14, 2) = 0, (14, 3) = 0, (14, 4) = 0, (14, 5) = 0, (14, 6) = 0, (14, 7) = 0, (14, 8) = 0, (14, 9) = 0, (14, 10) = 0, (14, 11) = 0, (14, 12) = 0, (14, 13) = 0, (14, 14) = 0})

(13)

lowerband := Matrix(S, S*dmax, [seq(-evalf(CCnew || n), n = 0 .. dmax-1)])

lowerband := Matrix(2, 14, {(1, 1) = -0.893343e-1, (1, 2) = 0.875645e-2, (1, 3) = -.878185, (1, 4) = 0.535149e-1, (1, 5) = -3.71351, (1, 6) = .145236, (1, 7) = -8.74318, (1, 8) = .224455, (1, 9) = -12.36499, (1, 10) = .206143, (1, 11) = -10.49630, (1, 12) = .104669, (1, 13) = -4.94965, (1, 14) = 0.224740e-1, (2, 1) = -0.864552e-2, (2, 2) = -.115162, (2, 3) = -0.776011e-1, (2, 4) = -1.10141, (2, 5) = -.282396, (2, 6) = -4.51408, (2, 7) = -.537873, (2, 8) = -10.26209, (2, 9) = -.568850, (2, 10) = -13.97156, (2, 11) = -.318088, (2, 12) = -11.39493, (2, 13) = -0.737097e-1, (2, 14) = -5.15750})

(14)

M3 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), lowerband])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -0.8933432215e-1, (1, 2) = 0.8756453086e-2, (1, 3) = -.8781848755, (1, 4) = 0.5351488825e-1, (1, 5) = -3.713506389, (1, 6) = .1452358276, (1, 7) = -8.743183575, (1, 8) = .224454831, (1, 9) = -12.36498974, (1, 10) = .2061431666, (1, 11) = -10.49629552, (1, 12) = .1046685287, (1, 13) = -4.949652167, (1, 14) = 0.2247395242e-1, (2, 1) = -0.8645515381e-2, (2, 2) = -.1151624586, (2, 3) = -0.7760109089e-1, (2, 4) = -1.101411529, (2, 5) = -.2823958356, (2, 6) = -4.514076892, (2, 7) = -.537872817, (2, 8) = -10.26209174, (2, 9) = -.5688498906, (2, 10) = -13.97155838, (2, 1...

 

``


Download Test_Block_matrix.mw

restart

with(LinearAlgebra):

interface(displayprecision = 5)

5

(1)

interface(rtablesize = 81)

10

(2)

S := 2

2

(3)

dmax := 7

7

(4)

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

CCnew0 := Matrix(2, 2, {(1, 1) = 336750255587/3769550688757, (1, 2) = -14853552191797/1696297809940650, (2, 1) = 665096091/76929605893, (2, 2) = 1328910382993/11539440883950})

(5)

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

CCnew1 := Matrix(2, 2, {(1, 1) = 49655436033349/56543260331355, (1, 2) = -75647656451147/1413581508283875, (2, 1) = 29849106694/384648029465, (2, 2) = 10591394356218/9616200736625})

(6)

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

CCnew2 := Matrix(2, 2, {(1, 1) = 299962512141959/80776086187650, (1, 2) = -1231816081155781/8481489049703250, (2, 1) = 155175716729/549497184950, (2, 2) = 260449617208489/57697204419750})

(7)

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

CCnew3 := Matrix(2, 2, {(1, 1) = 50445725001719/5769720441975, (1, 2) = -9065291388901/40388043093825, (2, 1) = 21111399914/39249798925, (2, 2) = 2819495262394/274748592475})

(8)

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

CCnew4 := Matrix(2, 2, {(1, 1) = 142685068141037/11539440883950, (1, 2) = -116560067351321/565432603313550, (2, 1) = 44654487647/78499597850, (2, 2) = 53741323977599/3846480294650})

(9)

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

CCnew5 := Matrix(2, 2, {(1, 1) = 60560690824604/5769720441975, (1, 2) = -88774498025543/848148904970325, (2, 1) = 12484906049/39249798925, (2, 2) = 65745570806567/5769720441975})

(10)

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

CCnew6 := Matrix(2, 2, {(1, 1) = 20398649489879/4121228887125, (1, 2) = -5446073753219/242328258562950, (2, 1) = 14465459393/196248994625, (2, 2) = 8502096238511/1648491554850})

(11)

S*dmax

14

(12)

M1 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), Matrix(S, S*dmax, [seq(-CCnew || n, n = 0 .. dmax-1)])])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -336750255587/3769550688757, (1, 2) = 14853552191797/1696297809940650, (1, 3) = -49655436033349/56543260331355, (1, 4) = 75647656451147/1413581508283875, (1, 5) = -299962512141959/80776086187650, (1, 6) = 1231816081155781/8481489049703250, (1, 7) = -50445725001719/5769720441975, (1, 8) = 9065291388901/40388043093825, (1, 9) = -142685068141037/11539440883950, (1, 10) = 116560067351321/565432603313550, (1, 11) = -60560690824604/5769720441975, (1, 12) = 88774498025543/848148904970325, (1, 13) = -20398649489879/4121228887125, (1, 14) = 5446073753219/242328258562950, ...

 

M2 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity)])

M2 := Matrix(14, 14, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 1, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (1, 10) = 0, (1, 11) = 0, (1, 12) = 0, (1, 13) = 0, (1, 14) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 1, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (2, 10) = 0, (2, 11) = 0, (2, 12) = 0, (2, 13) = 0, (2, 14) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (3, 5) = 1, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (3, 10) = 0, (3, 11) = 0, (3, 12) = 0, (3, 13) = 0, (3, 14) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 0, (4, 5) = 0, (4, 6) = 1, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (4, 10) = 0, (4, 11) = 0, (4, 12) = 0, (4, 13) = 0, (4, 14) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (5, 7) = 1, (5, 8) = 0, (5, 9) = 0, (5, 10) = 0, (5, 11) = 0, (5, 12) = 0, (5, 13) = 0, (5, 14) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0, (6, 7) = 0, (6, 8) = 1, (6, 9) = 0, (6, 10) = 0, (6, 11) = 0, (6, 12) = 0, (6, 13) = 0, (6, 14) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 0, (7, 5) = 0, (7, 6) = 0, (7, 7) = 0, (7, 8) = 0, (7, 9) = 1, (7, 10) = 0, (7, 11) = 0, (7, 12) = 0, (7, 13) = 0, (7, 14) = 0, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 0, (8, 5) = 0, (8, 6) = 0, (8, 7) = 0, (8, 8) = 0, (8, 9) = 0, (8, 10) = 1, (8, 11) = 0, (8, 12) = 0, (8, 13) = 0, (8, 14) = 0, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = 0, (9, 6) = 0, (9, 7) = 0, (9, 8) = 0, (9, 9) = 0, (9, 10) = 0, (9, 11) = 1, (9, 12) = 0, (9, 13) = 0, (9, 14) = 0, (10, 1) = 0, (10, 2) = 0, (10, 3) = 0, (10, 4) = 0, (10, 5) = 0, (10, 6) = 0, (10, 7) = 0, (10, 8) = 0, (10, 9) = 0, (10, 10) = 0, (10, 11) = 0, (10, 12) = 1, (10, 13) = 0, (10, 14) = 0, (11, 1) = 0, (11, 2) = 0, (11, 3) = 0, (11, 4) = 0, (11, 5) = 0, (11, 6) = 0, (11, 7) = 0, (11, 8) = 0, (11, 9) = 0, (11, 10) = 0, (11, 11) = 0, (11, 12) = 0, (11, 13) = 1, (11, 14) = 0, (12, 1) = 0, (12, 2) = 0, (12, 3) = 0, (12, 4) = 0, (12, 5) = 0, (12, 6) = 0, (12, 7) = 0, (12, 8) = 0, (12, 9) = 0, (12, 10) = 0, (12, 11) = 0, (12, 12) = 0, (12, 13) = 0, (12, 14) = 1, (13, 1) = 0, (13, 2) = 0, (13, 3) = 0, (13, 4) = 0, (13, 5) = 0, (13, 6) = 0, (13, 7) = 0, (13, 8) = 0, (13, 9) = 0, (13, 10) = 0, (13, 11) = 0, (13, 12) = 0, (13, 13) = 0, (13, 14) = 0, (14, 1) = 0, (14, 2) = 0, (14, 3) = 0, (14, 4) = 0, (14, 5) = 0, (14, 6) = 0, (14, 7) = 0, (14, 8) = 0, (14, 9) = 0, (14, 10) = 0, (14, 11) = 0, (14, 12) = 0, (14, 13) = 0, (14, 14) = 0})

(13)

lowerband := Matrix(S, S*dmax, [seq(-evalf(CCnew || n), n = 0 .. dmax-1)])

lowerband := Matrix(2, 14, {(1, 1) = -0.893343e-1, (1, 2) = 0.875645e-2, (1, 3) = -.878185, (1, 4) = 0.535149e-1, (1, 5) = -3.71351, (1, 6) = .145236, (1, 7) = -8.74318, (1, 8) = .224455, (1, 9) = -12.36499, (1, 10) = .206143, (1, 11) = -10.49630, (1, 12) = .104669, (1, 13) = -4.94965, (1, 14) = 0.224740e-1, (2, 1) = -0.864552e-2, (2, 2) = -.115162, (2, 3) = -0.776011e-1, (2, 4) = -1.10141, (2, 5) = -.282396, (2, 6) = -4.51408, (2, 7) = -.537873, (2, 8) = -10.26209, (2, 9) = -.568850, (2, 10) = -13.97156, (2, 11) = -.318088, (2, 12) = -11.39493, (2, 13) = -0.737097e-1, (2, 14) = -5.15750})

(14)

M3 := Matrix(S*dmax, S*dmax, [Matrix(S*(dmax-1), S), Matrix(S*(dmax-1), shape = identity), lowerband])

Error, (in Matrix) this entry is too tall or too short: Matrix(2, 14, {(1, 1) = -0.8933432215e-1, (1, 2) = 0.8756453086e-2, (1, 3) = -.8781848755, (1, 4) = 0.5351488825e-1, (1, 5) = -3.713506389, (1, 6) = .1452358276, (1, 7) = -8.743183575, (1, 8) = .224454831, (1, 9) = -12.36498974, (1, 10) = .2061431666, (1, 11) = -10.49629552, (1, 12) = .1046685287, (1, 13) = -4.949652167, (1, 14) = 0.2247395242e-1, (2, 1) = -0.8645515381e-2, (2, 2) = -.1151624586, (2, 3) = -0.7760109089e-1, (2, 4) = -1.101411529, (2, 5) = -.2823958356, (2, 6) = -4.514076892, (2, 7) = -.537872817, (2, 8) = -10.26209174, (2, 9) = -.5688498906, (2, 10) = -13.97155838, (2, 1...

 

``


Download Test_Block_matrix.mw

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 7 Page 3 of 9