как построить функцию (частную производную) в Maple с использованием 3D-графика

У меня есть следующая функция

u := proc (x, t) options operator, arrow; 
    50+sum((100*n*Pi*sin(n*Pi)+100*cos(n*Pi)-100)*cos((1/10)*n*Pi*x)*exp(-(1/100)*n^2*Pi^2*t)/(n^2*Pi^2), n = 1 .. 20) 
    end proc;

plot3d(u(x, t), x = 0 .. 10, t = 0 .. 20)

Я пытаюсь построить его в Maple, но не дает мне правильную трехмерную форму.


person Alqosh Hazkeal    schedule 15.02.2017    source источник
comment
Что вы подразумеваете под правильной формой?   -  person zhk    schedule 15.02.2017
comment
Производная по чему?   -  person zhk    schedule 15.02.2017
comment
Я исправил ваш код (отсутствуют операторы умножения и исправил pi на Pi). Он все еще не работает?   -  person Therkel    schedule 15.02.2017


Ответы (1)


restart:
u := proc (x, t) options operator, arrow; 
    50+sum((100*n*Pi*sin(n*Pi)+100*cos(n*Pi)-100)*cos((1/10)*n*Pi*x)*exp(-(1/100)*n^2*Pi^2*t)/(n^2*Pi^2), n = 1 .. 20) 
    end proc:

plot3d(u(x, t), x = 0 .. 10, t = 0 .. 20);

введите здесь описание изображения

plot3d(diff(u(x, t),x), x = 0 .. 10, t = 0 .. 20);

введите здесь описание изображения

person zhk    schedule 15.02.2017
comment
plot3d не является частью пакета plots, поэтому with(plots): здесь избыточен. @alqoshhazkeal не ответил на мой комментарий к вопросу, но я не уверен, что ошибка сохраняется после моего редактирования, т.е. проблема была из-за небольших ошибок в коде. - person Therkel; 16.02.2017
comment
@Therkel Какая ошибка? После того, как вы внесли изменения, я больше не вижу никаких синтаксических ошибок. Вы правы насчет plot3d. - person zhk; 16.02.2017
comment
См. первую версию. Этот код не запускается. - person Therkel; 16.02.2017
comment
@Therkel Вот что я говорю после того, как вы внесли исправление, больше нет проблем с синтаксисом. - person zhk; 16.02.2017
comment
Я полагаю, моя точка зрения заключается в том, что нет никакой разницы между вопросом и вашим ответом. :) - person Therkel; 16.02.2017
comment
@Therkel В заголовке написано, как построить partial diff. Кстати, вы не должны исправлять синтаксис кода, редактируя вопрос. Посетителям это не поможет. - person zhk; 16.02.2017