Функции в Maple

HaVk · 06.12.2009

Функции в Maple

Такая проблема, определяю нужные функции:

restart:
with(student):

# Задаем закон движения ведущего звена в виде функции времени
omega1 := t -> Pi * n / 30;

#Вводим временные функции углов
alpha := t -> alpha0 + int(omega1(t),t=0..t);
beta := t -> Pi / 2 + alpha(t, alpha0);

#Функция времени l(t)
l := t -> -(a*(sin(alpha(t))*cos(beta(t))-cos(alpha(t))*sin(beta(t))))/(sin(beta(t)));

#Cкорость
VBa := t -> diff(l(t),t);

Задаю конкретные значения:

a := 100:
n := 60:
alpha0 := 0:
alpha1 := 20 * Pi * 1 / 180:
t0 := 0:
t1 := alpha1 / ((1 / 30) * Pi * n):

Затем пишу:

plot(VBa(t), t=t0..t1);

и получаю график, НО на команду:

VBa(0.04));

получаю:

Error, (in int) wrong number (or type) of arguments

Т.е. не могу узнать значение функции в конкретной точке

Кто-нибудь знает в чем проблема?

Burnash · 06.12.2009

Я не специалист, но вот это

omega1 := t -> Pi * n / 30

выглядит как очень странное определение функции. Может вместо n нужно t ?

HaVk · 06.12.2009

Нет, здесь всё правильно, просто omega1 не зависит от времени.

К тому же график имеет ожидаемую форму, т.е. все записи до

plot(VBa(t), t=t0..t1);

были правильными.

Burnash · 06.12.2009

А какого типа должен быть аргумент t во всех определениях? Он точно может быть не целым (ну Вы поняли)?

HaVk · 06.12.2009

Код:

t0 := 0:
t1 := alpha1 / ((1 / 30) * Pi * n):

t меняется от 0(нуля) до 0.05555555556
т.е. t - не целое.

Burnash · 06.12.2009

VBa(0.04));

Глюк со скобками??

HaVk · 06.12.2009

нет, это я здесь опечатался...
в Maple написал правильно:

Код:

VBa(0.04);

Вот не могу понять... plot же как-то вычисляет значение этой функции, чтоб нарисовать график. При этом t меняется от 0(нуля) до 0.05555555556.
0.04 должно работать...

Burnash · 06.12.2009

Может VBa может принимать как параметр только переменную, но не константу?
Пробуй:
q:=0.04
VBa(q)
Что-то в этом роде.

HaVk · 06.12.2009

Попробовал. Результат тот же

Нашлось решение!!!
Нужно было вместо

Код:

VBa(0.04);

писАть:

Код:

subs(t=0.04,VBa(t));