restart:

eq1:=diff(A1(t),t$2)+2*diff(A2(t),t$2)+3*diff(A3(t),t$2)-4*diff(A1(t),t)+diff(A1(t),t)*diff(A2(t),t)+3*diff(A3(t),t);

eq2:=2*diff(A1(t),t$2)+3*diff(A2(t),t$2)-3*diff(A1(t),t)+2*diff(A2(t),t)+3*(A1(t));

eq3:=2*diff(A2(t),t$2)+3*diff(A3(t),t$2)-4*diff(A2(t),t)+2*diff(A3(t),t)+3*A1(t)+A2(t)+A3(t);

dsolve({seq}(eq||i,i=1..3));

restart:

eq1:=a1+3*a2+5*a3;eq2:=6*a1-12*a2+20;eq3:=8*a1-53*a2+23*a3;

fsolve({seq(eq||i,i=1..3)},{a1,a2,a3},a1=-5..0)

restart;

 asal := proc(n)
local asal,m,s,i;
 if n<2 then
 m :=n/2+1;
 s:=m;
 print(m);
 else
 for i from 2 to n do
 if irem(n,i)=0 then
 print("asal değil");
 elif i = m then
 print("asaldir");
 fi;
 od;
fi;
 end proc;

asal(1);

asal(5);

restart:

m:=10;# number of n[i]'s ,to put in the function #

P:=proc (M::list) local i ,S,f ;
for i to m do
f:=x->1+2*(1+x)^n[i];
S[i]:=plot(f(x),x=0..5);od:
return plots:-display(seq(S[i],i=1..m)):
end proc;

#example#

for i to m do
n[i]:=i:od:

M:=[seq](n[i],i=1..m);

P(M);

restart:

ode:=diff(y(x),x$2)+y(x);

dsolve(ode);

eval(ode,%);

ode1:=diff(y(x), x)+2*tan(y(x))*tan(x)-1;

dsolve(ode1);

eval(ode1,%);

restart:

P:=proc ( n) ; expand((a+b)^n) ; end proc ;

P(14);

P(5);

restart:

L:=[seq](i^2,i=1..20);

for i to 20 do
FileTools[Text][WriteFloat]( "example.txt", L[i] );
writeline("example.txt","    "):
od:
FileTools[Text][Close]( "example.txt" ):

restart:

P:=proc(n)  local f1,f2;
f1 := expand((a1*x^2+b1*x+c1)^n);

f2 := a1*x^2+b1*x+c1;

solve([coeff(f1, x, 4) = coeff(f2, x, 4),

coeff(f1, x, 3) = coeff(f2, x, 3),

coeff(f1, x, 2) = coeff(f2, x, 2),

coeff(f1, x, 1) = coeff(f2, x, 1),

coeff(f1, x, 0) = coeff(f2, x, 0)], [a1,b1,c1,d1]);

end proc:

P(5);

restart:

int(1/1-exp(a*x)*erfc(a*x),x); # indefinite integral

int(1/1-exp(a*x)*erfc(a*x),x=0..10);

evalf(%);

eval(%,a=5);

restart:

with(Finance):

with(Statistics):

X:=RandomVariable(Binomial(nn, pp));

PDF(X, t);

restart:

a := Vector(5, i -> i+1 );

b := Vector(6, i -> i^2 );

A1:=plot([seq]([m,a[m]],m=1..numelems(a)),color=blue,gridlines=false):A2:=plot([seq]([m,b[m]],m=1..numelems(b)),color=yellow,gridlines=false):

plots:-display(A1,A2);