a=1;
While[a<10,
Print["a= ",a];
a=a+1;
]

For[i=0,i<9,i+=1,
Print["i = ",i];
]

f[x_]:=Cos [x]+Sin[2x];
g[x_]:=x+1;
Plot[{f[x],g[x]},{x,-9,10}]

h[x_,y_]:=3x^2-2Sin[y];
g[x_,y_]:=2x^3-2Cos[y];

Plot3D[{h[x, y],g[x,y]},{x,-2,2},{y,-2,2}]

a={{1,1,1},{2,2,2},{3,3,3}}
MatrixForm[a]

(({{1, 1, 1}, {2, 2, 2}, {3, 3, 3}}))
b={({{1, 1, 1}, {2, 2, 3}, {4, 5, 6}})}
b

{{{1,1,1},{2,2,3},{4,5,6}}}


c=[({{1, 2, 3}, {2, 3, 4}, {5, 6, 7}})]
c

A=(({{1, 1, 1, 6}, {1, -1, 1, 2}, {2, 2, -1, 6}}))


MatrixForm[RowReduce[A]]

Demanda:=100;
Preco=5
S=2;
Cm=4;
CT=Demanda*Preco+Demanda*S/q+Cm*q/2;

CT[1]

{min,val}=Minimize[{CT,q≥0},q]

Plot[CT,{q,1,100},Epilog->{Red,PointSize[0.025],Point[{q/.val,min}]}]

c={-5,-2};
a=(({{-1, 0}, {0, -1}, {-1, -2}}));
b={-3,-4,-9};
sol=LinearProgramming[c,a,b]

constr=And @@MapThread[GreaterEqual,{a.{x,y},b}];
region=RegionPlot[constr,{x,0,10},{y,0,10},PlotStyle->{orange}];
objplot=DensityPlot[c.{x,y},{x,0,10},{y,0,10},PlotLegends->Automatic];
Show[objplot,region,Graphics[{PointSize->Large,Point[sol]}]]

f=-Cos[x]/(x^2+6);
Plot[f,{x,-10,10},PlotRange->All]