Ask Your Question
3

Is there any way to define an as-yet-unknown function?

asked 2012-03-23 14:11:38 +0100

jdpipe gravatar image

updated 2012-03-23 14:13:32 +0100

I'd like to know if there's a way of declaring functions in sage that are as-yet unknown. For example, let's say I have a function

p = R*T/v - a(T)/v/(v+b)

And I would like to be able to take a derivative like this

deriv(p,T)

and be given something back a partial derivative something like

R/v - diff(a(T),T)/v/(v+b)

However at present I can't seem to put an abstract function a(T) into my expression or find anything in the documentation that says how this is done.

As I recall there was a way to do this with wxMaxima, so maybe I just haven't found the trick in Sage yet.

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted
6

answered 2012-03-23 14:17:30 +0100

Shashank gravatar image

updated 2020-01-05 00:25:44 +0100

slelievre gravatar image

Answer from 2012. See further down for 2020 update.

You can define a function using the keyword function. The script below gives the answer you expect.

T = var('T')
a = function('a', T)
R = var('R')
v = var('v')
b = var('b')
p = R*T/v - a(T)/v/(v+b)
diff(p, T)

[2020 update] Nowadays, function no longer takes the variable as an argument.

Use a = function('a') instead of a = function('a', T).

Example:

sage: version()
'SageMath version 9.0, Release Date: 2020-01-01'
sage: b, v, R, T = SR.var('b v R T')
sage: a = function('a')
sage: p = R*T/v - a(T)/v/(v+b)
sage: diff(p, T)
R/v - diff(a(T), T)/((b + v)*v)
edit flag offensive delete link more

Comments

Great. seems to work fine. The output is a bit weird though, in typeset sage, the diff(a(T),T) is written as D [ 0 ] ( a ) ( T ), not particular intuitive. But it works! Thank you.

jdpipe gravatar imagejdpipe ( 2012-03-23 18:53:03 +0100 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2012-03-23 14:11:38 +0100

Seen: 2,155 times

Last updated: Jan 05 '20