torabi

70 Reputation

One Badge

3 years, 210 days

MaplePrimes Activity


These are questions asked by torabi

Hello,

How I can remove some problems in my code.

I wrote it before in  software, but in some lines, I have some problems with it.

How I can change it in maple to run well?


 

restart

M := 2; B := 1; C := 3; K := 6; f := 2*exp(-t); q0 := 1; qd0 := -1; qd20 := 1; alpha := 1/2; beta := 1/12; gammaa := 1/4; h := .1; Tf := 2

"N:=size(M,1)  T:=0:h:Tf;  lt:=length(T);  F:=subs(f,t,T);  F:=subs(t=T,f);    qd30:=M\(F(1)-B*qd20-C*qd0-K*q0);   q:=zeros(N,lt)  qd:=zeros(N,lt)  qd2:=zeros(N,lt);  qd3:=zeros(N,lt);  q(:,1):=q0  qd(:,1):=qd0  qd2(:,1):=qd20  qd3(:,1):=qd30"

Error, invalid function arguments

"N:=size(M,1)  T:=0:h:Tf;  lt:=length(T);  F:=subs(f,t,T);  F:=subs(t=T,f);    qd30:=M\(F(1)-B*qd20-C*qd0-K*q0);  q:=zeros(N,lt)  qd:=zeros(N,lt)  qd2:=zeros(N,lt);  qd3:=zeros(N,lt);  q(:,1):=q0  qd(:,1):=qd0  qd2(:,1):=qd20  qd3(:,1):=qd30"

 

"A:=M+alpha*h*B+gammaa*(h^2)*C+beta*(h^3)*K;  for n:=1:lt-1 do      D:=F(n+1)-B*(qd2(:,n)+(1-alpha)*h*qd3(:,n))-...          C*(qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n))-...          K*(q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n));      qd3(:,n+1):=A\D;      qd2(:,n+1):=qd2(:,n)+(1-alpha)*h*qd3(:,n)+alpha*h*qd3(:,n+1);      qd(:,n+1):=qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n)+gammaa*(h^2)*qd3(:,n+1);      q(:,n+1):=q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n)+...          beta*(h^3)*qd3(:,n+1);  end do"

Error, invalid function arguments

"A:=M+alpha*h*B+gammaa*(h^2)*C+beta*(h^3)*K;  for n:=1:lt-1 do      D:=F(n+1)-B*(qd2(:,n)+(1-alpha)*h*qd3(:,n))-...          C*(qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n))-...          K*(q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n));      qd3(:,n+1):=A\D;      qd2(:,n+1):=qd2(:,n)+(1-alpha)*h*qd3(:,n)+alpha*h*qd3(:,n+1);      qd(:,n+1):=qd(:,n)+h*qd2(:,n)+(0.5-gammaa)*(h^2)*qd3(:,n)+gammaa*(h^2)*qd3(:,n+1);      q(:,n+1):=q(:,n)+h*qd(:,n)+(h^2)/2*qd2(:,n)+(1/6-beta)*(h^3)*qd3(:,n)+...          beta*(h^3)*qd3(:,n+1);  end do"

 

plot(T, q, 'r-o', 'linewidth', 1.5)

Error, (in plot) unexpected options: [r-o, , 1.5]

 

``


 

Download

 

Please help me for removing an error in my code.

Thank you


 

"   restart:     NL:=960:    LL:=4:      BB:=0.6 :  T:=16:    MM:=NL*(LL+1):          for  NC  from 0 by 1 to MM-1  do       S:=BB+(2*Pi*NC/T)*I;          powerFactor := 0.05 * (sqrt(2.0) - 1.202081528);     powerTerm := sqrt(S)*sqrt(49*S^(2)+280*S+800)/sqrt(S+4);       preFactor := 1/((exp(sqrt(2.0)*powerTerm/10.0)-1.0)*S);       XS(NC+1):= preFactor * ((-exp(-powerFactor*powerTerm) + exp(powerFactor*powerTerm)) * exp(sqrt(2.0)*powerTerm/20.0));    end  do;    for  timeCounter from 0 by 1 to  NL-1  do     totalTerm:=0;     for KK from 0 to NL-1 do        for L from 0 to LL do           Term1:=Re(XS(KK+L*NL+1))+I*(Im(XS(KK+L*NL+1)));           Term2:=exp((timeCounter*KK*2*Pi/NL)*I);           totalTerm:=totalTerm+Term1*Term2;        end do     end do     XT(timeCounter+1):=((2.0/T)*exp(BB*timeCounter*T/NL)*(-0.5*Re(XS(1))+Re(totalTerm)));     TT(timeCounter+1):=(timeCounter*T)/NL;         end do    "

Error, invalid loop statement termination

"   restart:    NL:=960:  LL:=4:  BB:=0.6 :  T:=16:    MM:=NL*(LL+1):       for NC  from 0 by 1 to MM-1  do     S:=BB+(2*Pi*NC/T)*I;    powerFactor := 0.05 * (sqrt(2.0) - 1.202081528);   powerTerm := sqrt(S)*sqrt(49*S^2+280*S+800)/sqrt(S+4);    preFactor := 1/((exp(sqrt(2.0)*powerTerm/10.0)-1.0)*S);     XS(NC+1):= preFactor * ((-exp(-powerFactor*powerTerm) + exp(powerFactor*powerTerm)) * exp(sqrt(2.0)*powerTerm/20.0));    end  do;    for timeCounter from 0 by 1 to NL-1  do  totalTerm:=0;     for KK from 0 to NL-1 do        for L from 0 to LL do  Term1:=Re(XS(KK+L*NL+1))+I*(Im(XS(KK+L*NL+1)));  Term2:=exp((timeCounter*KK*2*Pi/NL)*I);  totalTerm:=totalTerm+Term1*Term2;        end do     end do   XT(timeCounter+1):=((2.0/T)*exp(BB*timeCounter*T/NL)*(-0.5*Re(XS(1))+Re(totalTerm)));   TT(timeCounter+1):=(timeCounter*T)/NL;      end do    "

 

``

``


 

Download

 

 

how I can write a program code for newmark method.

in this method time has 3 order derivation

I can  this equation.

CV.mw
 

restart; c__v := 1.2; `τ__q` := 8.5*10^(-12); `τ__T` := 90.0*10^(-12); rho := 1000; k := 10

1.2

 

0.8500000000e-11

 

0.9000000000e-10

(1)

k*(diff(T(x, t), x, x))+k*`τ__T`*(diff(T(x, t), t, x, x)) = rho*c__v*(diff(T(x, t), t))+(diff(T(x, t), t, t))*c__v*rho*`τ__q`+(1/2)*c__v*rho*`τ__q`^2*(diff(T(x, t), t, t, t))

10*(diff(diff(T(x, t), x), x))+0.9000000000e-9*(diff(diff(diff(T(x, t), t), x), x)) = 1200.0*(diff(T(x, t), t))+0.1020000000e-7*(diff(diff(T(x, t), t), t))+0.4335000000e-19*(diff(diff(diff(T(x, t), t), t), t))

(2)

Boundary condition:

T(0, t) = 300; T(10, t) = 300

#####################################

INITIAL CONDITIONS:

 

T(x, 0) = 300; (D[1](T))(x, 0) = 0, (D[2](T))(x, 0) = 0

(D[1](T))(x, 0) = 0, (D[2](T))(x, 0) = 0

(3)

``


 

Download CV.mw

 

 

How I can take Laplace Transform from equation.

Thanks

LAPLACE

1 2 3 4 5 6 7 Last Page 1 of 18