Координаты центра тяжести плоской фигуры. Maxima or SciLab

[grelik]

Определить центр тяжести полусегмента параболы $y^2 = ax$ , отсеченного прямыми
$x = a, y = 0, (a > 0, y > 0) .$

/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 13.04.2 ] */

/* [wxMaxima: input   start ] */
wxplot2d([y=sqrt(a*x),x=a,y=0],[x,-0.5,1.5],[y,-0.5,1.5])$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
s:integrate(sqrt(a*x)-(a),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
'integrate(sqrt(a*x)-(a),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
mx:integrate(x*(sqrt(a*x)-(a)),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
'integrate(x*(sqrt(a*x)-(a)),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
my:integrate((0.5*((sqrt(a*x))^2-a^2)),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
'integrate((0.5*((sqrt(a*x)-(a))^2-(a)^2)),x,0,1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
xc:mx/s;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
yc:my/s;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([[discrete,[[9/20,0.45]]],sqrt(a*x),a],[x,-1,2],[y,-1,2],[style,points, lines,lines])$
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$

Вообщем вроде нашёл формулу по которой считать всё подставил, но как то непонятно он константу "a" не принимает, как быть?

[grelik]

Неужели никто не подскажет?