File: desoln.mac

package info (click to toggle)
maxima 5.47.0-9
  • links: PTS
  • area: main
  • in suites: forky, sid
  • size: 193,104 kB
  • sloc: lisp: 434,678; fortran: 14,665; tcl: 10,990; sh: 4,577; makefile: 2,763; ansic: 447; java: 328; python: 262; perl: 201; xml: 60; awk: 28; sed: 15; javascript: 2
file content (21 lines) | stat: -rw-r--r-- 817 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
/* -*- Mode: MACSYMA; Package: CL-MAXIMA -*- */

load(pwilt);

desolve(eqns,vars):=block([teqns,tvars,ovar,lvar,flag,dispflag],
	modedeclare(flag,boolean),
	flag:false,
	if not listp(vars) then (eqns:[eqns], vars:[vars], flag:true),
	if length(ovar:vars[1])#1 then
		error(ovar,"contains more than one independent variable."),
	ovar:inpart(ovar,1),dispflag:false,
	teqns:map(lambda([z],laplace(z,ovar,lvar)),eqns),
	tvars:map(lambda([z],'laplace(z,ovar,lvar)),vars),
	/* if not freeof(nounify('diff),teqns) then return(false), */
	teqns:errcatch(solve(teqns,tvars)),
	if teqns=[] or teqns=[[]] 
	  then error("desolve can't handle this case.")
	  else teqns:first(teqns),
	if flag#true then teqns:first(teqns),
	teqns:map(lambda([z],pwilt(z,lvar,ovar)),teqns),
	if flag and length(tvars)=1 then teqns[1] else teqns)$