Оптимизация кода

[Falex]

Извините за такой вопрос.Изучаю Maple с сегодняшнего дня и вот какая проблема возникла:оптимизация кода.
Все дело в том,что код получается слишком горомоздким,хотя в других математических пакетах это делается достаточно
легко.Посоветуйте как можно избавиться от такой проблемы.
Вот файлик:

[aar]

Например:
> f := proc (zf) local j;
> add(abs((eval(find_cos, z = eval(zf, t = (j+1)/n)))-(eval(find_cos, z = eval(zf, t = j/n)))), j = 0 .. n)
> end proc:
> f(zf1)+f(zf2)+f(zf3)+f(zf4);
14.13556662

[Falex]

А массив функций нельзя сделать:чтобы не писать zf1,zf2,zf3,zf4 а например zf[1],...ну типа как-нибудь так?

[aar]

Можно массив, а можно и так:
add(f(zf||i), i=1..4);
14.13556662

[Falex]

Вот еще вопросик.написал код,который рисует 4 отрезка на различных plot'aх.А мне надо на одном.Как этого добиться?

Код: Выделить всё

line:=proc(t,j) if (j=1) then 1-2*t-1*I elif (j=2) then 1+I*(-1+2*t) elif (j=3) then 1-2*t+1*I elif (j=4) then -1+I*(1-2*t) end if end proc;
for j from 1 to 4 do plot([Re(line(t,j)),Im(line(t,j)),t=0..1]): end do;


[Irmantas]

line:=proc(t,j) if (j=1) then 1-2*t-1*I elif (j=2) then 1+I*(-1+2*t) elif (j=3) then 1-2*t+1*I elif (j=4) then -1+I*(1-2*t) end if end proc;
plots[display](seq(plot([Re(line(t,j)),Im(line(t,j)),t=0..1]),j=1..4));

[Falex]

Спасибо.А почему здесь выдает ошибку,хотя делаю я всё правильно:

Код: Выделить всё

> line:=proc(t,j) if (j=1) then 1-2*t-1*I elif (j=2) then 1+I*(-1+2*t) elif (j=3) then 1-2*t+1*I elif (j=4) then -1+I*(1-2*t) end if end proc;
> z1:=-0.999-0.999*I;z2:=-z1;n=1000;

> find_cos:=proc(z) arccos(((Re(z)-Re(z1))*(Re(z)-Re(z1))+(Im(z)-Im(z1))*(Im(z)-Im(z1))+(Re(z)-Re(z2))*(Re(z)-Re(z2))+(Im(z)-Im(z2))*(Im(z)-Im(z2))-(Re(z1)-Re(z2))*(Re(z1)-Re(z2))-(Im(z1)-Im(z2))*(Im(z1)-Im(z2)))/(2*sqrt((Re(z)-Re(z1))*(Re(z)-Re(z1))+(Im(z)-Im(z1))*(Im(z)-Im(z1)))*sqrt((Re(z)-Re(z2))*(Re(z)-Re(z2))+(Im(z)-Im(z2))*(Im(z)-Im(z2))))) end proc;
> variation:=proc(j)::float; local summ::float,i::integer;summ:=0;
>   for i from 0 to n do summ:=summ+find_cos(line(i/n,j)) end do;
>   return summ;
> end proc;
> print(variation(3));
Error, (in variation) final value in for loop must be numeric or character


[Irmantas]

рестартуй весь документ.
проедуру find_cos можно написать и папроще
find_cos:=z-> arccos(( abs(z-z1)^2+abs(z-z2)^2-abs(z1-z2)^2 )/(2*abs(z-z1)*abs(z-z2)));

[Falex]

неа.Все равно эту ошибку выдает.Вот файлик:

[Irmantas]

неа.Все равно эту ошибку выдает.Вот файлик:

забыли присваивание n.
n:=2000;

[Falex]

Вот блин:всего напрос-то точку с запятой забыл
А почему ответ то не пишет,а такую строчку (при n=6):

Код: Выделить всё

evalf(variation(1));
find_cos(1 - I) + find_cos(2/3 - I) + find_cos(1/3 - I) + find_cos(-I) + find_cos(-1/3 - I) + find_cos(-2/3 - I) + find_cos(-1 - I)


[Irmantas]

А почему ответ то не пишет,а такую строчку (при n=6):

Код: Выделить всё

evalf(variation(1));
find_cos(1 - I) + find_cos(2/3 - I) + find_cos(1/3 - I) + find_cos(-I) + find_cos(-1/3 - I) + find_cos(-2/3 - I) + find_cos(-1 - I)


при n:=6
variation(1) выводит резултат 11.64999743
сделайте рестарт

[Falex]

А как сделать рестарт?!
Просто выход из программы и запуск заново не поможет?

[Irmantas]

А как сделать рестарт?!
Просто выход из программы и запуск заново не поможет?

просто напишите в начале документа restart: и нажмите кнопку с "!!!"

[Falex]

Все равно не помогло (